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ABSTRACT 

A.l   Introduction 

It  is  our  purpose  to  establish  a  theory  describing  the  state  transition 
of  transistor  flipflops  making  use  mostly  of  phase  plane  techniques,  but  using 
also  a  time  or  even  frequency  domain  point  of  view  whenever  helpful.   Based  on 
such  a  theory,  we  further  wish  to  devise  practical  engineering  methods  of 
analysis,  design,  and  optimization  of  transistor  flipflops. 

We  shall  restrict  ourselves  to  considering  the  asymmetric  (Fig.  A.l) 
and  the  Eccles-Jordan  (Fig.  A. 2)  flipflops,  with  constant  current  T  fed  into 
the  common  emitters,  and  except  for  a  short  discussion,  we  shall  consider  only 
the  case  of  a  rectangular  trigger. 

A°2  The  Flipflop  Differential  Equation 

Based  on  the  differential  equations  relating  terminal  voltages  with 
the  charges  at  the  junctions  and  diffusion  tails  in  a  transistor  we  obtain  the 
transistor  pair'  characteristics  below. 

Wk  =  2[1  "  (-!)ktanh  XJ  (A. la) 


o    l  -  a 
Zk=Wk  +  ~^~wk  ^. lb) 


where  k  -  1,  2,  is  the  transistor  index  used  in  Figs.  1  and  2,  and  the  symbol 

o  _  dx 

x  -  t-  ;  besides, 

x  =  x1   ~   *2  =  normalized  base-to-base  voltage 


Throughout  this  volume  the  expression  "transistor  pair"  refers  to  the  pair 
of  identical  transistors  with  the  constant  current  I_  fed  into  the  common 
emitter,  as  in  Figs.  1  and  2.  h 
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qvk  t 

Xk  ~  2kT  ~  normalized  base  voltages  ' 


XCk 
Wk  ~  al~  "  normalized  collector  current 
E 


XBk 
k  _  ^Y~  ~  normalized  base  current 


S 


t  -  —  -  normalized  time  variable 


t'  =  real  time  variable 
t  =  collector  time  constant  of  the  transistors 

Equations  (A.l)  therefore,  relate  the  collector  and  base  currents  of 
the  transistor  pair  to  the  base-to-base  voltage. 

Analysis  of  the  feedback  networks  yields  a  differential  equation  in 
the  variable  xfc  for  each  value  of  k  (k  =  1,  2  for  the  Eccles -Jordan,  but  k  =  1 
for  the  asymmetric  flipflop),  respectively  related  to  trigger  plus  base  currents 
^k  +  ^k^  and  collector  currents  i  „}    £  f   k. 

We  shall  concentrate  our  attention  on  the  asymmetric  flipflop.   The 
result,  after  normalization,  is: 

T  T       T   +  T    +  T 

i  o  oo    ioioo  f        T-o         R 

—  -  x  + x  +  x  =  2P  Jb  +  w2  +  -i°  (ex  +  §x)  +  _£  (Qi  +   Zi)j 

(A. 2) 
where : 


f 

The  k  m  the  denominator  is  the  Boltzmann  constant. 
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T.    =   R.C. 

l  11 


T.       -    R.C 
lO  1     O 


T    .    -  R  C. 

Ol  O    1 


A 


y    --  time   constants 


T       =    R    C 
o  o  o 

J 


IJ     +  E 

■D      -  B     S  C 

aLR ~  biasinS  condition 

"TE  o 


\  ~   £~T~  -  normalized  trigger  current  actually 
T!   fed  into  the  flipflop  circuit 


Since  G±   is  the  trigger  current  actually  fed  into  the  flipflop,  it 
depends  on  the  nature  of  the  trigger  circuit,  and  depends  also  on  x  itself. 

In  general,  if  the  trigger  circuit  of  transistor  T   is  represented  by 
a  current  source  of  intensity  i   and  shunt  interval  conductances  G  ,  we  get 


K     =   s  -   G,  x 
k    k   k 


(A.3) 


where, 


-  tk         -,  •   a 
Sk  "  ccT     ~   nomallzed  trigger  current  source  intensity. 
E 


Use  of  (A.l)  with  (A. 2)  and  (A.3)  will  result  in  a  second  order  non- 
linear differential  equation  in  x  containing  terms  in  tanh  x  along  with  its 
first  and  second  derivatives,  and  under  a  forcing  function  s  (t),  and  its  first 
derivative  s^t).   This  equation  shall  be  called  "the  flipflop  equation." 
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A.3  Piecevise  Linear  Approximation 

A  general  solution  to  the  flipflop  equation  is  not  known.   Let  us 
consider,  however,  instead  of  (A.l),  the  characteristics  below: 


^k  =  |[1  -  (-D%(x)] 


o    1  -  a     /        . 
wk  +  — ^ —  wt,  (as  before) 


a 


(A. 4a) 

(AAb) 


where, 


r 

-1, 


<P(x)  =  <7x, 
+1; 


x  <  -  -  :   region  I 


|xj  <   -  :   region  II 


x     -  +  y    :   region  III 


(A. 5) 


and  7  is  selected  so  that  cp(x)  is  the  best  possible  piecewise  linear  approxi- 
mation to  tanh  xJ  We  are,  therefore,  considering  an  ideal  transistor  pair 
whose  behavior  approximates  that  of  the  real  transistor  pair,  but  whose  char- 
acteristics are  piecewise  linear  functions. 

Equations  (A. 2),  (A. 3)  and  (A.k)   will  furnish  equation  (A. 6),  i.e., 
a  second  order  differential  equation  in  x,  where  the  forcing  functions  s  (t) 
and  sx(t),  although  nonlinear  in  the  whole  range  of  x,  will  be  formed  by  three 
linear  equations,  respectively  valid  in  the  three  regions  as  in  equation  (A. 5) 
above.   The  second  derivative  of  cp(x)  is  a  pair  of  impulse  functions  occurring 
in  the  transitions  of  x  into  and  out  of  region  II. 


1 


With  e(x)  -  cp(x)  -  tanh  x,  both  a  criterion  of  getting  7  such  that  max  e(x) 

r»  ,  ,  x 

is  minimum  or  7  such  that  /   e(x)dx  =  0  will  lead  to  7  «  0.7. 

J0 
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00    -,  O  O  n  1 

ax  +  bx  +  ex  =  d  +  ms  +  ns  +  f  (x)  •  s(x  +  -)         (a. 6) 

where  a,  b,  c,  d,  m  depend  on  the  region  containing  x,  i.e.,  they  are  step 

functions  of  x,  discontinuous  at  x  =  +  - 

7 

A  solution  of  (A. 6)  can  be  obtained  by  solving  a  second  order  linear 
differential  equation  in  each  region,  using  the  final  conditions  in  each  region 
as  the  initial  conditions  of  the  next  region.   The  impulse  functions  can  be 
considered  as  discontinuities  in  the  value  of  x. 

A°^  The  Rectangular  Trigger  and  the  Phase  Plane  Portrait 

If  s1(t)  is  a  rectangular  function  of  time,  of  duration  T,  then  s  (t) 

ov  ' 

is  also  a  pair  of  impulse  functions  and  can  be  taken  account  of  as  discontinu- 
ities of  S,  at  t  =  0  and  t  =  T,  whereas  s^t)  itself  can  be  represented  as  a 
constant  equal  to  W  if  the  trigger  is  ON  and  0  (zero)  if  the  trigger  is  OFF. 

We  have,  so  far,  succeeded  in  reducing  the  analysis  of  a  flipflop  to 
the  solution  of  three  equations  like  (A. 7)  below: 


ax  +  bx  +  ex  =  D  (A. 7) 


d  +  W   if  the  trigger  is  ON 
D  = 


d     if  the  trigger  is  OFF 

in  three  regions  of  values  of  x,  matching  their  solutions  in  the  boundaries 
+  -  ,  and  considering  discontinuities  of  x  at  t  =  0,  t  =  T,  and  at 


x  =  7± 

7 


The  phase  plane  equation  corresponding  to  (A. 7)  is 


^Z  -  D  "  GX   ^ 

dx     ay   "  a  (A.o> 
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where 

o 

y  =  x 

The  singular  points  (points  of  y  =  y  =  o)  are  given  by  the  solutions  of 

D  -  ex  -  0  (A.9) 

i.e., 


D 
c 


x  =  "  (A. 10) 


and  the  nature  of  these  singular  points  can  be  analyzed  by  studying  the  natural 
frequencies  of  the  system  (A.ll)  below  in  a  neighborhood  around  each 
singularity. 


o 

x  =  y 


ode     b 
y  =  — x y 

a   a     a  J 

J 


r  (A.ii) 


This  analysis  shows  that,  if  the  trigger  is  OFF,  and  under  proper 
biasing  conditions,  these  singularities  are  two  stable  nodes, ^  one  in  each 
region  I  and  III,  and  a  saddle  point'"''  in  region  II. 

The  action  of  a  trigger  of  amplitude  W  is  essentially  to  shift  the 

stable  nodes  by  an  amount  Ax  =  ^  ,  and  the  saddle  point  by  an  amount  Ax  =  -  22 

u  c 

(parameters  always  calculated  in  the  correct  regions!). 


In  a  neighborhood  of  a  stable  node  both  natural  frequencies  of  the  system 
are  real  and  negative  (by  definition). 

In  a  neighborhood  of  a  saddle  point  the  natural  frequencies  of  the  system 
are  real  and  have  opposite  signs  (by  definition). 
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When  a  singular  point  is  shifted  out  of  its  proper  region  we  say  it 
has  become  "virtual"  because  its  nature  still  determines  the  behavior  of  the 
system  in  that  region,  but  it  does  not  really  exist.   There  is  a  value  of  W 
above  which  the  saddle  point  and  one  of  the  stable  nodes  become  virtual.   Then 
the  representative  point  P  of  the  system  will  describe  a  trajectory  towards 
the  remaining  stable  node.   When  the  trigger  is  turned  OFF  this  remaining  stable 
node  (and  also  the  other  two  singularities)  returns  to  its  resting  position; 
P  will  now  move  towards  this  point,  which  is,  in  effect,  one  of  the  two  stable 
states  of  the  flipflop.   Figures  A. 3a,  b,  c,  d  illustrate  a  transition. 

Clearly,  improper  biasing  may  result  in  other  singularity  configura- 
tions which  will  not  correspond  to  flipflop  behavior. 

Study  of  some  geometrical  properties  of  the  flipflop  phase  plane 
portrait,  expecially  the  study  of  the  separatrices  and  of  certain  properties  of 
the  trajectories  can  be  made  by  analytical  solutions  to  equation  (A, 8).   The 
usefulness  of  this  approach  is  limited,  however,  by  the  complexity  of  the 
algebraic  expressions  involved. 

A°5  Engineering  Methods  for  the  Solution  of  Flipflop  Problems 

These  methods  consist  essentially  of  the  drawing  of  figures  like  A. 3, 
i.e.,  phase  plane  portraits  which  are  good  approximations  to  the  true  phase 
plane  portrait  of  a  trajectory. 

The  time  durations  of  any  portion  of  a  trajectory,  say,  between  two 
points  Pa  and  P  is  given  by: 

[^   i 


X 

a 


But  if  y(x)  is  a  straight  line, 
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Various  procedures  can  be  devised  for  the  drawing  of  these  approxi- 
mate portraits  once  the  coefficients  of  (A.8)  are  known  in  every  region,  both 
for  trigger  ON  and  OFF.  From  any  approximate  form  of  y(x)  it  is  a  simple  matter 
to  obtain  any  waveform  of  interest,  such  as  x(t),  wjt),  z^t),  etc.,  by  very 
simple  graphical  constructions. 

Even  though  only  the  asymmetrical  flipflop  is  considered  in  this 
Abstract,  it  is  very  simple  to  show  that  the  symmetric  Eccles-Jordan  flipflop 
is  formally  equivalent  to  the  asymmetric  circuit . 

A.  6  Extensions 

A.6.1--Simple  formulae  for  fast  estimates  of  flipflop  transition  times  can  be 
derived  from  geometrical  properties  of  the  transition  phase  plane  portraits. 

A.6.2--The  general  (nonsymmetric )  Eccles-Jordan  flipflop  can  be  treated  by  a 

modified  phase  plane  which  takes  into  account  the  two  sets  of  variables  (x     v  ) 

v  k^k"'' 

■■   1,  2,  with  some  use  of  curves  in  the  time  domain  x,  (t),  and  of  a  plane 
(xi,x2).   In  some  cases  where  there  is  a  certain  relationship  between  its  para- 
meters a  nonsymmetric  flipflop  may  be  reduced  to  the  asymmetric  case. 

A.6.3-Some  degenerate  cases,  where  one  or  more  capacitances  are  zero  can  have 
y(x)  expressed  in  analytical  and  closed  form. 

A.6^-Flipflops  under  other  types  of  trigger  can  be  studied  by  other  techniques 
than  phase  plane  techniques,  such  as  frequency  or  time  domain  techniques. 
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Ao7  Practical  Results  and  Future  Developments 

Comparison  of  theoretical  with  experimental  results  indicates  that 
an  error  of  the  order  of  +  10  per  cent  is  to  be  expected  from  the  application 
of  this  theory. 

There  are  strong  indications  that  a  large  part  of  this  error  is 
caused  by  the  fact  that,  although  cp(x)  s  tanh  x,  this  is  not  so  for  the  first 
and  second  derivatives  of  these  functions.   However,  the  results  are  satisfactory 
and  the  theory  could  be  extended  in  future  to  consider,  for  example:   the  non- 
linearity  of  the  capacitances  involved,  inductances  in  the  passive  network,  or 
perhaps  the  fact  that  these  are  distributed  parameters.   In  another  direction, 
it  could  be  improved  by  a  formal  mathematical  investigation  of  the  piecewise 
linear  approximation  method  for  solving  nonlinear  differential  equations. 


1 .   INTRODUCTION 


Many  attempts  have  been  made  in  the  past  to  establish  a  theory  of 
bistable  systems .   The  nonlinear  nature  of  bistable  circuits  (f lipf lops ),  how- 
ever, makes  it  impractical  to  approach  the  analysis  of  such  systems  by  means 
of  mathematical  techniques  developed  specially  for  the  study  of  linear  systems. 

Nevertheless,  several  flipflop  theories  have  been  devised  making  use 
of  strictly  linear  techniques  [3,  16,  19,  23]. ^   However,  these  theories,  one 
way  or  another,  had  limited  objectives.   Typically,  each  of  them  aims  at  a 
specific  problem  among,  for  example,  obtaining  an  approximation  for  the  transi- 
tion time,  or  finding  minimum  values  for  trigger  duration  and  amplitude,  or 
studying  some  aspect  of  stability,  etc.   These  limited  objectives  could  be, 
and  were,  attained  with  the  linear  techniques  employed. 

Adopting  a  quite  different  point  of  view,  flipflop  theory  can  be 
reduced  to  the  study  of  nonlinear  bistable  equations  or  systems  of  equations  (in 
general  two  first  order  equations)  [26,  27],   The  phase  plane  is  the  mathematical 
tool  that  immediately  suggests  itself  for  this  type  of  problem  [l,  26]. 

In  the  literature,  authors  are  usually  not  concerned  with  the  analysis 
of  any  specific  physical  problem,  or  mathematical  model,  but  just  with  the 
development  of  the  mathematical  technique  itself.   Objectives  here  were  also 
limited,  since  the  establishment  of  a  general  theory  was  not  aimed  at,  but  it 
was  desired  to  show  that  some  results  could  be  obtained  by  applying  phase-plane 
methods  to  certain  types  of  equations  [26,  27]. 

It  seems  natural  to  attempt  a  general  analysis  of  bistable  physical 
systems  using  the  phase-plane  as  the  main  technique.   However,  a  general 
analysis  of  this  kind  would  be  of  very  little  use  (if  it  could  be  made  at  all) 


A  square-bracketed  set  of  numbers  refers  to  works  listed  under  "Bibliography' 
at  the  end  of  this  dissertation. 
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first  because  nonlinear  systems  do  not  lend  themselves  too  easily  to 
generalization  and  second,  because  whatever  generalization  we  would  achieve, 
we  would  pay  by  not  being  specific  about  the  most  important  of  all  electric 
bistable  systems  (we  consider  the  common  switch  as  being  a  mechanical  bistable 
system...)  and  the  only  one  in  which  we  are  really  interested,  namely,  the 

transistor  flipflop. 

Thus  the  analysis  of  transistor  flipf lops— and  of  these  only- -appears 
to  be  attractive,  since  specific  results  and  methods  can  be  obtained  and  applied 
immediately.   This  is  the  objective  of  this  report:   a  general  analysis  of 
transistor  flipf lops  by  phase-plane  methods.   It  is  general  in  the  sense  that 
any  information  about  the  flipflop  behavior  can  be  obtained  from  it,  and  also 
in  the  sense  that  the  general  (nonsymmetric )  Eccles-Jordan  flipflop  is  con- 
sidered.  Its  differential  equations  are  established  along  with  the  more  impor- 
tant cases  of  the  asymmetric  (one  base  grounded)  and  the  symmetric  flipf lops, 
and  we  go  as  far  as  suggesting  means  to  analyze  this  general  case  on  a  special 

phase -plane. 

Of  course,  phase-plane  treatment  practically  outlaws  any  other  but 
the  rectangular  trigger  (as  will  be  clear  in  the  sequel).   However,  the  equations 
established  do  not  demand  phase-plane  treatment;  in  fact,  for  other  types  of 
trigger,  numeric  solutions  may  be  necessary.   But  we  shall  concentrate  our 

attention  on  the  rectangular  trigger  case.   The  first  part  of  the  theory  is 

4 
concerned  with  the  establishment  of  the  general1  dynamic  transistor  flipflop 

equations. 


"f  In  a  common  emitter  coupling  configuration,  i.e.,  we  consider  a  general 
Eccles-Jordan,  and  the  asymmetric  flipf lops. 
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Next  the  problem  of  approximating  the  system  is  discussed;  a  set  of 
approximate  differential  equations  is  established  on  a  basis  of  piecewise 
linear  approximation  to  the  original  system. 

This  set  of  approximate  equations  is  then  analyzed  in  the  phase-plane 
using  the  time  domain  whenever  expedient;  this  analysis  includes  the  study  of 
existence  and  nature  of  singular  points  and  the  description  of  separatrices . 
Then  follows  a  discussion  of  triggering  conditions,  establishing  their  relation- 
ship with  the  system  parameters.   A  specially  detailed  discussion  of  trajectories 
and  their  relation  to  parameters  and  trigger  is  presented.   Some  theorems  about 
the  geometry  of  the  system  phase-plane  portrait  are  proved. 

On  the  basis  of  the  above  analysis  some  techniques  are  described  for 
the  actual  computation  of  trajectories,  transition  times,  waveforms,  or  general 
characteristics  of  a  given  flipflop  driven  by  a  given  triggering  circuit. 
Especially  simple  formulae  are  presented  with  the  purpose  of  better  under- 
standing the  qualitative  influence  of  the  various  parameters;  these  formulae 
are  also  useful  when  a  quick  estimate  of  some  transition  characteristics  is 
needed. 

Several  extensions  to  the  theory  are  also  discussed,  such  as  three 
degenerate  cases,  the  nonsymmetric  Eccles- Jordan,  and  other  types  of  trigger. 

Finally,  some  experiments  are  reported  on  and  their  results  compared 
with  the  theoretical  results. 


2.   THE  FLIPFLOP  DIFFERENTIAL  EQUATION 

2.1   Introduction 

We  shall  try  to  find  an  equation  which  will  not  only  describe  the 
flipflop  behavior  but  which  is  also  relatively  simple.   This  is  a  rather 
difficult  task,  since  the  number  of  equations  describing  such  a  circuit  is  not 
finite.   However,  if  we  add  the  further  requirement  that  the  state  variable  be 
a  natural  circuit  variable,  i.e.,  either  a  node-to-node  voltage  or  a  branch 
current,  then  the  possible  number  of  such  equations  is  reduced  to  about  a  dozen 

possibilities . 

A  comparison  among  the  most  promising  ones  leads  us  to  select  the 
base-to-base  voltage  as  the  circuit  state  variable.   The  ensuing  analysis 
illustrates  the  relative  simplicity  of  such  a  choice. 

2.2   The  Transistor  Pair  Transfer  Equation 

Analysis  of  the  transistor  pair  presented  in  Fig.  1  furnishes  the 
following  system  of  equations,  which  represents  an  entirely  general  description 
of  the  circuit  both  statically  and  dynamically  [2,  15]= 

hi  =   §E1  +  %1   +  JEl^^El  "  1}  -  aClJGl^XpT1VCl  "  X)     (2'1} 
^2  =  §E2  +  §B2  +  JE2(eXpT1VE2  "  l}  "  aC2J  C2(eXpT1VC2  "  l}     {2'2) 

V  -  -°%i  +  aErWexpTlvEi  '  x)  "  Jci(expTlvci  "  1}  (2*3) 

V   =    -§C2   +  aE2JE2^X^VE2    "    ^    "   J  C2  ^VC2    "    l}  (2J° 

*B1   +  °qCl+    {1~  ^El^El^^El 


+  8       +   (1  -  QL-  )JW  (exptlv-.    -  1)  -    (1  -  «cl)JC1(e^vci  -  1} 


(2.5) 
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FIGURE  1:    LOADED  COMMON  EMITTER  COUPLED  TRANSISTOR  RWR. 
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+  (1  -  oOjTO(exPTlv__  -  1)  -  (1  -  a&h&ieWcz  "  ^ 


^2  =  qE2  +  qB2  +  QC2  +  U  "  UE2yjE2V^'VE2   *'    ^   ~C2"C2' 


^1  +  *K>  *  ^ 


VE2  '  VE1  =  V 


VC1  =  EC1  +  ^Cl  "  vi 


VC2  =  EC2  +  ^02  "  V2 


Vl  "  V2  =  V 


(2.6) 
(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.H) 


In  the  system  above: 

a)  The  symbol  x  means  |p-;  and  t1  is  the  time  variable. 

€  =  electron  charge 

b)  r\  =  JL  f   with  <  k  =  Boltzmann  constant 

VT  =  absolute  temperature 

c)  R_  .  R   are  the  load  resistances. 
;   HL1'   L2 

d)  a   ,  q   ,  a„  are  the  charges  stored  in  transistor  T  capacitances: 
;  qEk>  4Bk'  4Ck  K 

q   at  emitter  diffusion  tail  and  depletion  layer  of  the  emitter-base 
Ek 

junction. 

q   at  the  base  diffusion  capacitance. 
Bk 

q   at  the  depletion  layer  of  the  collector  to  base  junction. 

UK. 

e)  i   ,  i   are  the  saturation  currents  of  transistor  T  ,  respectively  of  its 
'  JEk'  dCk  K 

emitter-base  and  collector-base  junctions,  measured  with  the  opposite 
terminal  short  circuited  to  the  base. 

f )  a      ,  a       the  normal  and  inverse  alphas  of  transistor  T  respectively  from 
'        Ek'   Ck 

emitter  to  collector  and  from  collector  to  emitter. 


-7" 

Obs.:      Notice  the  unusual  convention  of   signs,    which  has   been  adopted  in  this 
report  for  convenience   only. 

Comments:  The  system  has  11  equations  and  13  variables;  usually  variables  v 
and  v2  would  be  independent,  and  then  any  other  variable  can  be  expressed  as 
functions   of  them. 

However,    this  system  of  equations   can  be "considerably  reduced  by 
making  the  following  simplifying  assumptions    (see  Fig.    2). 

The  transistors   and  operation  points   are   such  that  the  following 
relations   hold,    accordingly  simplifying  equations    (2.l)  through   (2.6): 

JE1(exPT]vE1  -  1)  »  aclJcl(expTivcl  -  l)  *  -a^  (2.12) 

JE2(eXpllVE2    "    1}   >y  aC2JC2(eXpT1VC2    "   l}   *  "^02  <2'13) 

aE1JE1(exprivE1   -   1)  »  jcl(expr)vci   -   l)  *  -J^  (2,14) 

aE2JE2(exprivE2   -   1)  »  JC2(expr]vC2   -   l)  «  -J^  (2.15) 

(1  "  "ei^ei^^^ei  "  1}  y>  (1  "  aci)jci(expT]vci  '  1}  *  "(1  "  aci)jci 

(2.16) 

(1   -  aE2)JE2(eXpT1VE2    -!)»(!-  a C2h C2(eX^V C2   '   l}   ~  "^   "  ac2)jC2 

(2.17) 

For  each  transistor,  the  following  relations  hold,  accordingly 
simplifying  equations  (2.l)  through  (2.6): 

qEl  <<C   qBl    and    *C1  ^  qBl  ^2-18) 

qE2  <<C  qB2    and    ^C2  K<   qB2  (2'19) 
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+ 

*J)Eci 


+ 


C2 


0 


rift 


FIGURE  2:     UNLOADED  TRANSISTOR  PAIR,WITH  JUNCTION  CAPAClTAKCl 
NEGLECTED. 
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qBl  =  TCliCl    and    %>  =  TC21G2  (2-20) 

where  T Q.    is  the  "collector  time  constant"  of  transistor  T..   The  transistors 
are  such  that  the  following  relations  hold,  accordingly  modifying  equations 
(2.1)  through  (2.6): 

TC1  =  TC2  =  T  (2.21) 

°E1  =  aE2  =  a  (2.22) 

JE1  =  JE2  =  J  (2.23) 

T,  a,  j  independent  of  any  system  variable.  (2.24) 

Assumptions  (2.12)  through  (2.17 )  divide  the  system  (2.l)  through 
(2.11)  into  two  interdependent  systems,  namely  (2.1)  through  (2.8)  and  (2.9) 
through  (2.11).   The  first  has  eight  equations  and  nine  variables,  so  that  all 
variables  are  determined  if  one  of  them  is  given;  v  is  the  "natural"  independent 
variable.   The  second  system  has  three  equations  and  seven  variables,  so  that 
four  variables  must  be  specified;  in  general,  v±   and  vg  would  be  given  (thus 
specifying  v  by  (2.ll))  and  also  i^  and  i^  which  are  solutions  of  the  first 
system  of  equations. 

These  approximations  restrict  our  analysis  to  unsaturating  matched 
transistors.   Constancy  of  parameters  with  respect  to  system  variables  as  well 
as  negligible  junction  capacitances  are  fairly  strong  assumptions,  since  they 
cannot  occur  in  practice  [l5];  however,  we  feel  that  the  increased  complexity 
involved  in  trying  to  take  into  account  such  things  is  too  high  a  price  to  pay 
for  the  small  increase  in  accuracy  to  be  obtained. 
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We  obtain  the  system  of  equations  below  by  substituting  (2.12)  through 
(2.2*0  into  equations  (2.l)  through  (2.8): 


o     l  . 

XE1  "  T1C1  +  a  a"Cl 


o     1  . 
XE2  =  TlC2  +  a  1C2 


icl  =  ajEexpTlvE1 


ic2  =  aJEexprivE2 


1  -  a  . 
hi   =  T1C1  +  ~^T   xci 


1  -  a  . 
1B2  =  TlC2  +   a     '  XC2 


iEl  +  XE2    -""E 


VE2  '  VE1  =  V 


(2.25) 

(2.26) 
(2.27) 
(2.28) 

(2.29) 

(2.30) 
(2.31) 
(2.32) 


Here  we  have  neglected  all  terms  jQ^   ac^c^    i1  "   ack^Ck'  in  the 
same  spirit  we  have  used  the  approximation 

JEk'eXpT1VEk*  JEkfeXpT1VEk  "  X)  (2'33) 

which  is  not  strictly  valid.   Actually,  it  is  only  valid  as  long  as  vE1  is 
positive,  but  not  if  ^  is  negative.   However,  to  use  this  approximation  (or 
should  we  call  it  a  substitution!)  for  all  the  range  of  vEk  is  equivalent  to 
connecting  a  current  source  of  strength  j£k  in  parallel  with  the  emitter-base 
junction  of  T,  ,  in  a  direction  such  as  to  yield  input  current  zero  when 


v 

Ek 
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This  will  make  very  little  difference  as  far  as  results  are  concerned, 
but  will  considerably  simplify  the  algebra  involved.   To  be  sure,  it  will  tend 
to  compensate  for  the  fact  that  we  have  already  neglected  the  collector-base 
junction  saturation  current,  yielding  even  a  better  composite  characteristic 
for  the  transistor  pair. 

Equations  (2,9)  through  (2.1l)  are  now  irrelevant  since  we  seek  to 
express  the  currents  as  functions  of  v,  which  can  be  obtained  just  by  solving 
system  of  equations  (2.25)  through  (2.32). 

For  convenience,  we  apply  the  following  transformations  of  variables 
to  this  system: 


t=T  (2-34) 


XCk 
w,  = 


k   a 


V 


^Bk 
zk  -  aIE  > 


k  =  1,  2  (2.35) 


k  =  1,  2  (2,36) 


X  =  I  TJV  (2.37) 


RQ  is  the  output  resistor  as  shown  in  Fig.  3, 

From  now  on,  the  symbol  x  will  be  used  for  ^  ,  unless  otherwise 


specified.   We  obtain: 


rO     o 


a(w1  +  w2)  +  (wx  +  w2)  -  1  (2.38) 

W2  =  w1exp2x  (2.39) 

W2  =  (Wl  +  2xw  )exp2x  (2.40) 
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FIGURE  3:    ASYMMETRIC  FLIPFLOP 
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o    1  -  a 


o    l  -  a 
Z2  =  V2  +  "IT"  W2  (2.^2) 

from  which  we  get,  using  (2.38),  (2.39)  and  (2.40): 

^  +  j-  +  x(l  +  tanh  x)|  w1  =  I  (1  -  tanh  x)  (2.43) 

This  equation  can  be  written  in  the  form: 

Mdt  +  Ndw  =0  (2.kk) 

If  this  is  done,  then  we  have: 


\5^  "  5t  J  =  ja  +  x(l  +  tanh  x)j  "  «(*)  (2.45) 


where  g(t)  is  only  a  function  of  t,    i.e.,    it  is  independent  of  w  ;    here  we 
assume  that  x  ■■■■   x(t)  is  a  given  function  of  t,  independent  of  w  . 


Then  u(t)  =  exp 
(2.^3)  [8].   We  have  then: 


g(t)dt 


is  an  integrating  factor  for  equation 


u(t)  -  exp  -g  +  1(1   +  tanh  x)xdtf  (2.46) 


But 


j(l  +  tanh  x)xdt  =  /(l  +  tanh  x  )dx  =  x  +  SmJ  (cosh  x)  +  const.     (2A7) 
Hence,  ignoring  the  constant, 


t 
XQ!  1  -  tanh  x 


,                  exp  — 
u(t)  =  exp(-  +  x)  •  cosh  x  - ~L — "  (2.48) 


-Ik- 


The  differential  equation  can  now  "be  easily  solved, 
Let  F(-w  ,t)  =  0  "be  the  solution;  then: 


dF 
3w\ 


=  u(t)      (since  M  =  l) 


and 


£-*<*> 


SI 


(2.U9) 


We  get  from  (2.46) 


F(w^t)  =  wx  •  u(t)  +  f(t)  =  0 


(2.50) 


and  from  (2.^9)  and  (2.50) 


N|i(t)  =  wg  •  ji(t)  +  ?(t) 


(2.51) 


Therefore 


f(t)  =  -  |  (1  -  tanh  x)u(t) 


(2.52) 


So 


f(t)  =  -  |  /(l  -  tanh  x)  '  u(t)dt 


(2.53) 


From  (2.U8)  into  (2.53) 


f(t)  =  -  (§  exp  -  +  const.) 


2  c^  a 


(2.510 


From  (2.50),  explicitating  w  ,  we  finally  get 


w  =  I  (1  -  tanh  x)[l  +  Cexp(-  -)] 


(2.55) 
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where  C  is  an  arbitrary  constant  determined  by  some  initial  charge  stored  in 
the  bases  prior  to  the  moment  t  =  tQ  when  the  circuit  takes  on  the  configuration 
described  by  (2. 1+3).   Therefore,  for  all  practical  purpose  Cexp(-  -)  is  zero, 
since  we  can  assume  this  configuration  to  be  in  existence  for  an  arbitrarily 
long  time;  we  finally  have,  using  (2.56),  (2.39),  (2.4l)  and  (2A2): 


The  "transfer  functions": 


w  =  —  (l  -  tanh  x) 


2  =  -  (l  +  tanh  x) 


and  the  "input  functions": 


z.  = 


1  1  -  a 


1   2   a 


(l  -  tanh  x)  -  x  •  sech  x 


2   2    a 


For  completion,  we  have,  with  k  =  1,  2 


o    f      Nk  1 

Zk  =  (~^  2 


1  Jl  -  a  /,   .  ,   x   o     ,2 


(l  +  tanh  x)  +  x  •  sech^ 


1  -  a   o   00   ^/0\2  , 

x  +  x  -  2(x)  tanh  x 


a 


,2 

sech  x 


(2.56) 


(2=57) 


(2.58) 


(2.59) 


(2.60) 


which  are  the  characteristic  functions  for  a  transistor  pair  connected  as  in 
Fig.  2. 

As  a  corollary  of  this  solution,  we  notice  that,  at  all  times,  the 
following  relations  hold: 
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Wl  +  W2  =  1 


=>   ^1  +  ^2  =  aIE 


^^2-^r1    =>    H^h.  =  ^--Kj 


(2.61) 


which  check  with  (2.7):   iEl  +  iE2  ~  ^^ 

2„3  The  Asymmetric  Flipflop 

The  asymmetric  flipflop  shown  in  Fig.  3  has  a  feedback  network  whose 
configuration  depends  on  the  trigger  circuit,  as  shown  in  Figs,  ha   and  kh. 

However,  we  can  obtain  the  equilibrium  equation  for  this  network 
assuming  an  arbitrary  trigger  current  i,  which  can  be  replaced  by  its  value 
when  a  trigger  is  specified  (see  section  2.5). 

Analysis  of  the  circuit  of  Fig.  ^a  yields: 


Vo°v°i  +  (Ti  +  Toi  +  To)\  +  vi  =  R^  +  V„  +  Ec +  Vio(?  +  °bi'  ^  "s-  '  -a 


/O  O  /  .         .         \ 

. .    (i  +   i„.  )  +  R   (i  +  ip1  ; 
o    C2  B  s  C  o   iov 


(2.62) 


where 


a)       T     -  R.C. 


b)       To   =  RoCo 


c)     T.      =   R.C 
'        io  i   o 


d)     T    .    =   R  C. 
'       oi         oi 


e)       R     =  R.    +  R 

'  s  1  o 


(2.63a) 
(2.63b) 
(2.63c) 
(2.63d) 
(2.63e) 
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a)  FEEDBACK  NETWORK  AND  EQUIVALENT  CURRENT  SOURCE 
FEEDING  TRIGGER  CURRENT  it. 
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b)    THE  CURRENT  SOURCE  WITH  AN  INTERNAL  CONDUCTANCE 
6  IS  REPLACED  8Y  AN  IDEAL  VOLTAGE  SOURCE  WHOSE 
OUTPUT  IS  r. 


FIGURE  4-  TRIGGER    SOURCES 
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Notice  that  in  this  asymmetric  flipflop  (Fig.  3)  ^  a  v.   Let  us 
apply  to  (2.62)  the  same  transformations  of  variables  used  in  section  2.2,  i.e., 
equations  (2.3*0  through  (2.37),  and  also  (2.6U)  below  (a  similar  transformation 
for  the  trigger  current): 

(2.61+) 


a 


h 


dx 
The  result,  with  x  again  standing  for  —  ,  will  be 


T  T       T   +  T  .  +  T 

i  o  00    i    01    o  o 
x  + x 


+  x  =  p  {2B  +  2w2  +  2  ^   (8  +  §x)  +  2  ^  (0  +  Z]_)} 


2  T 

(2.65) 


where 


a)  B  = 


TBRs  +  EC  (2.66a) 


E  o 


h)     0   =   1  ~  a   -i    (not  explicit  above)  (2.66b) 


«0  p-^wft  (2-66c) 

If  we   expand    (2.65)  by  replacing   z      according  to    (2.58)  we  get: 


T    T  T       +    T     .     +    T 

i    OOO           1            01            O    O 
—  x+ x  +  x  ?) 

T 

T  R  "] 

-  p  {(1  +  2B  +  p)  +   (1  -   P)tanh  x  +  2  -^  (0  +  8^  +  2  ^  (S   -   zQ)j 


where 
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z  =  -  x  sech  x  (2.68) 


The  condition  given  by 

1  +  2B  +  p  =  0  (2.69) 


is  called  state-symmetry. 


11 


In  practice  the  effect  of  the  base  current  on  biasing,  represented 
by  the  term  pp, '  is  usually  small  (see  (2.66b)),  so  that  the  state-symmetry 
condition  can  be  approximated  by  setting  1  +  2B  =  0.  Of  course  the  larger  Ot 
is  (closer  to  l),  the  better  this  approximation  will  be. 

Equation  (2.65),  whether  in  this  form  or  its  more  explicit  form 
(2.67),  is  the  general  asymmetric  transistor  flipflop  equation. 


2.4.   The  Eccles -Jordan  Flipflop 

Analysis  of  the  two  feed-back  circuits  of  the  flipflop  of  Fig.  5 
leads  to  the  following  pair  of  differential  equations : 


T..T  T-   +  T     +  T 

ik  ok  00     1k    01k    ok  o 

O X1  +  Z XT     +   XT 

,2    k  t        k    k 


(2.70) 


=  Pk^2Bk  +  2w,  +  2!^(§k  +  §k)+2^(0   +z  ) 

ok 


rith  k  =  1,  2;  I  =  1,   2;  I  £   k;  t  as  in  (lU.l)  and  (ik.k)   and: 


a)  Tik  -  Rikcik  <2-71a) 


Notice  that  this  term  implies  that  the  base-to-base  voltage  i,v  (in  this  case 
identical  to  v  )  will  have  a  constant  component  equal  to   1(1  -  a)]L,R  . 

1  P  J±i  s 
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b)       T  ,    =  R  ,  C  .  (2.71b) 

ok         ok  ok 


c)     t  =  R     C  ,  (2.71c) 

iok  lk  ok 


d)     T    ..    =  R..  C,  (2.71d) 

'       oik         ok   lk 


e>     Rsk  "  Rlk  +  Rok  <2-71e> 


I_,n  R  ,    +  E_ 
f)  B=-^ ^  (2.71f) 

2  aIERok 


)  \  =  |  T)vk  (2.71g) 


h)        pk  =  J  iaIER0,k  (2.7lh) 


i)         0,    =  -4-  (2.711) 

E 


j)     w     and  z      are  given  by  equations    (2.56)  through    (2.59)  with 

x  replaced  by   (x      -  x    ).  (2.71j) 

Notice  that  i   is  the  trigger  current  actually  fed  into  the  flipflop; 
i   is  the  trigger  current  put  out  by  the  pure  current  sources  (see  Figs.  K). 
Executing  (2,71j)j  w  and  z  turn  out  to  be: 

it  K. 


Vk    = 


I  jl  +    ("l)k  tanh(X;L   -  x2)j"  (2.72) 


1   \l_Zm  a 
\      2  1     a 


j      ~a  a   [1  +    (-1)     tanh(x1   -  x2)]    +    (-l)    (x^   -  x2)sech    (x1   -  xg  )| 

(2.73) 
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4 
Taking  (2.72)  and  (2,73)  into  (2.70)  we  get: 


+  xk  =  ?k  {&  +  2Bk  +  Pk} 


TT         T    +  T  .,  +  T  , 

ik  okoo     ik    oik    ok  o 

X    + X, 

k   "   -  (2.7U) 

-  (-l)k(l  -  Pk)tanh(x1  -  xe)  +  2  ^  (g  +  §  )  +  2  !i*  (^  +  (-1)%) 
^  ok 

where,  with  k  =  1,  2: 

a)  p  =  I   -  a-  ^  (2.75a) 

a;  pk    a   R  . 
ok 

b)  zQ  -  i  (%  -  S2)seeh2(Xl  -  x2)  (2.75b)   ; 

Equations  (2.7*0  describe  a  general  Eccles -Jordan  transistor  flipflop. 
The  asymmetric  flipflop  equation  (2.67)  is  a  special  case  of  (2.7*0,  and  so  is 
the  symmetric  flipflop,  represented  by  equation  (2.77),  as  follows: 

Let  us  consider  the  case  of  a  symmetric  (but  possibly  unsymmetrically 
biased)  flipflop.   Subtraction  of  (2,70,  k  =  2)  from  (2.70,  k  =  l)  yields: 


T    T                  T.    +   T     .     +    T 
iooo          1           Ol           OO 
X   + X 


+  x  =  p  J2(B1   -  B2)  +  2(w1   -  w2) 


T2  T 

T  (2.76) 

+  2  ^  [(§-   -  e)  +  (§x  -  §2)]  +  2  ^  [(ex  -  e2)  +  (z,  -  Zg)]} 


This  furnishes  (2.76)  below,  which  can  also  be  obtained  directly  by 
subtraction  of  (2.7*+,  k  =  2)  from  (2.7*+,  k  -  l): 


Notice  that  (-l)   =  -(-l)  • 
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T  .T  T  .     +    T     .     +    T 

1    O  OO  1  Ol  O    O 

— —  x   + X   +  X 


T 


T 


T.  o  R 

2(B1   -    B2)   +  2(1   -   p)tanh  x   +  2  -^    [0    +   z]    +  2  ^   [©    -   2zQ] 


(2.77) 


where  the  indices  1,  2  have  been  dropped  from  all  parameters  P  whenever 
P  =  P  ,•  in  the  case  of  the  variables  we  have  as  before: 


a)  x  =  x.,  -  x. 


(2.78a) 


b)  e    = 


(2.78b) 


c  )   z 


Zl  -  Z2  = 


tanh  x  +  2z~ 

a  0 


j    from  (2.73) 


and  (2.78a) 


(2.78c) 


d j   z   =  —  x  sech  x 


(2.78d) 


Inspection  of  (2,7^-)  and  (2.75)  shows  that  the  symmetric  flipflop  is 
formally  equivalent  to  the  asymmetric'  one.   Furthermore,  if  it  is  biased  in 
such  a  way  that  B  =  B   (not  necessarily  symmetric  biasing),  then  it  is 
formally  equivalent  to  a  state-symmetric  asymmetric  flipflop. 

We  should  stress  the  importance  of  this  conclusion,  since  it  actually 
doubles  the  effectiveness  of  the  theory!   But  better  yet,  it  allows  us  to 
compare  the  performance  of  any  symmetric  flipflop  with  its  asymmetric  counter- 
part, for  we  will  actually  be  comparing  formally  equivalent  things. 


In  this  report,  the  word  "asymmetric"  is  reserved  for  the  flipflop  of  Fig.  3» 
The  Eccles-Jordan  with  different  parameters  on  either  side  will  be  called 
"nonsymmetric . " 
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We  can  anticipate  here  that  (as  will  be  seen  in  the  following)  that 
all  other  things  being  equal,  the  asymmetric  flipflop  besides  being  simpler, 
performs  better  than  the  symmetric  flipflop  with  the  same  parameter  values. 

2.5  Triggering 

The  equations  we  have  found  so  far  take  the  trigger  into  account  in 
a  rather  general  form.   Among  the  infinite  possible  trigger  waveforms,  however, 
we  shall  discuss  only  those  which  are  the  most  important  in  practice,  and 
study  them  in  some  detail. 

First  we  consider  the  type  of  trigger  circuit.   It  can  be  a: 

a)  Voltage  source  with  negligible  internal  resistance. 

b)  Current  source  with  negligible  internal  conductance. 

c)  Voltage  source  with  considerable  internal  resistance. 

d)  Voltage  source  with  considerable  internal  conductance. 

So  for  the  trigger  current  waveform,  the  following  are  three 
important  types  (see  Fig.  6): 

a)  Rectangular. 

b)  Impulse. 

c)  Pair  of  rising  and  decaying  exponentials. 

d)  Sinusoidal  wave  plus  constant,  simulating  usual  trigger. 

We  will  study  type  a)  in  detail,  and  discuss  types  b),  c)  and  d) 

(see  section  5-5 )• 

The  asymmetric  and  symmetric  flipflops  will  be  treated  together 
(since  they  are  formally  equivalent),  and  an  indication  will  be  given  as  to 
how  a  basically  similar  method  applies  to  the  more  general  nonsymmetric 
Eccles- Jordan  flipflop.   This  is  justified  because  of  the  much  greater  practical 
importance  of  the  asymmetric  and  symmetric  case,  compared  to  the  more  academic 
importance  of  the  nonsymmetric  case. 


a, 


•25- 
0=0 


to+T 


a)   RECTANGULAR      TRIGGER. 
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b)    IMPULSE      TRIGGER. 
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c)    PAIR  OF  RISING  AND  DECAYING  EXPONENTIALS. 
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d)    SINUSOIDAL    WAVE   PLUS  CONSTANT, 
SIMULATING    USUAL    TRIGGER. 
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In  the  treatments  of  both  voltage  and  current  triggering  sources,  we 
will  assume  that  the  trigger  has  a  definite  duration,  outside  of  which  current 
sources  have  zero  output  current,  and  voltage  sources  are  isolated  from  the 
input  of  the  flipflop  (say,  by  means  of  an  open  diode  gate). 

One  way  of  simulating  this  condition  is  to  require  that,  outside  the 
duration  of  the  trigger,  the  voltage  or  current  of  the  triggering  sources  are 
such  that  the  normalized  trigger  currents  9^   be  zero,  for  any  k. 

Let  us  consider  the  various  types  of  trigger:   type  a)  is  a  trivial 
case,  since  the  variable  x  is  specified;  case  b)  is  both  theoretically  and 
formally  a  special  case  of  d)  and  case  c)  is  shown  to  be  formally  equivalent 
to  d),  which  is  not  surprising  since  one  case  can  be  converted  into  the  other 
by  Thevenin  or  Norton  transformations. 

We  will  establish  expressions  for  9    in  each  of  the  special  cases 
mentioned  above  except  case  a). 

Case  b) :   See  case  c). 

Case  c):   Let  the  trigger  circuit  be  as  shown  in  Fig.  7a.   We  get, 

with  k  =  1,  2  : 


G,R  . 
k  ok 

tk  "  2p.   Xk 


_  _     G,  R  , 

o  o      k  ok  o 

9      -   9        -  —r x, 

k  tk    2p,    k 


(2.79) 


(2.80) 


where 


)   =  -^    is  the  normalized  trigger  current  put  out  by    (2.8la) 
E    the  ideal  current  source. 


Q     =  —  is  the  normalized  trigger  current  actually  fed   (2.8lb) 

k   aIE    into  the  flipflop  terminals. 
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FIGURE  7a:  TRIGGER  CURRENT  ik  PRODUCED  BY  CURRENT  SOURCE 
CLEARLY,  ifc—  *tk~"Gkvk- 


FIGURE  7b:  TRIGGER  CURRENT  ik  PRODUCED  BY  VOLTAGE  SOURCE. 
CLEARLY, 


lk"TT^tk"vK,"Hk'"Gkvk 


•Pifk'GkVtk    AND    Gk  =  ^ 


FIGURE  7:     TRIGGER  SOURCE 
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G     is  the  trigger  circuit  internal  conductance,       (2.8lc) 
k 

Nov  case  b)  is  obtained  by  letting  Gk  -  0  in  equations  (2.79)  and 

(2.80). 

Case  d):   Let  the  trigger  circuit  be  as  shown  in  Fig.  fb .  We  get, 

with  k  =  1,  2: 

R  R  , 

ok  „       ok  _  „  (2.82) 


x . ,  -  ~ — =-  x, 


k   2pA  tk   2pA  k 


gv  -  _i_  8   -  _i- £  (2.83) 


where 


Xtk   2  nVtk 


-  nv      is  the  normalized  source  voltage.       (2.84) 


He 


nee,  independently  of  the  nature  of  the  trigger  circuit  the  form 


of  6      is  the  same.   In  other  words,  the  voltage  source  x^  with  internal 
resistance  R  can  be  converted  into  a  current  source  of  intensity  6^  and 
internal  conductance  G,  given  by: 


P. 


ok  (2.85) 


x, 


tk   2pkRk  tk 


r*  _  _L  (2.86) 

k 


From  now  on  sk  will  be  used  meaning  either  9^   or  e*k,  whatever  the  case  may 
be. 
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Furthermore,  we  will  require  that  whatever  the  waveform  of  s  ,  it 

must  be  nonzero  only  for  a  time  interval  (t   ,  t   ),  being  zero  at  any  other 

aK   DlC 

time. 

From  now  on  we  will  use  only  s,  and  G,  to  describe  the  trigger. 

For  the  asymmetric  flipflop,  we  will  have  only  s,  and  G,  ,'  for  the 
symmetric  Eccles-Jordan  s   -  sp  =  s,  and  we  require  G,  =  Gp  =  G.   For  the  non- 
symmetric  Eccles -Jordan  we  will  keep  s   and  G,,  k  =  1,  2,  in  their  respective 
independent  equations . 

We  can  also  account  for  the  possibility  that  G  has  one  value  when 
the  trigger  is  ON  and  another  when  it  is  OFF  by  using  an  index  \x   to  indicate 
which  case  is  being  considered:  Then  u-  =  0  or  1  will  indicate  respectively 
trigger  OFF  or  ON,  and  G,   will  assume  the  appropriate  one  of  its  two  possible 

KU 

values,  G,  ~  or  G,  ,  . 
'   kO     kl 

With  this  notation,  the  general  flipflop  equations  are: 
For  the  asymmetric  case : 


T.T       T.  +  T  .  +  T   +  T.  R  G, 

1  O  OO     1     01     o     io  o  tu  O    /_,  _  _  % 

— r-  x  + =  x  +  (1  +  R  G  )x 

T2                    T                 .   v  s  \1J 


J  r-i        rs-n  \        /-,  w     i  io  oo      f{±  -  a)     io  s  \ 

pj(l  +  2B+p)+(l-p  )tanh  x  -     —  x  +  ( A— ^ — ~  ~  +  R~  ) 


„      io   ,0\2,      , 

2  — — -   (x)  tanh  x 


sech2x  +  20Hg  +f7)  }  (2.87) 


For  the   symmetric  Eccles-Jordan  case    (not  necessarily  symmetric 
trigger ) : 
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TT  T      +T.+T      +T.RG, 

i  QQ9   ,      i  01 2 1Q  °   ^  i  +    (1  +  R  G    )x 

2     x  +  T  s   u 


=  p 


2(B     -   B   )   +2(1  -   p)tanh  x   -  2  I  — 


Tio  oo       /"(l  -a)     io  s\o 

-IT-  —  +  rJ 


io    /0\2,      , 
—  (x)  tanh  x 


/T.  R        N 

2  r-/     io   o  s 

sech  x  +  2(  —  s   +  p-  s 


[2.88) 


where 


x   =  x      -   x2 


s     —    S,     -     Sp 


G     -  G,      -   G0 
u  l|i  2p 


(2.89a) 
(2.89b) 
(2.89c) 


All  other  coefficients   being  the   same  for  both  (2.89d) 

values   of  circuit   index  k,    the   index  has  been  dropped. 

For  the  general  nonsymmetric  Eccles -Jordan  flipflop: 


T..T 


ik'ok    oo        Tik  +  Toik  +  Tok   +  Tiok   ok  tk)i  o  ,±  +  R     Q      v 

— 2 —     xk  +  T  k        v  sk  k|i 


pk  |(1  +  2Bk  +  Pk)   -    (-l)k(l  -   Pk)   tanh  x 


("DJ 


"Tiokoo      /(l  -  a)  liok       ^sk\  o   _  lipk   (o}2tanh         . 


,2 
sech  x 


/\ok  o  sk 


k  =   1,    2 


(2.90) 
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2.6  The  Approximation  Problem 

We  have  thus  established  the  nature  of  flipflop  behavior  through 
equations  (2.65),  (2.70),  (2.74)  or  their  respectively  equivalent  forms. 

It  is  clear  that  the  chances  of  success  of  an  attempt  to  solve  such 
equations  exactly  are  very  slim  indeed  [18]. 

Therefore,  in  order  to  get  useful  results,  we  must  bow  and  try  to 
find  approximate  solutions  to  these  equations. 

There  are  two  possible  approaches  to  this  problem.   One  is  to  find  an 
approximate  model  to  the  circuit.   This  model  must  be  described  by  solvable 
equations.   We  then  consider  such  solutions  as  approximations  to  the  exact 
solution  of  the  original  problem.   The  other  approch  is  to  approximate  the 
original  equations  directly,  rather  than  the  model.   The  latter  approach  would 
probably  allow  more  accurate  results  to  be  obtained,  since  an  error  introduced 
in  approximating,  say,  iCk(v)  would  not  necessarily  propagate  through  ipk(v)  to 
iBk(v),  and  then  to  ^(v).   That  is  to  say  that  each  function  would  be 
calculated  exactly  for  our  original  model  and  then  each  one  of  them  would  be 
approximated  as  well  as  possible. 

However,  as  we  will  show,  this  question  does  not  necessarily  affect 
the  nature  or  even  the  complexity  of  the  equation  or  system  of  equations  we 
must  analyze.   For  in  fact,  since  an  exact  solution  is  not  to  be  found  but 
instead  we  must  be  content  with  an  approximation,  we  might  as  well  linearize 
the  problem.   This  would  involve  dividing  the  x-axis  into  regions  where  the 
function  of  x  appearing  on  the  right-hand  side  of  the  equation  can  be  approxi- 
mated by  a  linear  function  of  x.   Let  N  be  the  number  of  such  regions.   Then 
we  would  reduce  our  problem  to  that  of  solving  N  linear  second  degree 
differential  equations  valid  in  their  respective  regions  and  match  the  solutions 
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at  the  "border  between  every  pair  of  adjacent  regions.   On  both  the  phase  plane 
and  time  domain,  these  regions  are  N  strips  separated  by  (N  -  l)  lines  of 

x  =  constant. 

One  of  the  great  advantages  of  this  type  of  approximation  is  the 
relative  simplicity  of  its  applications,  whether  in  the  problem  of  analyzing 
the  properties  of  the  circuits  involves,  or  in  the  actual  computation  of 
transition  times  and  waveforms. 

If  this  direction  is  chosen,  then  we  have  decided  to  pay  a  price  in 
exactness  for  the  advantages  of  simplicity  and  usefulness.   Then  it  is  a 
question  of  inspection  to  see  that  the  most  promising  way  to  arrive  at  con- 
sistent results  is  to  approximate  the  transistor  pair  model,  i.e.,  equations 
(2.56)  through  (2.59)  must  be  the  basis  for  the  definition  of  a  "reasonable" 
piecewise  linear  model  [26]. 

Of  course  the  problem  centers  on  how  to  linearize  tanh  x.   There  is  a 
good  degree  of  arbitrariness  here,  but  we  will  select  the  following  three 
region  approximation  (see  Fig.  8). 


tanh  x  w  <p(x) 


(2.91) 


-l; 


cp(x)  -  J.    7x 


+1, 


x 


7 


-  <   x  <  +  - 
7  -       ~       7 


+  i<x 
7   - 


(2.92) 


where  7  is  a  factor  which  depends  on  the  criterion  for  the  approximation.   In 
general  7   will  be  somewhere  between  0.5  and  1.0, especially  if  a  minimum 
integral  error  is  sought.   In  this  case,  we  would  have; 
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tanh  x  a    <£(x) 
A 


tonhx 


TANGENT  LINE  TO   tanh  x    AT  x=0 
tonhx 


-OX 


a)  THE  PIECEWISC 
APPROXIMATION  OF 
tanhxBY^U) 


-OX 


to  THE  RRST  DERIVATIVE 
OF$(x) 


oX 


c)  THE  SECOND 
DERIVATIVE  OF*(x) 


FIGURE  8:     THE  APPROXIMATION  OF  tanhx  BY$(X) 


-3^- 


-j  /   (tanh  x  -  7x)dx  +  /   (tanh  x  -  l)dxj  =  0 


lim 

0 

7 


$mJ  cosh  |  -  ~  -!+yj  ^^-^-O 


i.e. 


7  =  —tF-^       =>   7  -  0.721 

This  criterion,  of  course,  is  quite  arbitrary,  and  there  is  nothing 
to  prove  that  it  is  the  best;  however,  it  also  reduces  the  maximum  absolute 
difference  between  the  two  functions  to  about  0.120  at  the  worst  point.   This 
is  not  a  minimum;  if  7  is  selected  to  minimize  this  maximum  absolute  error 
rather  than  the  integral,  then  it  would  be  7  ~  0.71^  and  the  maximum  absolute 

error  would  be  about  0.115. 

So,  for  simplicity,  we  could  for  example,  take  7  =  0.7,  or  alter- 
nately (and  perhaps  better)  take  7  =  7. 

Whatever  the  criterion  may  be,  the  intention  of  using  such  a  factor 
7  is  to  reduce  somewhat  the  error  by  which  the  solution  to  the  differential 
equation  is  affected  due  to  the  piecewise  linearization  process.   This  would  be 
the  ideal  criterion  if  it  were  not  impractical. 

A  five  region  approximation  could  also  be  used  with  perhaps  better 
accuracy,  but  increased  labor  involved  both  in  computation  and  analysis.   In 
general  the  accuracy  can  be  improved  by  increasing  the  number  of  regions,  which 
results  also  in  increased  labor. 

Since  the  technique  does  not  change  in  essence,  we  shall  use  the 
three  region  approximation  in  this  dissertation. 
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Equations  (2.56)  through  (2.59)  become; 


w. 


=  |{i+  (-i)V*)} 


(2.93) 


z,  = 


k   2   a 


o    /  xk  1 

zk  =  ("1}  2 


j^-2  [1  +  (-l)S(x)]  +  (-l)k  Sep' (x)}        (2.9^) 


a 


a  o  00 

X  +  X 


cp'(x)  +  (S)V(x) 


(2.95) 


with  k  =  1,  2.   The  prime  indicates  differentiation  with  respect  to  x.   Also, 


1  o  ,  /  \ 

zQ  =  -  xcp'(x) 


(2.96) 


And  obvious ly, 


0,      x  <  - 


cp'(x)  =  <   yt 


-i<x<+i 
7  7 


°>  +   T  <  x 

'      7 


(2.97) 


Differentiation  of  cp'(x)  will  clearly  consist  of  two  impulse  functions, 
since  cpf(x)  is  constant  everywhere  except  for  two  discontinuities.   We  get: 


<P"(*)  =  7{z(x  +  j)   -   6(x  -£)]• 


(2.98) 


where  7£>(x)  is  an  impulse  function  of  strength  7  occurring  at  x  -  0. 

The  flipflop  equations  (2.87),  (2.88)  and  (2.90)  will  be  approximated 
by  the  following  equations,  respectively: 

For  the  asymmetric  case,  (2.87)  is  approximated  by: 
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TT  T+T.+T+T.     RG 

i   °9?  ,    _i 2i 2 1Q  °  ^  $  +    (1  +  R  G    )x 

2  1  s   M- 


=  P 


(1  +  2B  +  p)  +  (1  -  p)cp(x)   -     —  x  + 


T  /,  T.  R 

io  oo      f  1  -  a.     io         s  \  o 


T      +R-,X 

o 


cp'(x) 


^(8)VW  +  2^^ 


(2.99) 


For  the   symmetric  Eccles-Jordan  case    (also  not  necessarily  symmetric 
trigger),    (2.88)   is   approximated  by: 


T.T 


T-+T.+T       +   T.     P    G 


L2  9°+  -i 2i 2 !2_°_fc!:  £  +   (l  +  R  G   ) 


s    |_l 


p  {2(Bl   -   B2)   +  2(1   -   pMx)   -  2  [>°x°+  (i^2  !±2  +  Qx}p'(x) 


2  ^2  (g)V(x)  ♦  2(^2  8  +  ^  .) 


(2.100) 


with 


X   —   x_    -   x_ 


s   =   Sl   "   S2 


Gu   "   Glu   =   °2u 


(2. 101a) 
(2.101b) 
(2.101c) 


All  other   coefficients   being  the   same   for  both  values  (2.101d) 

of  circuit   index  k,    the   index  has   been  dropped. 
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For  the  general  nonsymmetric  Eccles-Jordan  flipflop,  (2.90)  is 
approximated  by: 


ik  ok  op     ik  +  oik  +  Tok  +  TiokRokGtku.  o    ,  , 


=  Pk{(l  +  2Bk  +  pk)   -    (-l)k(l  -   pk)<p(x) 


+  (-if 


iokoo        1  -  a  Tiok       Rsk     o 

ok   J 


cp'(x) 


(2.102) 


+  (-I)*  ^  (x)V(x)  +  2  Tl°k  9 


T   Sk  +  R^Sk 


Therefore,  equations  (2,65)  and  (2.77)  will  be  approximated  by  an 
equation  of  the  type; 


V°+  V*  +  V  =  dv  +  f^H8(x  +  i) 


-  8( 


x  - 


1m     o 

yjj  +  mvs  +  ns 


(2.103) 


where 


V  -  1,  2,    ...,    N,  is  the  region  index,  and 
V-  -   0,  1         is  the  trigger  index. 


The  system  of  equations  (2.70)  can  be  similarly  approximated: 
\v\   +  \vA  +   °kVUXk  =  akv°X°+  bkVX  +  CkVX  +  dkV  +  fk(*)tB(x  +  h    -    6(x-  i)] 


+  \vsk  +  nksk 


(2.104) 
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where 


x  -  x  -  Xp,    k  -:  1,  2, 


V  and  u  as  for  (2.103)  above 


Obs.:   Notice  that  C.  and  C  may  vary  with  the  state  of  the  transistors.  We 
will  recognize  the  possibility  of  one  value  in  each  region,  and  so,  strictly 
speaking,  T  ,  T  ,  TQi  and  tq  depend  on  V,  the  region  index.  When  necessary 
we  shall  indicate  this  dependence  explicitly,  but  not  otherwise. 

Tables  1.1  and  1.2  contain  the  values  of  all  the  parameters  of 
equations  (2.103)  and  (2.10*0  for  the  case  where  the  number  of  regions  is  three. 

Therefore,  a  detailed  study  of  such  equations  is  useful;  it  informs 
us  about  flipflop  characteristics  and  also  serves  as  a  basis  for  analysis  and 
design  procedures,  optimization  of  trigger,  and  study  of  interaction  with 
adjacent  circuits. 

2.7  Summary 

In  this  chapter  we  have  analyzed  a  general  transistor  flipflop  circuit 

and  discussed  its  equilibrium  equations. 

a)  The  characteristic  equations  for  the  transistor  pair  with  both  emitters 
connected  to  a  common  constant  current  source  is  obtained. 

b)  Analysis  of  the  asymmetric  flipflop  leads  to  a  second  order  nonlinear 
differential  equation. 

c)  Analysis  of  the  general  Eccles -Jordan  flipflop  leads  to  a  system  of  two 
second  order  nonlinear  differential  equations.   It  is  shown  that  both  the 
asymmetric  flipflop  and  the  symmetric  Eccles -Jordan  flipflop  can  be 
represented  by  a  single  differential  equation.   In  the  first  case,  because 
one  of  the  equations  does  not  arise  since  its  corresponding  would-be 
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variable is  forced  to  be  zero  at  all  times.   In  the  second  case,  because  the 
difference  between  the  two  original  equations  is  taken  leading  to  another 
equation  of  the  same  type.   Therefore  the  formal  equivalence  between 
symmetric  and  asymmetric  flipflops  is  established. 

d)  Triggering  is  discussed  and  the  representation  of  triggers  in  the  flipflop 
equations  is  presented. 

e)  Finally,  the  approximation  problem  is  discussed  and  a  course  of  action  is 
decided  upon,  which,  although  somewhat  arbitrary,  seems  to  minimize  the 
unavoidable  error. 


3.   STUDY  OF  THE  FLIPFLOP  EQUATION  FOR  THE  CASE  OF  A  RECTANGULAR  TRIGGER 

3.1  Introduction 

We  have,  so  far,  broken  down  the  problem  of  flipflop  analysis  and 
design  into  the  analysis  of  two  possible  cases  a)  and  b)  below: 

a)  Problem  of  the  fourth  degree:  the  general  nonsymmetric  Eccles- 
Jordan  flipflop,  described  by  a  system  of  two  second  order  nonlinear  differential 

equations . 

b)  Problem  of  the  second  degree:   the  asymmetric  flipflop  and 
symmetric  Eccles-Jordan  flipflop,  which  can  both  be  described  by  the  same 
single  second  order  nonlinear  differential  equation.   This  will  be  called  "the 
basic  flipflop  equation,"  because  its  analysis,  besides  its  greater  practical 
importance,  turns  out  to  be  fundamental  for  the  analysis  of  the  general  case, 
since  the  two  equations  describing  the  latter  are  closely  related  (formally)  to 
the  single  equation  describing  the  former.   It  is  also  basic  (in  a  sense)  in 
obtaining  the  equation  for  the  following  special  case. 

c)  Problem  of  the  first  degree  (internally  restricted  flipflop): 
We  are  considering  only  the  second  case  where  the  transistors  are  identical. 
In  this  case,  if  the  circuit  capacitances  are  negligible,  no  matter  what  the 
symmetry  of  the  flipflop  may  be,  it  can  always  be  described  by  a  single  first 
order  nonlinear  equation  which  can  have  its  exact  solution  presented  as  a 

quadrature  formula  0 

In  this  chapter  the  problem  of  the  second  order  (case  b))  under  a 
rectangular  trigger  will  be  analyzed  mainly  from  a  phase  plane  point  of  view, 
but  the  time  domain  will  be  used  whenever  convenient,  to  complement  such 

analysis. 

In  Chapter  5  we  will  extend  the  theory  to  treat  a  general  Eccles- 
Jordan  flipflop  (case  a))  and  also  the  internally  restricted  flipflop  (case  c)) 
which  will  be  treated  as  a  problem  in  its  own  right. 
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The  possibility  of  different  kinds  of  trigger  will  also  be  discussed 
in  Chapter  5,  and  the  methods  available  for  the  treatment  of  such  problems  will 
be  examined. 

Our  main  concern  is  the  phase  plane  analysis  of  the  second  order 
flipflop  under  a  rectangular  trigger,  but,  nevertheless,  we  shall  try  to 
approach  the  most  general  problem  as  much  as  possible . 

3-2  Phase  Plane  Analysis  of  the  Basic  Flipflop  Equation 
3.2.1  General  Remarks 

Before  discussing  our  specific  problem,  we  should  briefly  describe 
the  various  possible  types  of  singularities  on  the  phase  plane  for  a  system  of 
equations  such  as  (3-l),  where  P  and  Q  are  polynomials  with  no  common  factor. 

X  =  P(x,y)  1 

I  (3'1} 

y  =  Q(x,y) 

The  plane  (x,y)  is  called  the  phase  plane  (by  extension,  even  when 
y  f   x);  on  the  phase  plane  the  points  (xQ,  yQ)  such  that  y  =  £  =  0  are  called 
"singular  points"  or  "singularities"  of  the  system  (points  of  "velocity"  zero, 
"equilibrium  points,"  or  "states"  of  the  system)  [l,  h,    17,  22]. 

The  singularities  are  essentially  "point  paths"  in  the  phase  plane, 
and  therefore,  as  a  consequence  of  the  existence  and  uniqueness  theorem  a  path 
may  tend  to  a  singularity,  but  will  never  reach  it. 

The  nature  of  a  singularity  PQ  is  given  by  the  behavior  of  the  system 
in  an  arbitrarily  small  neighborhood  of  PQ.   This  behavior  is  the  same  as  that 
of  the  "system  of  the  first  approximation  about  P  ." 

If  PQ  has  coordinates  (xQ,  yQ),  the  system  of  the  first  approximation 
about  PQ  is  the  linear  system  obtained  by: 
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a)  making  a  transformation  of  variable 


X  =  x  -  x, 


y  =  y  -  yr 


>  so 


P(x,y)  =  P(x,y) 


Q(x,y)  -  Q(x,y) 


J 


(3.2] 


b )  taking 


A±1x  +   Algy 


A21X  +  A22y 


> 


(3.3) 


where  the  A   are  the  coefficients  of  the  terms  of  first  order  of  P  and  Q,  in 
the  order  indicated. 

The  solution  of  (3.3)  in  parametric  form  with  parameter  t  will  con- 
tain exponentials  of  i^t  and  \  t,   where  \Q  and  \p  are  the  two  characteristic 
frequencies  of  the  system,  i.e.,  the  two  eigenvalues  of  the  matrix  (\^)t    given 
by  the  solutions  of 


Au-X 


21 


A 


12 


A^-X 


=  0 


(3.M 


By  definition,  the  singularity  P0(x0,y0)  of  the  original  system  is 

called: 

(i)     Stable  node  if  \     and  \a   are  real  and  negative,  i.e.,  x(t)  and  y(t) 
contain  only  damped  exponentials . 
(ii)  Unstable  node  if  \  and  \     are  real  and  positive,  i.e.,  x(t)  and  y(t) 
contain  only  growing  exponentials . 
Ciii)  Saddle  point  if  \     and  \_  are  real  and  have  opposite  signs,  i.e.,  both 
x(t)  and  y(t)  include  one  growing  exponential. 
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(iv)   Stable  focus  if  \q   and  ^  are  complex  with  negative  real  part,  i.e., 
x(t)  and  y(t)  undergo  damped  oscillatory  motions, 
(v)  Unstable  focus  if  ^  and  ^  are  complex  with  positive  real  part,  i.e., 
x(t)  and  y(t)  undergo  growing  oscillatory  motions. 

The  study  of  singularities  is  important  mostly  because  if  a  system 
maintains  the  same  qualitative  properties  in  a  relatively  large  neighborhood  of 
a  singularity,  then  the  nature  of  the  singularity  will  give  us  considerable 
information  about  its  behavior  in  this  neighborhood. 

3-2.2  Existence  of  Singularities 

Taking  into  account  equations  (2.91)  and  the  group  (2.93)  through 
(2.96),  equation  (2.103)  was  obtained  from  (2.65)  and  (2. 77),  with  parameters 
as  described  in  Table  I;  for  convenience,  we  repeat  (2.103): 

V  +  V  +  %uX  =  dV  +  f^)t8(x  +  i)  -  5(x  -  £)]  +  ravs°  +  ns    (2.IO3) 

where  V  =  I,  II,  Illf    is  the  reglon  lndeXo 

We  define  a  new  variable  y  =  %f    and  express  (2.103)  as  a  system  of 
two  first  order  equations: 


0 

x  =  y 

0       D 
y  =  _ 

V 

Vu 

-   — ^  X 

a 

V 

^ 


b 

— — t-  v 

a  J 

V 

J 


> 


(3.5) 


where 


D  =  %  +  t(i)[b(X   +i)  -  6(x  _i)]  +mvg  +  ns 
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Notice  that  d  is  a  constant  in  each  region,  and  that  the  impulses 
occur  in  the  borders  between  regions  I  and  II  and  between  regions  II  and  III 
(this  is  a  direct  consequence  of  the  way  <p(x)  is  defined). 

Let  us  assume  that  s  is  a  rectangular  function  and  that  therefore  s 
is  a  pair  of  impulses : 


0, 


t  <  t, 


=  j  v,    te  <  t  <  te  +  t9  y 


\. 


0, 


tQ   +  TQ  <  t 


=  W(8(t  -  O  -  S(t  -  t   -  T  )} 


(3.6) 


(3^7) 


Convention  on  Impulse  Amplitude:   The  impulse  contained  in  §  has  its  strength 


proportional  to  y2/a,  but  occurs  at  a  point  where  both  "y"  and  "a"  are  discon- 
tinuous.  In  that  case  we  will  take  those  values  of  "y"  and  "a"  adjacent  to  the 
discontinuities  but  inside  region  II. 

It  is  easy  to  see  that  any  impulse  of  amplitude  A  contained  in  D  will 

Ao  .   o 
correspond  to  a  discontinuity  Ax  m  x: 

In  the  time  domain: 


Ao(t  -  t  )  -  AS  =  lUn(£(t  +  O  -  x(ta  -€)}-  — 


G-^0 


(3.8) 


In  the  phase -plane 


This  is  not  so  arbitrary  as  it  may  seem.   In  fact,  we  are  trying  to  analyze  a 
nonlinear  differential  equation  by  means  of  a  piecewise  linear  approximation. 
It  is  very  easy  to  show  that  any  other  choice  for  y  and  a  between  the 
respective  extremes  of  their  discontinuities  will  lead  to  the  existence  of 
crossed  paths  in  the  phase  plane,  which  certainly  will  not  be  a  portrait  of 
any  physical  system. 
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A5(x  ■■  xj  -  Ax  =  lim{x(x^  +  T))  -  x(x   -  Tl)}  =  - 


a' 


T)-0 


a 


a 


a 


(3-9) 


Therefore,  whenever  an  impulse  occurs,  it  can  be  taken  into  account 
by  assuming  new  initial  conditions  after  the  impulse  with  the  discontinuity 
in  x  described  by  (3.8)  or  (3.9)  and  x(t)  being  continuous. 

With  this  device  we  completely  eliminate  the  impulse  functions  from 
our  equation,  absorbing  them  into  suitable  discontinuities  of  x. 

Now  we  can  express  (2. 101+)  as  two  separate  cases: 


o 

x  =  y 


d 
o    vu 

y   a 
V 


r 


vu 


vu 
a 


y 


(3.10) 


Therefore,  the  phase-plane  equation  is 


■j    &,     -   c  x   b 
dy  _  Vu    vu  _   Vu 


dx 


a-  y 

VJ 


(3.11) 


where 


dVu  ~  dv  +  VnW'  whe^e 


if  trigger  is  OFF,  u  =  0 
if  trigger  is  ON,   u  =  1 


Obs,:  Notice  that  the  definition  of  d   has  a  meaning  only  in  the  case  of  a 
rectangular  trigger. 

The  singularities  x   of  (3.1l)  have  the  coordinates  (x  ,y   ),  where 
"the  x   are  the  roots  of 


d,   -  c  x  =  0 
vu    vu 


(3.12) 
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We  get 


x   =  Si       and   y   =  0  (3.13) 

Vu    C^  *VU 


where 


x   =  (x   ,0)  is  in  region  V,  and  u  is  the  trigger  index. 
Vu     Vu' 

If  x   is  not  in  region  V  for  a  particular  value  of  Vu,  then  we  refer 
Vu 
to  this  point  as  a  "virtual  singularity/'  since  it  has  an  effect  on  a  repre- 
sentative point  P  only  if  P  is  in  region  V,  but  ceases  to  have  any  effect  over 
P  if  P  leaves  region  V.   Since  x   itself  is  outside  region  V,  there  is  a 
neighborhood  of  x   in  which  P  is  not  affected  at  all  by  x^. 

Table  II  lists  all  possible  values  of  x^,  and  also  the  effect  of  the 
trigger;  the  important  special  case  of  state -symmetry  and  trigger  OFF  is 

presented  in  Table  III. 

It  is  not  superfluous  to  point  out  some  of  the  important  results 
presented  in  these  tables.   Let  us  consider  Tables  II. 1  and  II. 2  and  III. 

Notice  that  if  I  <  p7t  holds ,  then: 
u 

a)    *1[X  =  "  7  '  then  xnu  =  xiu' 

b)  XIIu  =  °>    thenXIH  =   "Xm^ 

c )  x  =  +  —  ,    then  xTT     =  x   TT    . 

d)  cp  +  un'W  is  such  that  X   is  virtual,  then  x^  is  also  virtual,  but  x^ 

is  real. 

e)  cp  +  un'W  is  such  that  X^  is  real,  then  x^  and  x^  are  also  real. 

f)  9  +  un'W  is  such  that  x^  is  virtual,  then  x^  is  also  virtual,  but  x^ 
is  real. 
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TABLE  II.   SINGULARITIES  (SEE  FIG,  9) 


TABLE  II. 1 


Singularity 


Iu 


T|_L         c 


in 


drr 
Up. 


"IIu       c 


IIu 


-  dnin 


"IIIu       c 


Illu 


Value 


un'W  -   Ejf  +  cp 


H 


un?W  +  cp 
I     -   Hp7t 


U-n'W  +  Ety  + 


Necessary  and  Sufficient 
Condition  to  be  Real 


I     <  p7(Ht  -  cp  -   un'W) 


P7  |cp  +  un'wj   <    |l     -   Hp7^ 


I     <  p7(H)Jr  +  cp  +  jjin'w) 

p* 


TABLE  II .2 

o 

•H 
-P 

•H 

o 

O 

Assumptions   Rectangular  trigger:   s,  = 

Asymmetric  flipflop:   s  =  ^w 

Symmetric  flipflop:    s  =  s   -  s   =  p,( 

uW 

V 

W2)  =  |iW 

Definitions 

dV|i  =  p(nn'W  +  i]/ 

0 

L+l^ 

+  9) 

V 

Cvn 

Cvu  =  Xu  -  HP^ 

-o- 

1 

Lo.j 

I 

-  1  +  R  G 
s  u 

1  -  a   Rs 

a   R 
o 

*  =  1  -  1  "  a  Rs 
a   Rq 

cp  as  defined  in 
Table  I 
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TABLE  III.   SINGULARITIES  (SEE  FIG.  9) 


Singularity 


10 


x 


1 10 


"IIIO 


Value 


■Hp  f- 


+Hp  f- 
0 


Necessary  and  Sufficient 
Condition  to  be  Real 


IQ  <  Hp7^ 


always  real 


I  <  Hp7^ 


Obs.:   This  is  a  special  case  of  Table  II  for  the  case 

of  cp  =  0  and  |i  =  0,  i.e.,  assuming  an  untriggered, 
state  symmetric  flipflop. 


g)  Therefore,  the  possibilities  are: 

i)  either  three  real  singularities 

ii)   or  one  of  the  extremes  (x   or  x^)  are  real,  the  other  extreme 

fx     or  xT  )  and  the  median  (x    )  are  virtual. 
v  IIIu     I|i  J--41 

However,  if  I  >  p7*, 
h)   Only  one  (any  of  the  three)  singularities  is  real,  the  other  two  being 

virtual. 

Figure  9  illustrates  these  conclusions. 
We  will  now  prove  that  whenever: 
i)  three  singularities  are  real  (case  g-i),  the  extremes  are  stable,  the 
median  is  unstable; 
ii)  only  one  singularity  is  real  (cases  g-ii  or  h),  it  is  stable. 


3.2.3  The  Nature  of  the  Singularities 

We  have  discussed  the  existence  of  singular  points,  and  established 
the  importance  of  several  relationships  among  parameters  upon  the  position  of 


JS      REAL     STABLE     NODE 
>JS     VIRTUAL   STABLE  NODE 
MS      REAL     SADDLE     PT 
>JS    VIRTUAL  SADDLE  PT. 
)ERLINE      SITUATION 


EUR 


ow  we  should 
Lons   (310)   can 


Ox 


(3.14) 

point 

P  and  Q  by 

ns   about 


(3.15) 


(3.16) 


(3.17) 


FIGURE  9:   EFFECT  OF  VARIATIONS  d^     ,    OF 

g(x)=  dw^-c^x,  ON  THE  POSITION 

(AND  EXISTENCE  !)  OF  SINGULARITIES. 
SEE  TABLES  1,111,11. 
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CURVE      POSITION 


t  in* 


m 

SYMBOLS 


o  MEANS     REAL     STABLE     NOOE 

®  MEANS    VIRTUAL   STABLE  NODE 

X  MEANS      REAL     SADDLE     PT 

X  MEANS    VIRTUAL  SADDLE   PT 

•  BORDERLINE      SITUATION 


o« 


FIGURE   9    EFFECT  OF   VARIATIONS  d^  ,    OF 
g(x)-d„M-c^«,ON  THE  POSITION 
(AND  EXISTENCE  !)  OF  SINGULARITIES 
SEE  TABLES  H.TJI.IS. 
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ILLINOIS  LIBRARY 


-52- 


singular  points,  and  therefore  upon  their  real  or  virtual  nature.   Now  we  should 
like  to  study  the  nature  of  these  singularities.  The  system  of  equations  (3I0) 
be  written  in  the  form  of  (3.1 ): 


can 


x  =  P(x,y) 


y  -  Q(x,y) 


(3.14) 


The  "system  of  the  first  approximation  about  the  singular  point 
^XVu'yvu)"  has  been  defined  as  being  the  system  formed  by  replacing  P  and  Q  by 
their  respective  first  terms  of  corresponding  Taylor  series  expansions  about 

vu 

Let  us  define  the  new  variables : 


X  =  x  -  x„.  ,•   Y=y-v 


The  system  of  equations  (3.14)  becomes 


(3.15) 


* m  hi****,,**  +  K^r^y* 


y"   vu^vu' 


Comparison  of  (3.1*0  with  (3.10)  shows  that: 

P(x,y)  =  y 


Q(x,y)  -554  -S±x   _S± 


av    av 


y 


(3.16) 


(3.17) 


Therefore 
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Px(VyVn}   =  °> 


VVyvn'  =  x 


c  "b 

V|i'-Vn'  av      '        V   VH,3V  "   "  a 


ScKu'^J 


(3.18) 


Let 


D  -  -  -»»       aM  =  _  _Vu. 

a..  n 


(3.19) 


Then    (3.l6)  becomes 


X  =  0   •   X  +  1    •    Y 


Y  =  D   •   X  +  E    •    Y 


(3o20) 


The   characteristic  equation  of  this   system  is 


■\  1 

D       E  -  X 


=  0,        i.e.,        \     -  EX.  -  D  =  0 


(3.21) 


So 


^p  =  l{E±^2  +  ^'} 


(3.22) 


i.e. 


=  \  IE  -  7e2  +  kB  | 


a 


V^E  + 


f 


Je2   +  k~D 


(3.23) 
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and  we  have  the  following  rules : 


fi)  \  and  \_  real  and  negative       =>  stable  node. 
v  '   a     p 

(ii)  \     and  \„  real  and  positive       =>  unstable  node. 
v  '       a  P 

(iii)  \  and  \„  real  and  opposite  signs  =>  saddle  point. 
v   '  a  p 

(iv)  \  and  X,  complex  conjugates 

LX  P 

with  negative  real  part  r>  stable  focus . 

(V)  A.  and  \n   complex  conjugates 
v  '       a  P 

with  positive  real  part  =>  unstable  focus. 
Application  of  these  rules  to  equation  (3-22)  or  (3.23)  is  straight- 
forward; replacement  of  D  and  E  by  their  expressions  in  terms  of  parameters 
(through  equations  (3-19)  and  Table  i)  will  produce  the  results  summarized  in 
Table  IV,  as  can  be  shown  by  the  following  analysis: 
Theorem  1.   If  x   or  x     exist  they  are  stable  nodes. 
Proof  :      From  definitions  of  E  and  D,  and  from  Table  I  for  regions  I  and  III 


we  have 


V 


av 


T  +T.+  T+T.  RG 

_  _  _i_oi_o_io_oja  <  Q  (3o2l+) 

(t.t  /t) 

V=I,III  V  1  O 


D 

av 


V=I,III     TiTo/x2 


1  +  R  G 

-^<0  (3.25) 


Therefore,  |Ve2  +  kl)    |  cannot  be  greater  than  |e|;  we  will  show  next  that 
E2  +  i+D  >  0:   in  fact,  from  the  expressions  above  for  E  and  D,  we  have 

f 

,2 


(T   +  T  .  +  T   +  T.  R  G  )   -  k(l    +    R  G  )T  T 

E2  +  l^D  =   x    Q1    ° 1Q  °  V a         (3.26) 


(t.t  hy 


1  o 


Consider  the  numerator  M  of  the  fraction  above: 
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TABLE  IV. 


A  SUMMARY  OF  THE  NATURE  OF  THE  SINGULARITIES 
IN  ALL  POSSIBLE  SITUATIONS  (SEE  FIG.  9) 


Possible  Situations 


Exists 


'iu 


'IIu 


'IIIu 


'Iu 


'IIu 


"III|a 


'iu 


'IIu 


'IIIU 


'in 


'IIu 


'IIIu 


•J 


Does  Not  Exist 


V 


V 


n/ 


V 


V 


V 


sj 


V 


s/ 


^ 


Nature 


Parameter  Condition 


Stable  Node 


c    >  0   and   -^  <  -  I 

llfi  cTT      7 

IIu 


or 


Stable  Node 


CIIu  <  °   and 


IIu 


7 


cnn  >  ° 


and 


IIu 


IIu 


7 


CIIu:     and   c^:'  +  7 


Stable  Node 


IIu 


or 


Stable  Node 


Saddle  Point 


c    <  0   and   -S± 

1J-n         c 

IIu 


<-i 


cTT  <  0 
IIu 


and 


Stable  Node 


IIu 


IIu 


< 
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M  =    (T       +   T     .    +   T       +  T.     R   G    f    -    Ml   +    R   G    )T   T 

v    i  oi  o  10  o  n  B  u     I  < 


■  Ti  +  2Vo  +  To  +  2(Ti  +  To)Toi  +  Toi  +  2(Ti  +  Toi  +  'o^ 


+   T2R2G2    -    Mr  .T       -    i+T  .T    R   G 

o   i(i  i  o  i  o  s   u. 


■  A  -  2TiTo  +  To  +  2(Ti  -  To)Toi  +  Toi  +  UVoi  -  2<T1  +  Toi  -  To)ToRl^ 


2   2  2 

O    1    |i 


2   2 


=   (\t  *  \  -  \f  -  SEiVo^ol  +  Ti  -  To)  +  \\%  +  Toi  +  Vol 


So 


M  =  [T  .  +  T.  -  T  (1  +  T  R.G  )]2  +  T   [T    +  k   T  ]  >  0       (3-2?) 

oi    i    o      o  l  u       oi  oi      o 

This  implies  that  E2  +  kV  >   0,  and  so  \Q  and  \     are  real  and  negative  in  both 
regions  I  and  III,  which  means  that  x   and  Xj^j  if  they  exist,  are  stable 
nodes,  which  was  to  be  shown. 

Notice  that  even  if  one  of  them  (or  both)  is  virtual,  its  action  upon 
its  corresponding  region  will  be  that  of  a  stable  node.   This  follows  from  the 
proof  of  Theorem  1  and  from  the  nature  of  the  coefficients  of  equation  (2.103) 

Theorem  2.   If  p7>l>  >  I   and  if  x    exists,  then  it  is  a  saddle  point. 


Parameters  as  defined  in  Table  II;  this  is  a  necessary  condition  for  bistable 
behavior  (but  not  sufficient),  for  otherwise,  either  x  or  x^  will  exist, 
but  not  both!   This  is  stated  in  an  equivalent  way  in  conclusion  h,  page  W. 
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Proof :    From  definitions  of  E  and  D,  and  from  Table  I  for  region  II  we  have; 


E  - 


I  in 
LII 


T  .  +  T 
1     O 


.  +  T   +  T   R  G   +  p7T 

1      o      lO  O  n    *' 


1  -  a 
a 


T.     R 

10     s 

"T"  +  R~ 
o-1 


T.T 
1  O 

T      *'     io 


<  0 


(3.28) 


D  = 


:IIp:  .    T^u  -  P7^ 


a 


II 


T  .T 
1  O 

—   +  P7Tio 


>  0 


(3.29) 


Therefore  \Je2   +  4d'  |  >  Ie^  \q  and  \p  are  then  real  with  opposite  signs,  and 
consequently  x    is  a  saddle  point  as  was  to  be  shown. 

The0rem  3°   If  ^  <   V   th^n  the  existing  x^  will  be  a  stable  node. 


Proof : 


Using  the  calculations  made  for  the  proof  of  Theorems  1  and  2 


we 


get,  if  region  II  is  considered: 


E     +  k-D  = 


TN 


"T  .T 
1    O 

.    T  10 


(3.30) 


where 


N  -  M  +  2(t      +  t        +  T      +  t.     R   G    )p7 
1  01  o  10    o    n  -^ 


1  -  a     10         s 
L     a        t     +Fj 


+     ipy 


1  -  a     10         s 
a        t     +  r 

0J 


(3.31) 


/T.T 


/T    T 

N'   =  M'   -   Ipyr.      +1    (  -i-2  +  pn 

(-T        IO  U  V      T  T    P/     iq 


(3.32) 


1 


No  bistable  behavior   is   possible! 
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where  M'  stands  for  the  sum  of  the  first  three  terms  of  N,  and  p7ty  <  I  was 
used.   Then 

T  .T 
N'  =  M«  +  I  -~^  (3-33) 

U.   T 

since  M  >  0,  M*  >  0  and  N1  >  0.   Therefore  N  >  0;  as  a  consequence,  \^   and  \ 
are  real  and  negative  and  if  x    exists  it  is  a  stable  node.   Of  course  the 
proof  of  Theorem  1  is  independent  of  EZ*  ,  and  therefore  it  is  valid  under  the 
present  hypothesis  of  p7^  <  I  .   Thus  the  proof  of  this  Theorem  3  is  complete. 


3 . 2 . k     Diagonalization  of  the  Characteristic  Matrix  of  the  System 

The  eigenvalues  and  corresponding  eigenvectors  of  the  characteristic 
matrix  are  given  by: 


(3-34) 


Therefore,  I.    =   \.k±Z    I  =  OL,    p. 

Let  k.  =  1  arbitrarily,  and  I     =  \^. 

Column  i  can  be  normalized  by  multiplying  it  by  a  factor  /-_ 

SIX.   - 

but  it  is  more  convenient   not  to  normalize  it. 
The  polar  matrix  of  r  is 


1  > 

\7  +  i 

i 


(3-35) 


Its  inverse  is 
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•1     1 


^  a 


(3.36) 


The  diagonal  matrix  of  r  is 


\  ° 


0     K 


and  the  diagonal  form  of  r  is  [28] 


(3.37) 


r  =  r  •  r,  •  r~- 

p    d    p 


(3.38) 


Therefore : 


0    1 


D    E 


(3.39) 


The  system  of  equations  (3.20)  can  be  expressed  in  matrix  form 


X 


o 


0    1 


D    E, 


(3.40) 


whose  solution  is: 


=  e 


0  1 
D  E 


r 


with 


T  =  t  -  t, 


X0  = 


x(o) 

Y(0) 


(3.41) 
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where,  by  definition,  given  a  diagonalizable  matrix  A,  with  polar  matrix  A^ 
and  corresponding  diagonal  matrix  A^, 


A   A      d  .  A-l 
e  =  A   '  e   "A 
P         P 


(3^2) 


and  for  any  diagonal  matrix  Q 


qx     0 


Q 


02 


if    Q 


*2 


(3.^3) 


Therefore,  from  (3-39),  (3-M),  (3-^),  ($M),   we  get: 


'X 


(3M) 


and  this  turns  out  to  be  % 


\-xa 


a  v     .  P 


V    v 

We  ' e      " e 


\  T  \0T 

YQl-ea     +eP 


\  T  \6T- 

WV     +V 


(3^5) 


In  another  form 
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xoX3  -  Yo  V  -Vt,  +  Yo  V 


-  -  Vg  '  Y0  ,  V   'Va  +  Y0  .   V 


(3.46) 


and  we  recall  the  definitions  of  X  and  Y: 


X  =  x  -  x 


VU>  Y=y  ' 


=  v.   since  y   =0 
V|_t  ■  /J '  JV\i 


(3.15) 


Equations  (3.^6)  furnish  the  integration  constants  A  and  B  of  the 
differential  equation  as  functions  of  natural  frequencies  and  initial  values 
Substitution  of  (3-15)  into  (3. 46)  yields  the  general  equation  for 
the  flipflop  trans ition,  if  we  keep  in  mind  that: 

(i)   Changing  from  one  region  into  another  changes  all 

parameters;  therefore  one  must  be  careful  in  calculating 
new  initial  conditions,  new  singularity  position  (which 
can  be  virtual!)  and  new  natural  frequencies. 
(ii)  The  same  thing  occurs  if  there  is  a  change  in  trigger 
level,  in  particular  when  it  is  turned  ON  or  OFF. 
(iii)   Impulses  produce  discontinuities  in  y  (and  therefore  in 
Y),  but  everything  else  remains  unchanged.   Remember 
that  the  effects  of  such  impulses  depend  on  a  }    which 
changes  from  one  region  to  another;  the  impulses  due 
to  effects  of  base  current  occur  when  a  is  changing 
from  a  to  a   and  also  when  it  is  changing  from  a   to 

-Li  JJ_ 

aIII'  in  each  case  we  wiH  take  the  average  of  the  two 
adjacent  values  of  a  . 
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The  results  are  presented  in  equation  (3-56),  section  3.3. 

However  equation  (3.K6)   is  also  a  very  practical  form;  discontinuities 
of  X  and  Y  can  be  calculated  from  equations  (3-15)  whenever  they  occur. 

Some  more  information  can  be  obtained  from  the  "canonical  system  of 
the  first  approximation";  this  is  defined  as: 


(3.47) 


where 


p    a 


(3-W) 


Of  course,  system  (3.V7)  is  equivalent  to  system  (3.UO)  under  the 
transformation  of  variables  (3.^8).   Such  a  system  clearly  has  a  unique 
singularity  at  the  origin.   Solution  of  (3-^7)  is: 


(3A9) 


i.e. 


A 


$  =  $  e 


X   T 

a 


> 


X   =  XQe 


V 


(3.50) 


from  which  we  find 
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t  -  t, 


1          $ 
- —  Smi  x, —  = 

ii,I 

\           <J> 

a          0 

p          0 

(3.51) 


and  therefore, 


(3.52) 


which  is  the  equation  for  the  phase  plane  trajectory  in  this  canonic  system. 

If  the  singularity  is  a  stable  node,  then  \     <   0  and  \     <  0,  and  we 
assume  either  |\J    |\R|  or  |\J    |\  |  (which  is  irrelevant,  since  the 


P 


a' 


ordering  of  X     and  \     is   arbitrary!), 
From   ( 3  o  52  )  we  get : 


(Xp/\^) 


froA     ) 

*.    P    a 


$ 


(3.53) 


and,  by  differentiation, 


ax  _ 

V 

d$ 

\ 

a 

(Xp/A^) 


$ 


$ 


v  P'  cr 


(3.54) 


Therefore,  (3.53)  defines  a  family  of  curves  of  the  parabolic  type  in  the  plane 
($,X),  for  x   =  f(*)  [1]. 

From  (3-51)  it  is  seen  that  the  direction  of  movement  of  a  repre- 
sentative point  P  over  any  of  these  trajectories  is  clearly  towards  the  origin. 
Figures  10  show  the  two  types  of  parabolic  curve  according  to  whether 


(a)   |Xa|<|X0|, 


04> 


(b)    IXalHXpl 

FIGURE  10:     CANONIC  SYSTEM  TRAJECTORIES  WHEN  THE 
SINGULARITY  IS  A  NODE.    IF  IT  IS  A  STABLE 
NODE,  i.e.,  IF  Xa  8  X£  ARE  NEGATIVE,  THEN 
THE  REPRESENTATIVE  POINT  P  MOVES  TOWARDS 
THE  ORIGIN  (SINGULARITY)  WITH  INCREASING  TIME. 
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|\  I  <  |\R|  or  |\ J  >  |\  I .   If  the  singularity  is  a  saddle  point,  then  we  take 
X .     <   0  <  X     (again  no  loss  of  generality)  and  inspection  of  (3.53)  and  (3.54) 
shows  that  (3-53)  represents  a  family  of  hyperbolic  curves  with  the  coordinate 
axes  as  asymptotes  [l]. 

From  (3.5l)  it  is  seen  that  the  direction  of  movement  of  the 
representative  point  P  along  any  of  these  trajectories  is  found  to  be  towards 
smaller  values  of  |$|  and  larger  values  of  |x|.   See  Figs.  11. 

We  wish  to  know  how  these  curves  transform  into  the  (X,  Y)  plane,  i.e., 
the  system  of  the  first  approximation,  which,  although  not  canonic^  also  has  a 
unique  singularity  at  the  origin. 

The  best  way  to  see  this  is  to  find  the  lines  in  the  (X, Y)  plane  that 
correspond  to  the  Oandx  axes.   From  (3.48)  (definition  of  $  and  X.)  we  get: 


$  axis :x=0    =>    Y  =  X  X 

a 


►  (3-55) 

X  axis  :   <J>  =  0    =>  Y  =  X0X 

P  J 

We  realize  that  in  the  case  of  a  node  with  |\  I  >  \Xn\ ,    i.e.,  with 

1  a'    '  p1 '     ' 

X     <   Xg  <  0,  the  trajectories  have  the  <£  axis  as  the  direction  of  their  axes  of 
symmetry,  and  that  they  are  tangent  to  the  X  axis  at  the  singularity  point 
(origin),  as  shown  in  Fig.  9h.   In  the  case  of  a  saddle  point,  the  curves  are 
asymptotic  to  both  axes. 

The  linear  transformation  conserve  all  these  properties,  and  also  the 
direction  of  motion  of  the  trajectories  with  increasing  times. 

Therefore  we  can  get  a  fairly  good  picture  of  the  family  of  trajec- 
tories in  the  (X, Y)  plane  if  we  know  them  in  the  ($,x)  plane.   The  results  are 
qualitatively  illustrated  in  Figs.  12  and  13. 

Figure  14  qualitatively  shows  some  portraits  of  the  original  system 
(x;y)  from  what  we  have  found  so  far  [l,  22]. 
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i>* 


(o)  \a>0,  \0<Q 


»  > 


(b)  Xa<0,  X^O 

COMMENT:    THE  REPRESENTATIVE  POINT  P  MOVES  IN  THE  DIRECTION  OF  HEAVY  ARROWS 
WITH  INCREASING  TIME. 

FIGURE  I!'       CANONIC  SYSTEM  TRAJECTORIES  WHEN  THE  SINGULARITY  ISA  SADDLE  POINT. 
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COMMENT: 


SLOPE  OF  a  a'  =  X, 


SLOPE  OF  bb'  =  X 


SLOPE  OF  mm' 


x„x 


aA/3 


xa+x/3 


OX 


FIGURE   12:   TRAJECTORIES  IN  THE  SYSTEM  OF  THE  FIRST  APPROXIMATION, 
CASE  OF  A   STABLE   NODE  WITH    |xJ>|Xo|   -(CORRESPONDS 
TO  FIGURE  10b). 
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COMMENTS  :      SLOPE  OF    aa'  =  Xa 
SLOPE  OF    bt>  =  Xg 

SLOPE   OF    mm=XaXft 

THE    CURVES   ARE    HYPERBOLIC,   ALL    ASYMPTOTIC 
TO    LINES    aa  abb', WITH    MAXIMA  AND  MINIMA 
ON    LINE    mm';    ALL    CURVES     CROSS  THE  X-AXIS 
PERPENDICULARLY.      HEAVY    ARROWS    INDICATE 
DIRECTION    OF    MOVEMENT    OF     P    WITH   INCREASING 
TIME. 

FIGURE    13:    TRAJECTORIES  IN   THE     SYSTEM    OF 

THE    FIRST    APPROXIMATION,    CASE    OF 
A   SADDLE     POINT,  WITH         fxa>  0 


V 


mm~^ — ■>* 


Ox 


(c)  NEGATIVE   BIAS  ON  gU)=   dj.  <0 

but   di/x/cffM<l/y«  cnM  <0 


Ox 


Ox 


cter  of 
ponding 


le 


:o  a 


could 


>ond 

or  to  a 


her  from 


dg 
dx 


(f)  STABLE    SYSTEM:     cv    :   >0: 


FIGURE     14:  g(x)=    dWft-c^x;  dV/i  a    cu^ 
AS    IN    TABLE    I 


refore, 


(o)  BALLANCED  g(x)   dj    =0;    but    c^       <0 


Ox 


Ox 


O    « 


O    « 


(b)  POSITIVE    BIAS    ON    g(x):  d^      >0     but     |dK    /Cjj| 

<"'r.      o-    <o 


I 

A 

I 
A 

K 

A 

Y 

y 

A 

Y 

V  O  « 


O    « 


Ox 


WMONOSTABLE      SYSTEM  IN    I:    iL    /Cj.     >  +l/y 


(e)  MONOSTABLE      SYSTEM   IN  E:     dnM/cII/i  <_l/y 


FIGURE    14   glx)=   d^-c^x,  d„M  8    cu^ 
AS    IN    TABLE    I 
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It  is  worth  reminding  that,  due  to  the  piecevise  linear  character  of 
our  complete  system  in  (x,y),  in  each  region  the  solution  of  the  corresponding 
system  of  the  first  approximation  is  exactly  equal  to  the  solution  of  the 
original  system. 

3-2.5   Comments  on  Figs .  l^t- 

a)  Balanced  g(x):   d    =  0;  this  situation  would  usually  correspond  to  a 
state  symmetric  (d^  =  0)  untriggered  (a  =  0)  system.   However,  we  could 
have  this  same  situation  if  d^  /  0,  nW  =  -d^,  and  u  =  1,  since 

diiu  =  dno  +  MnW>  in  this  case  dm  =  °- 

b)  and  c)  Nonzero  bias  on  g(x):   a    /  0;  this  situation  would  correspond 

either  to  a  state  asymmetric  (d^  =  0)  untriggered  (u  =  0)  system,  or  to  a 
partially  triggered  (no  matter  what  kind  of  state  symmetry)  system 

(din  =  dno  +  nW  *  °)- 

d)  and  e)  Monostable  systems  in  either  region  I  or  III; 


lip 


>  — ;  either  from 


v    •  Cl1^ 

a  large  bias  in  an  untriggered  system,  or  from  an  adequate  trigger  in  any 
potentially  bistable  system. ' 
f )  Stable  system:   c    >  0;  the  system  is  not  a  flipflop. 

General  Comments 

1)  Notice  the  relationship  between  d    (bias  or  triggering  or  both, 

dIIl  =  dII0  +  nW^  and  the  Position  of  the  singularities:  the  value  of  ^ 

dx 

at  a  singular  point  and  its  nature: 


n 


Potentially  bistable  system  is  any  system  which  would  exhibit  bistable 
behavior  if  adequately  biased,  i.e.,  any  system  for  which  cTT„  <  0:  therefore, 
the  system  of  Fig.  13f  is  not  potentially  bistable.        II0 
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dg 
dx 


>  0  ->  stable  singularity  (in  our  case,  stable  node) 


dg 

-^  <  0  =>  unstable  singularity  (in  our  case,  saddle  point) 

Of  course,  the  exact  nature  of  the  singularity  also  depends  on  a  and  b 

V        VJ! 

2)  The  position  of  the  singularity,  if  it  is  real,  must  be  inside  its  corre- 
sponding region;  otherwise,  it  is  virtual. 

3)  The  singularity  nature  and  position  is  just  a  translation  of  the  coefficient 
of  the  differential  equation;  this  is  one  way  to  interpret  the  "action"  of  a 
singularity  over  its  corresponding  region  no  matter  where  it  happens  to  be 
and  justifies  calling  it  "virtual"  if  it  happens  to  be  outside  its  region  of 
influence:   it  does  not  really  exist,  but  it  would  exist  if  the  parameters 
CV|_i  aM  dvu.  Were  the  same  throughout  the  phase  plane  as  they  are  inside  the 
proper  region.   In  this  sense,  this  would-be  singularity  adequately  trans- 
lates the  coefficients  of  the  differential  equation  inside  its  proper  region. 

h)     We  have  used  a  coded  tag  to  describe  the  singularities  in  this  figure: 


SN  =  Stable  Node 
SP  =  Saddle  Point 


R  =  Real 
V  =  Virtual 


The  last  symbol  of  the  tags  on  the  singularities  is  one  of  I,  II  or  III,  and 
designates  the  proper  (corresponding)  region  of  the  singularity;  so,  if  the 
position  of  the  singularity  agrees  with  this  symbol,  it  is  real;  otherwise 
it  is  virtual.   The  symbol  R  or  V  is  therefore  redundant,  but  has  been  used 
for  clarity. 
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3„3  Trajectory  Equations 

Solution  to  equation  (2.103)  is 


o 
v   =  x 

JVU     V|i 


(3.56) 


where 

a)  V  refers  to  the  region:  V  =  I,  II,  III 


u  refers  to  the  trigger:   u  ==  0,  1=>  0  trigger  OFF,  ON 


h )     A  A     c   are  parameters  corresponding  to  a  given  value  of  the  pair 

J      avu'  pvu'  vu    * 


V|a. 

c)  \  \         are  the  natural  frequencies  corresponding  to  a  value  of  the  pair  V\x, 
'       avp/   pv^i 

d)  t   is  the  instant  of  time  when  the  pair  of  indices  assumes  a  new  value. 
;   Vu 

e)  The  most  common  sequence  of  index  pair  values  is  V^  =  10,  II,  III,  III1, 

IIIO. 

f )  Let  us  use,  in  general,  a  prime  to  mean: 

a1,  V,  c',  d!,  differential  equation  parameter  values  in  the 
previous  V\i   condition 

solution  parameter  values  in  the  previous  V|i 
condition 

natural  frequencies  in  the  previous  Vu  condition 
values  of  t,  x,  y  at  the  end  of  previous  Vu 
condition 


A*,  A'   C! 
or   p 


t',  x',  y' 


Actually  c   =  x   is  the  abscissa  of  singular  point  X  --and  here  X^  has 

nothing  toV^do  wife  X  =  X  -  x   used  in  the  previous  section. 

yyx 


■73- 


Unprimed  symbols  will  refer  to  parameters  and  variables  throughout  the 
present  V|i  conditions,  which,  whenever  necessary,  will  be  indicated. 


tQ  =  f 


xQ  =  x' 


U 

yo =  y +  a 


with 


U  =  magnitude  of  any  impulse  occurring  during  the  transition 
from  the  previous  value  of  vn  to  the  present  value;  it  is 
zero  if  no  impulse  occurs  at  this  point.   It  may  be  a 
function  of  yQ  or  y\    as  in  Table  V. 

t0>    x(y  y0  =  values  of  t,  x,  y  at  the  beginning  of  present 
V|a  condition. 


h)  A  _  +  (xo  -  x*\  -  yo 


(xn    -    X*)\       -    V, 

A  _   __0 J  a       Jo 

0  "        \_  -  X 
P   a 


C  -  £  =  x- 

c 


where 


-d 


(x*,0)  -  (-, 0)  =  location  of  singularity  corresponding  to 
present  V(i  conditions. 
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IABLE V.   IMPULSE  VALUES  FOR  CHANGES  IN  Vu 


V\x   Condition  Change 


d 
o 

■H 

-P 

•H 

w 

s 

a) 
U 

M 

H 


CI 
O 
•H 
-P 
•H 
W 

d 

EH 

H 

O 
P 


From 


10 


II 


III 


III1 


IIIO 


III1 


III 


II 


To 


II 


III 


nil 


IIIO 


nil 


in 


n 


10 


Magnitude  of  Impulse  U 


10 


+2p  —  |W 


-Hp7 


'io  2 


T   J0 


+Hp7  -=-  y* 


•2P  ^  |„, 


-2P  -f  |W| 


+Hp7 


'io  2 

y^ 


T  J0 


T 

-Hp7  —  y1 


+2p-^  |V| 


Comments 


W  >  0 


y  is  initial  value 
or  y  in  region  II 


y'  is  final  value  of 
y  in  region  II 


W  >  0 


W  <  0 


y^  is  the  initial 

0 
value  of  y  in  region 

II 


y'  is  the  final  value 
of  y  in  region  II 


W  <  0 


H  is  the  symmetry  factor;  H  = 


1  for  asymmetric  flipflop. 

2  for  symmetric  flipflop. 
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i) 


x    .1 

a       2a 


-b  + 


k  fb  -  ^~^\ 


j )  And  for  complet 


eness  we  repeat  equation  (2.IO3)  in  the  new  notati 


on 


00   _o 

ax  +  hx  +   ex  =  d.   for 


a  given  Vfi 


Now  equation  (3.56)  becomes 


x  =  A  e  a    °    a   3    0J 

+  A  e         +  x 


>s 


a 


> 


y  =  A  X   e  a    °  +  a  A  e  3    ° 
a  a         3  0 


^ 


In  the  phase  plane,  from  (k)  and  (h); 


(3.57) 


1        (x  -  X*).\   -  y 


0      ^       T^0  -  x*)\ 


y- 


a 


^0 


(3.58) 


°r,  in  another  form: 


eVp(t-to)  .  />  -  x*)\, 


y  1  a 


x~x**a  -  yl 


(x  -  x*)\p  -  y/p 


(*0  -  x*;xp  -  yc 


(3.59) 


Notice  that  t  and  X   are  normalized  and  have  no  dimensions! 
Obs.:   Equations  (3.58)  and  (3.59)  are  phase  plane  trajectory  equations,  and  the 
calculation  of  transition  time  is  hased  on  them.   Also  based  on  them  is  the  phase 
Plane  graphical  construction  which  simplifies  not  only  transition  time  calcula- 


tions but  also  the  analysis  of  waveforms 
and  trigger. 


and  design  optimization  of  both  circuit 
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In  this  analysis  the  trigger  is  assumed  rectangular,  so  the  term 
ns  =  unW  is  imbedded  into  d,*  the  impulse  terms  ms  -  (jmW[s(t  -  tQ) 
_  6(t  _  t  _  T  )]  and  also  f(x)[g(x  +  h   -  o(x  -  ^)]  are  both  considered  in 
g),   with  their  magnitudes  represented  by  U.   This  is  summarized  in  Table  V. 

Therefore,  expansion  of  g)  yields  for  this  special  case: 


Discontinuity 


y0  =   y1  -  hQ 


y0  =    y?  +  V 


yn  =    y*  +  hcc 


y0  =  i1  -  h 


,2 


Transition 

Through  Line 

I  -*  II 

x  = 

1 

7 

II  -*  III 

x  = 

1 
+  7 

III  -*  II 

X  = 

1 
+  7 

II  ->  I 

X  = 

1 
"  7 

(3.60a) 
(3.60b) 
(3.60c) 
(3-60d) 


where 


.)     I  - 


^II7 


2a 


II 


b)  m       and  a       as    in  Table   I 

c)  H  is  called  the  "symmetry  factor" 


H  = 


1  for  the  asymmetric  flipflop 

2  for  the  symmetric  flipflop 


We  can  write  equations  (3.60)  in  more  explicit  forms,  and  also  obtain 
the  inverse  functions : 
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Discontinuity 


M     y0  =  ^{^i  +  hh>'  -  1} 


> 


(ii)      y'    =    (1  +  iyj 


0^0 


(i)        7n  -    (1  +  iy') 


y 


(ii)  y'  =  ^{wrr^  -  i}J 


(i) 


=  ~H-  >/i  -  W1'  + 


yo  ~  2£ 


r    +  1 


> 


(ii)         y'     =      (1     -     i^)^ 


(i)         yQ   =    (1    -    iy«)y« 

(ii)  y«  -^{-N^n^r  +1j 


J 


Transition    Through  Lines 


I  -»  II 


II  -*  III 


III  -  II 


II  -  I 


x  =  -  i    (3.6la) 


1 
x  =  +  — 

7 


1 
x  =  +  - 

7 


1 

x  =  -  — 

7 


(3.6lb) 


(3.61c ) 


(3.6ld) 


Notice  that  in  cases  (a)  and  (b)  we  have  the  I  to  III  transition,  i.e., 
7  >   0  through  the  transition,  whereas  in  cases  (c)  and  (d)  we  have  the  III  to  I 
transition,  i.e.,  y  <  0  through  the  transition.   The  signs  of  the  square  roots 
are  selected  from  this  physical  consideration. 


3«^   Separatrices 

We  are  using  the  word  separatrix  in  a  somewhat  looser  sense  than  its 
strictly  mathematical  meaning  [17].   We  shall  call  a  "separatrix"  any  phase 
Plane  trajectory  which  divides  the  phase  plane  portrait  of  the  system  into 
qualitatively  distinct  families  of  trajectories. 

The  concept  of  "qualitatively  distinct"  is  purposely  vague;  this  means 
bat  a  separatrix  will  be  so  with  respect  to  some  stated  qualitative  distinction 
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between  the  two  families  of  trajectories  in  which  it  divides  the  phase  plane 
portrait  of  the  system. 

In  order  to  define  a  separatrix,  all  we  need  to  do  is  to  find  one  of 
its  points  Q,  let  Q  =  (x  ,y  ),  and  enter  xQ  and  yQ  instead  of  xQ  and  yQ  into 
equation  (3-58)  or  (3. 59).   Also  x*  must  be  known  (x*  as  defined  in  (3.56I1)). 

We  have  a  special  interest  in  defining  two  types  of  separatrix  over 
the  whole  portrait : 

(i)  "Transition  Separatrix"  which  divides  the  portrait  into  two 
families  of  trajectories:  those  that  cross  through  region  II  and  those  that  do 
not.  This  separatrix,  as  illustrated  in  Figs.  15  and  l6,  is  composed  of  four 
branches:   A,  B;  C,  D;  A,  B;  C,  D. 

These  lines  will  be  given  special  names : 

a)  ABCD:   (I  to  III)  transition  separatrix. 

b)  ABCD:   (ill  to  l)   transition  separatrix. 

c)  ABBA:   end-point  separatrix. 

d)  DCCD:   initial-point  separatrix. 

Since 

a)  ABCD  separates  the  lines  which  are  trajectories  from  region  I  to  region  III 
from  the  lines  which  are  not  so. 

b)  ABCD  separates  the  lines  which  are  trajectories  from  region  III  to  region  I 
from  the  lines  which  are  not  so. 

c)  ABBA  divides  the  portrait  into  two  sets  of  trajectories:  those  that 
terminate  in  x   and  those  that  terminate  in  XTII  • 

d)  DCCD  divides  the  portrait  into  two  sets  of  trajectories:   those  that 
originate  in  region  I  and  those  that  originate  in  region  III. 

(ii)   "Critical  Separatrix,"  whose  importance  will  be  later  explained, 
divides  the  portrait  into  two  families  of  trajectories:   those  that  lie  partly 
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in  region  II  and  those  that  lie  exclusively  in  region  I  or  III.  Fig.  15 
illustrates  the  separatrices  and  the  notation  used  in  this  section.   Notice 
the  following: 

a)  The  Transition  Separatrix  Determining  Points:  Obviously,  in  region  II, 
the  separatrix  itself  consists  of  the  two  asymptotes  whose  equations  are 
(see  Fig.  lk): 


yp  =  y(x-x»)i   let   yc  =  yp(+  ±),    y-  =  ^(.  I) 


> 


where 


x  is  in  region  II 

(x*,0)  is  the  singularity  corresponding  to  region  II 

\^   and  *£  are  the  natural  frequencies  inside  region  II 


We  find 


rA»  yD>      yA'        yD 


respectively  from 


TV  V       ^B'       ^C 


through  equation  (3.59): 


(a,ii),   (b,i),   (c,ii),     (d,i), 


after  the  replacements; 


yA"*y'  yjD-v  ^A^y''    ^-y0 
yB-.v  yc"y''  ^v    yc^y 


(3.62) 
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We  get  the  following  equations 


-   (1  +  hB)yB  (3.63a) 


yn  =  (1  +  hnhr  (3  =  63b) 


D        v  JC,,7C 


y-A  =  (l  +  £yg)ys  (3-630 


yB  =  (1  +  ^c)yc  (3o63d) 


b)  The  Critical  Separatrix  Determining  Points:   Of  course  those  trajectories 
which  cross  the  x  axis  at  points  x  =  -  j   and  x  =  ~   are  the  two  branches  of 
the  critical  separatrix. 

c)  The  Separatrices :   The  two  separatrices,  as  illustrated  in  Fig.  15,  are 
given  by  the  general  trajectory  equation  (3.58)  or  (3-59)  (they  are  entirely 
equivalent).   Since  here  we  are  more  concerned  with  geometrical  properties, 
we  prefer  the  latter,  which  we  repeat  here  as  equation  (3-6U),  leaving  out 
the  exponential  of  time,  and  replacing  (xQ,y0)  by  (xQ,yQ).   So  the 
separatrices  are  given  by: 

f  (x-x«V-yfo   f  (x-x*)^-y|XP  } 

t(*Q  -  **>\rTJ    -  lTxQ  -  x*)\p  -  yj 

and  Tables  VI  and  VII  furnish  the  values  of  the  parameters. 

Besides  the  separatrices,  three  more  lines  at  each  stable  node  are 
important  in  the  qualitative  study  of  trajectories  in  the  phase  plane.   These 
lines  are:   (by  definition)  assuming  \Q   <  \^   (see  Fig.  12): 
a)  The  "tangent  line":   y  -  \p(x  -  x*),  which  is  a  tangent  to  all  trajectories 

at  point  (x*,0). 
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TABLE  VI.   PARAMETERS  OF  THE  SEPARATRIX  EQUATION  (3.64) 


Calculate  by 

-- 

Equation   (3.63) 

Table   II 

Equation   (3„56i) 
and  Table   I 

Parameter 

XQ 

(depend  on  fi) 

X* 

X. 

1 

i  =  a,  p 

Branch 

Value 

Transition 
Separatrix 

Branch  A 

1 
7 

yA 

in 

iIji 

Branch  D 

1 

+  — 

7 

yD 

XIIIH 

N.IIIH 

Branch  A 

1 

7 

yA 

XIII(i 

xinm 

Branch  D 

1 
7 

yD 

XT 

N.in 

Critical 
Separatrix 

Branch  E 

1 

7 

0 

I|I 

ii|j. 

Branch  E            +  — 

i    7 

0 

XIII^ 

^illlH 

— 

Notice  that  the  index  (i  merely  indicates  if  the  trigger  is  ON  or  OFF;  the 
actual  direction  and  effects  of  the  trigger  (if  ON)  must  be  computed  through 
Tables  I  and  II.  B 


TABLE  VII.   RESULTS  OF  EQUATIONS  (3.62)  FOR  THE  BRANCHES  OF  THE 
TRANSITION  SEPARATRIX  INSIDE  REGION  II 


Branch 

Equations    (3.62) 

Furnish 

Transition 
Separatrix 

Branch  B_,    B 

ycrV(x  "  xn} 

yB  = 

yB  = 

ya(+  £) 

Branch  C.,    C 

yp  =  y(x  -  xn) 

yC  = 

V+  7) 

Tq     = 

y  (--) 

y^\    7) 

Notice  that  both  y.'s  and  X.'s,  I  =  a,    &,    depend  on  ,1,  so  all  y's  depend  on  p, 
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b)  The  "direction  line":   y  =  \  (x  -  x*),  which  does  not  cross  any  trajectory 
except  at  point  (x*,0),  where  it  crosses  all  trajectories. 

c)  The  "max-min  line":  v  =  .   p.   (x  -  x*)  which  crosses  all  trajectories  at 

a    p 
(x*,0)  and  also  every  trajectory  at  its  point  of  maximum  or  minimum  value 

of  y(x). 

d)  The  line  described  in  (t )  also  applies  to  the  case  of  saddle  points,  except 
that,  at  the  singularity  (x*,0)  it  crosses  only  the  asymptotes,  since 

the  other  trajectories  do  not  pass  through  the  singularity. 


3.5  Trajectories  and  the  Action  of  zhe   Trigger 

The  siinplest  trigger  is  the  rectangular  current  trigger.   In  fact  the 
action  of  a  trigger  of  a  different  waveform  will,  in  general,  differ  in  detail, 
hut  not  in  principle  from  the  action  of  the  rectangular  'rigger.   For  this  reason 
we  have  considered  it  important  enough  to  he  the  oasis  of  this  work,  and,  in 
this  section,  we  will  discuss  its  action  in  a  qualitative  manner. 

3.5.1  Turning  the  Trigger  ON  and  Possibility  of  "Under-Triggering" 

..gures  16  show  various  possibilities  of  trigger  action  upon  the 
phase  plane  portrait—singularities  and  separatrices--and  also  the  corresponding 

initial  value  of  y. 

Let  us  assume  that  the  flipflop  is  in  stable  state  I,  i.e.,  the 
representative  point  P  is  at  stable  node  I,  when  a  rectangular  trigger  of 
amplitude  W  is  applied  (a  positive  trigger). 

The  immediate  effect  of  the  trigger  upon  the  phase  plane  portrait 
to  shift  the  stable  nodes  to  points  respectively  £kx  and  ^IT1   "to  the  right,  and 
to  shift  the  saddle  point  by  £x       to  the  left,  where  Acy  =  xyl  -  xvQ,  as  given 
by  Table  II. 
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At  the  same  time  P  is  shifted  upwards  to  a  point  at  (x  ,y  ),  i.e., 
the  x  coordinate  does  not  change,  but  y  goes  to  an  initial  value 

y0  =  ^W  (3-65) 

Suppose  that  the  singularities  have  not  been  shifted  out  of  their 
proper  regions  (Fig.  l6b ) .   The  transition  separatrix  in  region  I  approaches  the 
critical  separatrix  (whose  shape  changes!),  thus  reducing  the  minimum  initial 
Value  yQmin  °f  y  necessary  fo?   a  complete  transition  to  occur. 

Let  X^  denote  the  singularity  (x  ,0),  and  let  y^^  be  the  least 
initial  ordinate  at  xIQ,  for  which  a  I  to  III  transition  occurs. 

There  will  be  two  possibilities: 

a^  y0  <  yOmin'  then  P  wil1  follow  some  parabolic  trajectory  and  tend  towards 
xj]_'   Here  there  are  yet  three  possibilities. 

lo   y0  <  :\/ScI°   P  wil1  not  cross  the  x  axis,  moving  towards 
X  ,  in  an  overdamped  manner. 

2°   y0  =  ^a^l1   P  w111  not  cross  the  x  axis,  following  a 

straight  path  towards  X  ,  approaching  it  in  a  critically 
damped  way. 

3»   yQ  >  ^a^Xj:   P  will  cross  the  x  axis,  moving  towards 
X   in  an  underdamped  way. 

In  case  3  we  might  still  distinguish  the  two  possibilities  of  P  going 
through  region  II  (entering  it  under  point  y  )  or  not. 

Si 

y0  >  yQmin'  then  P  wil1  als0  follow  a  parabolic  path  towards  X   in  the 

X[J. 

underdamped  way,  as  in  (a. 3)  above,  but  it  enters  region  II  above  point  y 
before  reaching  the  x  axis.   Then  a  transition  occurs.   This  effect  will  be 
called  "under-triggering"  (Fig,  19 ). 
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One  might  then  expect  that  the  minimum  value  W  .   of  W  necessary  to 

mm  J 

cause  a  transition  would  be  slightly  less  than  the  value  W     of  W  necessarv 

crit  J 

to  shift  X^  and  X^  out  of  their  proper  regions  (and  into  virtual  existence!). 

But  this  is  usually  not  so,  as  it  will  he  proved  at  the  end  of  this  section 
that,  under  certain  (usual)  conditions,  of  all  possibilities  mentioned  above 
only  (a.l)  occurs,  all  others  being  impossible  for  the  type  of  circuit  under 
consideration.   Therefore,  under-triggering  usually  does  not  occur  for  this  type 
of  circuit.   However,  it  is  a  possibility,  especially  in  a  general  equation, 
whose  coefficients  were  related  in  some  other  manner.   We  will  later  study  this 
effect  in  some  detail. 

3-5.2  Virtual  Singularities  and  the  Trajectory 

We  establish  that  W^  =  W^^  =   WQ;  suppose  that  W  >  WQ;  then  the 
portrait  becomes  as  in  Fig.  l6d.   There  is  only  one  real  singularity,  and  this 
is  the  stable  node  Xim;  the  end-point  separatrix  vanishes*  (since  all  lines  must 
now  terminate  at  Xmi);  the  critical  separatrix  and  part  of  the  initial-point 
separatrix  also  vanish,  and  the  remaining  part  of  the  initial-point  separatrix 
loses  its  meaning. 

From  its  initial  position  at  (xJ0,y0),  P  "sees"  only  the  virtual 

stable  node  Xn  (somewhere  in  region  II  or  III)  and  starts  to  move  on  a  parabolic 

path  towards  it.   However,  before  reaching  x__,  P  crosses  the  line  x  =  -  -  where 

11  7 

y  suffers  a  discontinuity  (-^y),  '  and  enters  region  II  where  now  it  "sees"  only 
the  virtual  saddle  point  Xm  (somewhere  in  region  i)  and  changes  its  trajectory 
into  a  hyperbolic  path  asymptotic  to  the  line  y  =  \in«(x  -  Xj   )  (a  remaining 
part  of  the  initial-point  separatrix,  and  here  we  see  why  it  is  meaningless); 


t  o 

See  footnote  on  page  90. 
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f inally,  P  crosses  the  line  x  =  +  -  ,  where  y  suffers  a  discontinuity  (+\j), 
and  enters  region  III. 

Once  in  region  III,  P  will  "see"  only  the  real  stable  node  x  m, 
towards  which  it  will  start  to  move  in  a  parabolic  path  (in  an  overdamped 
manner,  as  will  be  shown). 

3.5.3  Turning  the  Trigger  OFF 

The  next  event  with  the  trajectory  of  P  is  the  turning  OFF  of  the 

trigger. 

Rigorously,  the  trigger  may  be  turned  OFF  as  soon  as  P  has  progressed 
far  enough  into  region  II,  i.e.,  to  a  point  where,  after  the  negative  jump  of  y 
caused  by  the  trigger  turning  OFF,  P  finds  itself  at  side  III  of  the  end-point 
separatrix  (of  the  [X   =  0  system,  of  course). 

On  the  other  extreme,  we  could  leave  the  trigger  ON  for  an  indefinite 

amount  of  time. 

We  are  interested  in  establishing  a  criterion  with  which  to  judge  the 

adequacy  of  the  trigger  duration. 

3.5.U  Discussing  Trigger  Duration 

Assuming  the  trigger  is  sufficiently  long  to  produce  a  transition, 
we  recognize  five  possibilities  for  the  trigger  (see  Figs.  17  and  18): 
(i)  Too  short:   if  it  is  turned  OFF  while  P  is  still  in 
region  II. 
(ii)   Short:   if  it  is  turned  OFF  long  before  P  reaches  the 

line  x  =  x    ,  but  after  P  is  in  region  III  (Figs.  17a,  b 

and  l8a,  b ) . 

P 

Discontinuities  caused  by  impulses  +  Ay     £ —  respectively. 

dx 
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(iii)  Fair:   if  it  is  turned  OFF  approximately  as  P  crosses 
the  line  x  =  x    . 
(iv)   Long:   if  it  is  turned  OFF  long  after  P  crosses  the 
line  x  =  xIII0,  but  before  it  approaches  X 
(Figs.  17c,  d  and  l8c,  d). 
(v)  Too  long:   if  it  is  turned  OFF  after  P  is  already 
close  to  Xim. 
Of  course,  these  definitions  can  be  formalized  and  made  exact:   so,  a 
sufficient  rectangular  trigger  starting  at  tQ  =  0  and  having  duration  T  is  said 
to  be  : 

(i)  Too  short:   if  x(t)  <  -  - 

(ii)  Short:   if  +  -  <  x(t)  <  x     -  e 

7    v  J  IIIO    0- 

(iii)  Fair:   if  ^  -  eQ_   <  x(T)<  x^  +   eQ+ 
(iv)  Long:   if  x^,  +  6Q+  <_x(T)  <  xJnjL  -   ^ 

(v)  Too  long:   if  x     -  e  <  x(t) 
where  eQ_,  eQ+,  ^  are  positive  numbers  such  that  all  the  inequalities  above 
can  be  satisfied. 

No  matter  which  case  occurs,  y  will  suffer  a  discontinuity  equal  to 
(-yQ);  from  the  point  (x(t),  y(T)),  P  will  jump  to  (x(t),  y(T)  -  yQ),  and  then 
it  will  move  towards  XIII0  by  some  parabolic  path,  thus  completing  the 
transition. 

NOTE:   Of  course,  P  never  reaches  *1II0>   but  given  an  arbitrary  neighborhood 

eino  of  xino^  there  is  a  time  Tnio(eiiio)  such  that  for  any  t:>  Tino(eiiio^ 

P  is  in  €III0.  We  can  therefore,  arbitrarily  select  a  neighborhood  N     of 
XIII0,  and'  by  definl"tion,  say  that  a  transition  is  complete  after  P  enters  N 
for  the  last  time. 
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3.5.5  The  Concept  of  "Optimum  Trigger  Duration" 
Call: 


and 


PT_  the  point  (x(t),  y(T>) 


> 


PT+  the  point  (x(t),  y(T)  •■  yn) 


(3.66) 


etc.,  be  vectors  with  the  same  coordinates  as  the  points 
designated  by  the  symbols  under  the  arrows. 
By  definition,  let 


and  let  PT_,  PT+, 


RIII0(W,T)  =  |PT+  -  Xmoi 


(3.67) 


Given  a  rectangular  trigger  of  amplitude  W,  we  define  "optimum  duration 
T*"  as  that  value  of  T  for  which  Riiio(t)  is  a  minimum.   That  is: 


Rmo(W,T*)  =  min  RIIIQ(W,T) 


(3.68) 


Since  for  each  value  of  W  there  is  a  corresponding  value  T*,  we  con- 
clude that,  for  a  given  circuit,  T*  is  a  function  of  W: 


T*  =  e(w) 


(3.69) 


This  function  9(w)  is  a  characteristic  of  the  complete  circuit,  i.e.,  flipflop 
and  triggering  circuit  together. 

Notice  however  that  the  definition  of  T*  is  somewhat  arbitrary,  and 
probably  there  is  no  ideal  criterion  on  which  to  base  a  definition  of  optimum 
duration.   It  certainly  depends  under  what  criterion  we  would  like  to  have  the 
transition  optimized. 
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A  more  practical  definition  is  as  follows:  given  a  rectangular 
trigger  of  amplitude  W,  we  define  "optimum  duration  T*M  as  the  value  of  T  for 


which  x(T)  =  x    .   That  is 


X(T*}  =  XIIIO  (3.70) 


The  discussion  presented  in  this  subsection  can  be  applied  to  a 
III  to  I  transition,  if  we  make  the  necessary  (and  obvious)  changes. 

The  criterion  for  the  definition  of  T*  expressed  by  equation  (3.70) 
is  the  most  useful,  and  will  be  used  throughout  this  report.   So,  unless  other- 
wise specified,  the  expression  of  "optimum  duration"  or  the  symbol  T*  will 
imply  "as  defined  by  equation  (3.70)." 

3-5.6  Possibility  of  "Back-Triggering" 

We  have  said  that  after  the  trigger  is  turned  OFF,  if  it  is  suffi- 
ciently long,  P  (whose  y  coordinate  has  suffered  a  negative  discontinuity  equal 
to  (-y0))  "sees"  only  Xmo  towards  which  it  moves  by  some  parabolic  path. 

However,  we  must  ask  ourselves  if  this  is  always  true.   There  seems 
to  be  nothing  in  the  nature  of  the  equation  to  warrant  this  assumption.   The 
objection  is:   "The  position  (x(T-),  y(T-))  of  P  at  the  moment  of  turning  the 
trigger  OFF  might  be  such  that  the  new  position  of  P,  (x(T+),  y(T+))  (where 
y(t+)  =  y(T-)  -  yQ)  would  be  under  the  branch  A  of  the  transition  separatrix, 
and  therefore  P  would  return  to  X^,    rather  than  going  to  X    » "  This  effect 
will  be  called  "back- triggering. " 

In  fact,  the  possibility  of  back-triggering  is  small,  unless  the  coef- 
ficients of  the  differential  equation  were  not  related  by  the  circuit  parameters 
(representing  some  other  analogous  type  of  bistable  device). 
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We  will  prove  that,  under  certain  (usual)  conditions,  hack-triggering 
is  not  possible  for  the  circuits  under  consideration.   The  conditions  that 
make  back-triggering  impossible  are  the  same  that  make  under-triggering 
impossible. 

Actually,  these  two  characteristics  are  closely  related.   We  will 
presently  discuss  these  effects  in  some  detail,  explain  their  interrelation,  and 
find  the  conditions  that  make  them  impossible  to  occur. 


3.5.7  Trajectory  After  the  Trigger  is  Turned  OFF 

Let  us  assume  the  trigger  duration  is  sufficient  and  that  no  back- 
triggering  occurs.   Figures  17  and  18  illustrate  the  four  cases  as  considered 
below: 

(i)  x(T)  <  xIII(),  y(T+)  >  0 
(ii)  x(T)  <  xIII0,  y(T+)  <  0 
(iii)  x(T)  >  xIIIQ,  y(T+)  >  0 
(iv)  x(T)  <  xIII0,  y(T+)  <  0 
Figure  19  illustrates  the  three  possibilities  in  the  case  of  optimum 
triggering:  x(T)  =  xIII(y 
(i)   y(T+)  >  0 
(ii)  y(T+)  =  0 
(iii)  y(T+)  <  0 

We  should  point  out  that  x(t)  is  quite  arbitrary,  since  we  have 
absolute  control  of  the  trigger  duration  T;  but,  for  a  given  circuit,  y(T+)  is 


'  We  point  out  again  that  back-triggering  refers  only  to  the  case  of  sufficient 
trigger,  i.e.,  there  must  be  an  interval  of  time  Tmin  <  T  <  Tmax  for  which  a 
normal  transition  would  occur.   T  <  T^  means  insufficient  trigger  duration; 
T    would  be  imposed  by  back-triggering. 

m  a  v 


max 
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a  function  of  x(t)  given  by  the  phase  plane  portraits  of  the  differential 
equations  (with  trigger  ON  and  OFF).   This  means  that,  for  a  given  differential 
equation,  selection  of  T  =  T*  will  lead  to  a  certain  value  of  y(T*+)  for  which 
y(T*+)  <  0;  the  equation  and  also  the  trigger  amplitude  W  will  determine  which 
relationship  holds . 

It  would  he  interesting  to  know  how  the  coefficients  of  the  equation 
and  the  circuit  parameters,  as  well  as  the  trigger  amplitude  affect  the  curve 
y(T)  versus  x(t).   It  would  also  be  useful  to  know  how  T*  and  y(T*+)  depend  on 
W  for  a  given  equation. 

In  the  following  chapters  these  questions  will  be  considered. 

3.6  Under-Triggering  and  Back-Triggering 

Let  us  analyze  the  possibilities  of  under-triggering  and  back- 
triggering  for 

(i)  a  given  differential  equation  of  type  (2.103)  with 
coefficients  defined  in  the  three  regions 
(ii)  ignoring,  in  this  section,  the  relationships  established 
in  Table  I,  but  still  assuming 
(iii)  a  rectangular  trigger  (i.e.,  (3*6)  and  (3.7)  hold) 
and  that,  as  before 
(iv)  the  function  f(x)  =  f(y),  the  magnitude  of  the  impulses 
occurring  at  x  =  +  -  ,  are: 
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r 


e-O 


,   at  x  =  +  i 


f(y)  =  ^ 


undefined  for  any  x  f  +  - 


(3.71) 


with 


m. 


m.. 


I     _   _  r11  or  —  ,  according  to  whether  the  flipflop  is 
II   2a-[I    a^   asymmetric  or  symmetric. 


3.6.1  Under-Triggering 

Suppose  a  trigger  of  amplitude  W  <  WQ  is  applied  to  the  system;  assume 
P  was  at  X   =  (xI0,0)  before  the  trigger  was  turned  ON. 

P  jumps  to  a  point  PQ  =  (x1Q,yQ)>   i*  ?  ls  above  branch  A  of  the  transi" 
tion  separatrix  of  the  triggered  system,  then  a  transition  will  occur,  with  P 
going  to  X  m,  rather  than  to  X^  (under-triggering).   Fig.  19  illustrates 
this  effect.   Therefore  a  transition  will  occur  if  and  only  if 


y^  >  yQ 


10 


where  v    is  the  ordinate  of  branch  A  of  the  transition  separatrix  of  the 
Jal0 

triggered  system  (u  =  l)  at  x  =  x^. 

The  problem  is  to  find  the  value  Wmin  of  W  that,  for  a  given  system, 
will  cause  the  point  (xI0>y0)  to  be  on  branch  A  of  the  transition  separatrix 
of  the  triggered  system  (i.e.,  u  =  l). 


We  could  keep  the  expression  f(y)  =  -&y     at  any  other  point  x  {   +  -,  taking 
I  =  Z„    ,  Vu  indicating  region  and  state  of  trigger;  but  this  would  be  irrel- 
evant, since  the  function  is  multiplied  by  zero  at  these  points  anyway. 
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As  could  be  expected,  we  will  see  that  it  is  not  possible  to  solve 
this  problem  analytically,  but  only  by  a  graphical  or  iterative  numerical 
procedure.   In  fact,  we  have  a  set  of  formulae,  repeated  here  for  convenience, 
that  can  be  used  to  solve  this  problem: 

From  (3.13)  and  the  Obs .  at  the  end  of  section  3.2.2, 


d     d 
v^   cvu   cvu    V 


(3.72) 


from  which  we  get 


"10       c 


10 


I  n     w 

11  ""'  CI1    "  CI1 


■II 


"110        c 


1 10 


d-r-T 

x         =  +  ■ w 

III       clzl       clzl 


>      (3*73) 


0  <  w  <:  w„  =>  x™  <  xTn  <     i       0  <  w  <  wQ  =>  0  >  xin  >  -  - 


'0  10  II  7 


From   (3.62),    using  the  notation  \  and  X         as  before 


pv^x 


■  1 


yBl  =:  "Xalll"(7"fxni)- 


W  <  WQ  =>  yB1  >  0 


(3-7*0 


yAl  =  yBl^1  +  ^yBl^  Co.0-x*jr        jA1-    jb1 


ilearly       yA1   >  yT 


(3.75) 


with 


2  an ' 


H  = 


1  for  the  asymmetric  flipflop 

2  for  the  symmetric  flipflop 


Also,  from  equation  (2.103)  itself  (see  section  3.5.1): 
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mI 


(3.76) 


And  finally,  the  general  trajectory  equation,  expressed  by  (3.59) 
be  rewritten  as : 


can 


(xi  -  X*K  -  yj-fo  .  f(\  -  **frp  -  y . A 


(x.  -  x*)\  -  y 


(xj  "  x*^  -  57- 


(3.77) 


where  (x^y^  and  (x^y^)  are  any  two  points  over  the  same  trajectory,  in  the 
same  region,  and  x*  is  the  abcissa  of  the  singularity  and  \     and  \     are  the 

OL  (3 

natural  frequencies  corresponding  to  that  region;  so,  in  (3.77)  let: 


(x.,y.)  =  (x  ,y ■  );     x*  =  x 


"> 


II 


UyyP  =    (-  pyAb  \a  and  ^  -  ^  and  ^ 


►       (3.78) 


J 


Considering   (3-73),    (3-7*0,    (3-75),    (3-76)  and   (3.78),    (3.77)  yields 


mT  /d     +  nW       d     N 

ai  aI1V    cn         cia 


>, 


all 


dII  +  nWN 
Xalll  (7  +  — c 


I 


III 


,        d        +  nW 


x,        d     +  nWN 
+  \  ,„   f  -  + 


all  \  7  c 


II 


y 


=  < 


m  /d     +  nW       d     N 


alii  \y 


1  __  dn  +  nW 


111 


311 


£\ 


1  +  dn  +  nW\ 

in   y 


alii  V  7  c 


/,        d     +  nW^ 


^ 


(3-79) 


which  can  be  numerically  solved  for  W. 
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If  a  real  and  positive  root  W  x  can  be  found  for  (3-79)j  such  that 

W   <  W  ,  then  we  have  under-triggering  whenever  the  trigger  duration  W  satisfies 
rl    0' 


w  „  <  w  <  wA  (3.8o) 

rl       0 


In  this  case, 


W  .   =  W  T  (3.8l) 

mm    rl 


If  no  such  W   exists,  then  no  under-triggering  can  occur,  and  we 


rl 

have : 


w  .  =  wn  (3-82) 

mm    0 


Expansion  of  the  numerator  and  denominator  in  the  fractions  appearing 
in  (3=79)  furnish: 


K 


W  +  V      I"11  .   f    KBl'W  +  V      I  PI1  (3„83) 


(w2  +  p   w  +  ^  kv2  +  Pplw  +  QpI 


al 

,H+>,  v    m    P    Q    ¥         M  .  P   9  Q~T  as  in  Table  VIII. 1. 
with  KaI,  MaI,  Par  Har    ^pI,  iipr  rpT,    *pX 

An  entirely  analogous  analysis  could  be  made  for  the  case  of  a  III  to 
I  transitions   reverse  the  signs  of  all  coordinates  and  trigger,  and  change 
subscripts  I  to  III,  A  to  A,  B  to  B. 

3„6.2  Back-Triggering 

We  will  consider  only  the  possibility  of  back-triggering  in  the  case 
of  too-long  trigger,  i.e.,  P  is  assumed  to  be  at  Xim  =  (xnil,0)  when  the 
trigger  is  turned  OFF. 
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TABLE  VIII. 1.   PARAMETERS  OF  (3.83)  AS  FUNCTIONS  OF  THE  PARAMETERS  OF  (3.79) 


K 


al 


M 


a  I 


"al 


<k 


K 


PI 


M 


PI 


.ai  +  XqI1  cI1, 


&f\ 


n 


alii  c 


III 


AL 


\ 


aI1Vcn     cio, 

An 

sS +  X(X11  °^L 


in 


tj\ 


a  I II 


2iV  TT.  (■ \  +  J2 

alii  \  7       c 


_^  +  Wii: 
III/        Xalllcll 


-  -  1 


III 


n2i\ 


7   +   c 


II 


a  III 


III. 


^     TT,      A    +    — 

alll  \y       c 


-  1 


IZL 


aI1  V     cn. 


aT  +     pil  c 


11. 


i  U 


J  J 


n \2 

III- 


k 


pilvgn   <w 


in 


pi 


X2^\ 


^Siii  r"aInV 


i  +  !n 


a. 


'III. 


Will 
XaIIlcll 


Q, 


ill 


Pi 


nl\ 


a  III 


cm 


a  TT1 A  +  !s.>) 

aII1v   W 


- 1 


♦»m(l-£ 


ii. 
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P  jumps  to  a  point  PQ  =  (xIII:L*  -7Q)f    if  p0  is  under  branch  A  of  the 
transition  separatrix  of  the  untriggered  system  (u.  =  0),  a  transition  will 
occur,  with  P  returning  to  X  Q,  rather  than  going  to  XIII0  (hack-triggering) . 
Figure  20  illustrates  this  effect.   Therefore,  a  transition  will  occur  if  and 
only  if : 

-y0<ya01    °r    y0>lya01l  {3'Qk) 

where  v    is  the  ordinate  of  branch  A  of  the  transition  separatrix  of  the 
'aOl 

untriggered  system  (u  =  0)  at  x  =  Xj-q^- 

The  problem  is  to  find  the  value  W    of  W  that,  for  a  given  system, 
will  cause  the  point  (xmi>  -yQ)  to  be  on  branch  A  of  the  transition  separatrix 
of  the  untriggered  system  (i.e.,  u.  =   0). 

As  for  the  case  of  under-triggering,  we  will  see  that  this  problem 
cannot  be  solved  analytically  either,  but  only  by  a  graphical  or  numerical 
iterative  procedure . 

A  set  of  formulas  similar  to  the  one  used  in  the  case  of  under- 
triggering  can  be  used  to  solve  this  problem. 

We  had  obtained  (3.72)  which  we  repeat  here  for  convenience: 


x   =_YH  =  _v_+  n-iL-w  (3.72) 

v^   cvu   cvu     V 


We  get . 
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110   c 


II 
110 


kIII0   c 


"III 
IIIO 


xTT1  =  -£-  +  -2-  w 
111   CII1   CII1 


0  <  W  <  W_  =>  0  >  xT„  >  -  - 

0         III     7 


III      n   TT 


1111    ciin    ciin 


0  <  W  <  W„  =>  x  <x 

0  IIIO  III1 


>    (3^85) 


<  (2  "  7)xmo 


J 


Again  from   (3.62 ),    retaining    the   notation  \         and  \ 

avu     pvu ' 


yB0   XaII0'^7+xii0^ 


(3*86) 


yA0  =  yBO'^-^Io) 


(3.87) 


From  equation  (2.103)  itself: 


m 


■y^  = 


in 


0     a 


w 


III 


(3.88) 


(we  have  implicitly  taken 
true ) . 


LIII  "  ai  and  mHl  =  mi>  but  they  may  be  not  strictly 


And,  for  convenience,  we  repeat  the  general  trajectory  equation  (3.59), 
in  the  form  (3.77); 


(xi  -  **K  -_iAa    r(xi    x*)xe  ■■  yil^ 

T^^TX-Ty/ 


"(XJ  - **h  -  yJ 


(3.77) 


with  (xi,y±),  (x  ,y   ),  x*,  \  and  \  as  before 


J   J 
Now  let 


a      3 
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■> 


(x^y.)  =  (xim,-y0): 


(x^y^)  =  (+  -,  y^0)i? 


X*    =    X 


IIIO 


(3.89) 


xa  and  \  =  Xaliio  and  Vno 


J 


Considering   (3.85),    (3-86),    (3.87),    (3.88)   and   (3.89),    (3.77)  yields: 


r 


m 


J III_ 

01110  V      cno/l 


,  d^^  +  nW       dTTT 

III  tt  (     m  m 

—  W  +  X  TTT.    — - 

a allxOV        c. 


"s 


'III1 


'iiia 


alllO 


iX 


i  +  d^ 


alio  V  7       c 


IICL 


A     dm  ^ 
aino^     cmoy/ 


V 


v. 


m 


-i 


III 
aHI 


W  +   X 


'dIII  +   nW        dIII 


PIIIOV       c  cina 


X 


pino 


> 


A 


i    dII 

+ 


"alio  \7       c 


IIC 


ix 


fl  ,  !ll^    .  ! 


bII°  V  "  cna 


+  x 


A 


III 


PIIIO  \7        cina 


(3.90) 


which  can  be   numerically  solved  for  W. 

If  a  real  and  positive   root  W     n  can  he  found  for   (3-90),    then  we 
have  hack-triggering  whenever  the  trigger  duration  W  satisfies 


W  TTT  <  W 
rill 


(3.91) 


In  this   case, 


W  =  W   TTT 

max         rill 


(3.92) 


and   in  order  to  obtain  a  permanent  transition,   we  must  have 


W  <  W 


max 


(3.93) 


If  no  such  W     exists,  then  no  value  of  W  will  cause  back-triggering, 
rill 
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Expansion  of  the  numerator  and  denominator  in  the  fractions  appearing 
in  (3.90)  furnish: 


L(W 


'  +  painw  +  %aiJ  '  V  +  VnW  -  QPII/  (3'94) 


with  K 
a 


in'  Maiir  paiir  %air  Kpin>  Mpm>  ppnr  Qpm  as  in  Table  VIII°2° 

An  entirely  analogous  analysis  could  be  made  for  the  case  of  a  III  to 
I  transition:   reverse  the  signs  of  all  coordinates  and  trigger,  and  change 
subscripts  III  to  I,  I  to  A,  B  to  B. 

3.6-3  Discussion 

It  is  clear  from  inspection  of  either  Fig.  20  or  equation  (3.79)  and 
of  either  Fig.  21  or  equation  (3.90)  that  the  necessary  and  sufficient  conditions 
for  the  impossibility  of  occurrence  of 

under-triggering : 


m_ 


'W    '    IVll   >a~W^        °<W<W0  (3.95) 


back-triggering  s 


X 


HI 

1    '     ' Villi    >a— :W'         0<W<WQ  (3.96) 


where 


Vv  =  xvi^  '  xvo  (3.97) 


Expansion  of  (3-97)  yields 
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TABLE  VIII. 2.      PARAMETERS  OF    (3 .94)   AS  FUNCTIONS   OF  THE  PARAMETERS  OF    (3-90) 
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y  'O 


.  BOTH  ARE  DEFINED  ONLY 
J  IF   t0*t(x|0»' 


FIGURE  21:    DEFINITION  OF  TIME  INTERVALS  OVER  A  TRAJECTORY  IN  THF  flMfi 
DOMAIN,  RELATED  TO  THF  PHASE  PLANE.    (8EeTaBLE». 


n    w/v    - 

dv 

n    w  - 

,     Si  "  cvo 

Si        Si 

So 

Si 

v    SoSi 
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Ax  (3-98) 


Notice  that  if  either  of  conditions  (3-95)  and  (3-96)  is  not  met, 
then  there  will  exist  real  positive  solutions  W^  or  Wrm  to  (3-79)  and  (3-90) 

respectively. 

We  clearly  see  that  conditions  (3-95)  and  (3-96)  are  formally  the 

same . 

a)  If  we  make  the  further  assumption  that  the  coefficients  of  equation  (2.103) 
except  d  are  the  same  in  regions  I  and  III  (a  realistic  assumption!),  then 
the  only  difference  between  them  is  the  effect  of  different  values  of  the 
coefficients  with  the  index  u  (neglecting  variation  of  capacitances). 

b)  If  we  make  the  further  assumption  that  the  coefficients  of  equation  (2.103) 
are  invariant  with  u  (i.e.,  trigger  circuit  is  fixed!),  then  the  two 
conditions  are  identical. 

This  shows  a  close  relationship  between  these  two  effects,  i.e.,  they 
have  the  same  intrinsic  nature  and  origin  from  the  circuit  parameter  point  of 
view,  and  if  any  distinction  exists  between  them,  it  is  due  to  the  fact  that 
the  circuit  itself  is  not  the  same  in  each  case  (unless  the  two  conditions 

above  hold ! ) . 

One  further  assumption  leads  us  to  an  interesting  point:   assume  that 


G  =  0,    u  =  0,  1 


(3-99) 


i.e.,  consider  the  case  of  a  trigger  circuit  that  closely  approximates  a  true 
current  source.   Then  we  have  the 

Theorem  k.      If  equation  (2.103 )  describes  either  an  asymmetric  or  a  symmetric 
flipflop  (coefficients  as  in  Table  i),  and  if  the  trigger  circuit  satisfies 
(3.99)  above,  then  no  under-triggering  or  back- triggering  can  occur. 
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Proof.    Conditions  a  and  b  can  be  expressed  as  a  single  condition: 


m 


V  '  n>i  (3.100) 


since 


cvu  =  1   for   v  =  ^  XII>  n  =  o,  1 


and  with 


Ka  -   -(b  ♦  Vb2  -  Ua' )  i 


T.T 
1  O 

a  = 


T2 


T.  +  T  .  +  T 
b  =        01 


T 

L    = 

=  2p 

T  . 

io 

T 

n 

=  2p 

R 
s 

R 
o 

as  in  Table  I,  for  regions  I  and  III.   Therefore,  (3. 100)  becomes,  after 
expansion  and  simplification: 


2T.      < 

10 


{(T1  +  Toi  +  V  +  ^Ti+Toi+To)2-UTiTo'}  T  (3.101) 


This  becomes 


1  <  2    ^  +  7)  \(°  +  cr  +  r)  +  ^(c  +  cr  +  r)2   -  kcr   [  =   f(c,r)      (3.102) 
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where 


C  R„ 

c  =  — ,   and   r-T 

o  i 


Now,  let  us  make  a  change  of  variable,  as  follows 


r  =  r 


c  =  Sr 


so  that  f(c,r)  becomes  g(&,r): 

g(&,r)  =  \   (1  +  \)   |(6r  +  5r2  +  r)  +  V(5r  +  &r2  +  r)2  -  4sr2 '  j 


or 


g 


(6,r)  -  |  (r  +  l)  {(1  +  5  +  Br)  +  J(l   +  6  +  Br)2  -  ^'}  1   (3.103) 


Condition  (3.102)  has  become: 

1  <  g(B,r)  (3.1CA) 

Let  us  find  the  partial  derivative  of  g(8,r)  with  respect  to  r: 

k^ai("(1+6+5r)  W(l+6+Br)2   -kb     +   (rflhlu    f,2{l  +  \PT\} 
&r  2  l_v  I        V(i+  6  + Br)     -  h&   J 

(3.105) 


Then 


'  The  discriminant  always  positive,  since  (l  +  6  +  &r)  -  *+&  -  (l  -  B  +  Br) 
+  U&2r  >  0  so  g(6,r)  is  a  real  positive  number. 
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%^>o, 


for  all   6  >  0,  r  >  0 


(3.106) 


Therefore 


lim  g(6,r)  <  g(s,r) 

r-O 

rr  >  0 

s-*o 

8  >  0 

(3.107) 


But 


lim  g(8,r) 

r^O 
8-^0 


=  1 


(3.108) 


Sub 


stitution  of  (3.108)  in  (3.IO7)  yields 


1  <  g(8,r) 


r  >  0 
8  >  0 


(3.109) 


and  therefore  condition  (3.100)  holds,  and  the  theorem  is  proved. 

Comments : 

1)  Notice  that  the  theorem  is  very  strong,  in  the  sense  that  condition  (3.109) 
is  strong;  in  fact,  1  is  not  only  a  lower  bound  for  g(5,r),  but  it  is  an 
infimum  of  g(s,r),  i.e.,  its  greatest  lower  bound!   So,  anything  at  all 
upsetting  the  assumptions  made,  may  invalidate  the  theorem,  and  in  this 
case  either  under-triggering  or  back-triggering,  or  both,  at  least  in 
principle,  could  occur! 

2)  Notice  also  that  (3.IO6)  is  weak.   In  fact, 


dg(6,r) 

or 


>  lim  hi^il 

5>0     8^ 


(3.110) 
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So,  in  fairness,  we  should  point  out  that  the  possibility  of  occurrence  of 
under-triggering  or  back-triggering  in  practice  is  not  so  great,  since  any 
"reasonable"  value  at  all  of  c  and  r  ought  to  satisfy  condition  (3.100) 
with  a  good  margin. 

3.7  Summary 

In  this  chapter  we  have  analyzed  the  phase  plane  characteristics  of 
the  basic  flipflop  equation  in  the  case  of  a  rectangular  trigger. 

Conditions  related  to  the  existence  and  nature  of  singularities  were 
discussed  and  three  theorems  were  proved  with  respect  to  this  point. 

Some  properties  of  the  system  were  also  established  by  diagonalization 
of  the  characteristic  matrix  of  the  system  (in  effect,  considering  two  possi- 
bilities, respectively  for  the  two  possible  types  of  singularity). 

A  general  trajectory  equation  was  established,  and  the  geometry  of 
the  separatrices  was  discussed,  as  well  as  the  action  of  a  trigger  upon  the 
system  phase  plane  portrait,  with  special  attention  to  the  effects  of  turning 
the  trigger  ON  and  OFF.   Here  the  possibilities  of  under-triggering  and  back- 
triggering  were  discussed,  and  a  theorem  on  the  conditions  for  such  a  possibility 
to  exist  was  proved  for  an  important  special  case. 


h.      ANALYSIS  AND  DESIGN  TECHNIQUES 

U.l   Introduction 

At  this  point  we  would  like  to  utilize  this  information  we  have  about 
the  bistable  system  represented  by  equation  (2.103)  to  the  purpose  of  developing 
some  analysis  and  design  techniques .   Specifically,  we  wo  ild  like  to  find 
effective  methods  to  solve  the  following  problem:   given  a  flipflop,  its 
loading  and  triggering  circuits,  find  the  transition  wave  forms  of: 
(i)  base  currents  and  voltages 
(ii)   collector  currents  and  voltages 

We  would  then  be  able  to  find  the  optimum  trigger  duration .  Further- 
more, knowledge  of  the  base  and  collector  currents  and  voltages  as  functions  of 
time  would  help  to  improve  the  design  of  the  overall  system  [l,  k] . 

And  finally,  if  the  influence  of  the  circuit  parameters  upon  wave  form 
characteristics  is  known,  we  would  have  a  means  of  optimizing  the  design  of  the 
system  towards  approaching  some  transition  requirements  [7], 

Lastly,  if  the  transistors  are  given  (t  is  given)  and  also  the  trigger, 
but  if  the  circuit  (resistors  and  capacitors)  is  abritrary,  then  the  lower  limit 
in  transition  time  can  be  calculated,  and  a  convenient  figure  of  merit  for 
transistors  describing  their  performance  in  switching  circuits  can  be  defined, 
and  would  certainly  be  useful  in  the  selection  of  transistors  for  switching 
applications  (see  Chapter  5)  [10,  20] . 

^.2  Definitions  of  Time  Intervals 


We  have  divided  the  range  of  the  variable  x  into  three  parts,  which 
were  called  regions  I,  II  and  III,   Remember  that  x  is  a  normalized  form  of 


By  effective  we  mean:   "a  sensible  compromise  between  accuracy  and  ease  of 
application." 
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the  base-to-base  voltage.   Let  us  consider  the  variables  wfc,  which  are  normalized 
forms  of  the  collector  currents  (of  transistors  1  and  2  f or  k  =  1,  2,  respec- 
tively), as  given  by  equation  (2.93). 

It  is  clear  that,  as  long  as  x  is  in  region  I,  w1  =  0  and  wg  =  1; 
whenever  x  is  in  region  III,  this  situation  is  reversed,  with  v^  ■  1  and  wg  =  0; 
in  both  these  cases  one  of  the  transistors  is  cutt  off,  and  the  other  is  con- 
ducting a  fixed  current,  i.e.,  the  transistors  are  inactive;  they  are  active 
only  when  x  is  in  region  II. 

Def .  k.l.      So,  in  a  I  to  III  transition,  from  the  point  of  view  of  collector 
currents,  the  time  during  which  x  is  crossing  region  I,  from  x(tQ)  towards  -  -, 
is  really  a  delay.   It  will  be  called  "delay  time"  and  designated  by  1^.        If 
the  circuit  were  settled,  x(tQ)  =  x^  otherwise,  it  may  happen  that  x(tQ)  /  x^. 

1  ,    1   . 
Def.  4.2.   The  time  interval  when  x  is  in  region  II,  going  from  -  -  to  +  -  ,  is 

characterized  by  activity  of  the  transistors,  and  variation  of  collector  cur- 
rents.  It  will  be  called  "active  time"  and  designated  by  1^. 
Def.  4.3.   And  finally,  for  the  time  interval  when  x  is  already  in  region  III, 
from  +  -  until  final  settling  in  a  neighborhood  Nm0  of  *1110,    the  collector 
currents  are  constant  (having  reached  their  final  values)  and  the  transistors 
are  again  inactive.   It  will  be  called  "complementary  time"  and  designated  by 

T  . 
C 

However,  many  things  can  happen  while  x  is  in  region  III. 
Def.  k.k.      The  time  interval  in  which  x  goes  from  +  -  to  xIII0,  with  trigger  ON, 
is  called  "balanced  time"  and  designated  Tg. 

Def.  k.5.      The  time  interval  between  the  moment  the  trigger  is  turned  OFF  and 
the  moment  x  settles  inside  Nm0  is  called  "settling  time"  and  designated  Tg. 


This  time  interval  is  often  called  "discrimination  time." 
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Def'  k°6°      The  time  lnterval  between  the  moment  x  reaches  x^  (or  wo  ad  reach 
XIIIO  lf  the  trigger  were  kept  ON)  and  the  moment  the  trigger  is  turned  OFF 
is  called  "trigger  excess  overtime"  and  designated  T^.   We  will  make  the 
convention  of  using  negative  values  of  T^  if  the  trigger  is  turned  OFF  before 
x  reaches  x^,  i.e.,  given  the  function  t(x),  then  T^  =  TQ  -  t(xmo)  where 
t(xlII0)  is  calculated  assuming  a  sufficiently  long  trigger  (or  measured!)  and 
T0  is  the  trigger  duration. 


Def°  h°7°      The  time  interval  between  the  instant  x  reaches  +  i  and  the  instant 
the  trigger  is  turned  OFF  is  called  "trigger  overtime"  or  simply  "overtime/1 
and  designated  T^  i.e.,  TQV  =  TQ    -  t(+  i). 

With  our  criterion  for  optimum  trigger  duration  T*  (see  equation 
(3.70))  we  will  have: 

Te  =  t^xnio)  (^.ia) 

and  if  TQ   =   T|,  then  from  Def.  k,6, 

TE0  =  °  C^.lb) 

And  we  define  "optimum  overtime"  T*  : 

Tov  =  T?  "  *C+7)  =  TR  (^ic) 

-ef°  k°Q°      In  the  case  of  too  long  a  trigger,  we  define  the  "long  settling 

time,  '  TLg,  as  the  time  interval  between  the  moment  x  reaches  +  -   and  the  moment 

it  settles  inside  a  given  neighborhood  N     of  x 

III1     1111° 

Def\_ih9*   Finally,  we  call  "l-III  transition  time,"  TTR,  the  total  time  interval 
between  the  moment  the  trigger  is  turned  ON  and  the  moment  x  crosses  the  line 
XIIIO° 
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The  following  relations  are  obvious  from  the  definitions: 

TTR  =  TD  +  TA  +  V    "   X(t0)  =  XI0  (U-2a) 

Tc  -  TQV  +  Ts  (U.2b) 

T0V  "  TB  +  TE0  (4-20) 

TEQ  =  0   implies   TQV  -  T*y  =  Tg  (U.2d) 

Naturally,  all  the  above  definitions  apply  equally  well  to  a  III- I 
transition  by  replacing  III  by  I  and  I  by  III  everywhere.   The  symbol  T^ 
denotes  the  "ill- I  transition  time." 

Table  IX  contains  these  definitions  of  time  intervals,  which  are 
also  illustrated  in  Fig.  21. 

U„3  Calculation  of  a  Time  Interval  Over  a  Trajectory  by  an  Iterative  Formula 

It  is  not  possible  to  explicit  y  in  (3.59).   Therefore,  given  an 
initial  point  P^x^r  and  the  abcissa  xf  of  another  point  Pf,  in  order  to 
find  the  other  coordinate  y  of  Pf  such  that  PQ  and  Pf  are  on  the  same 
trajectory,  we  must  apply  an  iterative  procedure  to  (3-59);  the  time  intervals 
can  be  found  from  (3. 58). 

Equation  (3.59)  can  be  expressed  as 


■M    yi  a 


vr=fv^ 


Read:   Pn  "whose  coordinate  are"  (xQ,y0). 
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TABLE  IX.   DEFINITIONS  OF  TIME  INTERVALS  OVER  A  TRAJECTORY  (SEE  FIG.  21 ) 


Symbol 


Name 


T 


D 


B 


EO 


"OV 


TR 


I>* 


OV 


Delay  Time 


Definition 


t(-y) 


-  t. 


Comments 


Active  Time 


Complementary 
Time 


t(+£)  -t(-i) 


-  t(-  i) 
7 


If  circuit  was  settled  at 

XI1  at  V  ^xiq\*  V 
otherwise  t(x   )  f-   t  .   t 

is  the  instant  the  trigger 
is  turned  ON. 


Balance  Time 


t(xIII0,  u  =  1)  -  t(-  i) 


Settling  Time 


Trigger  Excess 
Overtime 


Trigger 
Overtime 


t  -  t, 
s    t 


tg  is  the  instant  when  P 
enters  N     for  the  last 
time. 


t(x- 


=  l)  is  the 


'IIIO'  M' 
instant  P  crosses  x=  x 
assuming  the  trigger  stays 
ON  (|i  =  1)  all  the  time  ' 


t(xIII(),  u  =  l) 


Transition 
Time 


Trigger 
Duration 


Optimum 

Trigger 

Duration 

and 

Overtime 

respectively 


*fl  "  ^   7) 


s    0 


tg  is  the  instant  the 
trigger  is  turned  OFF. 


t, 


0 


Defined  only  if 

*o  =  t(xio^ 

to  letting 


t.Q  -  t(x   ),  equivalent 


Notice  that  it  is  measured 

since  t   (not    t(xm)!) 

until  t  .  1U 

s 


Of  course,  only  approxi' 
mately  realizable. 


JQ 


t(xni0,  ix  =   1), 
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for  region  V  and  trigger  condition  |i.   Parameters  MQ,  M  ,  Qq,  Q^  are  gi 
Table  X. 

We  obtain  either  one  of  two  formulae 


lven  in 


y  n  =  M  -  Q  i 
^n+l    a    U 


a 


{h.ka) 


or 


a 


K 


£rv 


rM  -y  " 
a  Jn 


n+1    P 


V 


(l*.lrt>) 


The  only  difference  between  them  is  a  question  of  convergence. 

In  fact,  given  two  implicit  functions  f  and  g  of  y,  the  equation 

f  (y)  =  g(y)  (^5) 

can  be  solved  by  an  iterative  procedure  by  means  of  a  formula  such  as 

f(yn+1)  -  s(yn)  ^°6) 


i.e, 


n+1 


=  f_1[g(yj] 


(+.7) 


Let  y„  be  the  solution  of  equation  (U.5) 

The  iterative  formula  (U.6)  will  converge,  i.e.,  lim  yn  =  yf,  if  and 

n->oo 

only  if  there  is  a  number  e  >  0  such  that,  If  |yQ  -  yf  |  <  6,  there  is  a  positive 
number  K  such  that : 


No  loss  of  generality,  since  the  symbols  f  and  g  can  be  interchanged. 
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TABLE  X.   DEFINITION  OF  THE  PARAMETERS  OF  EQUATION  (4.3) 


Parameters  of  (4.3) 


M 


a 


%c 


K 


Q, 


P 


Expressed  as  Functions  of  the  Parameters  of  (3.59) 


(x  -  x*)\ 


a 


(xA  -  x*)\       -    v 

0     'a      J0 


(x  -  x*)k 


(3 


(x0  "  x*)^  -  yQ 


Comments:      a)     Region  V^    trigger  condition  u 

b)     x*  =  x 

Vu 

C)     V  =  W    V  W 

d)  (xo'yo^  =  any  Slven  Point  on  the  branch  of  trajectory 
under  conditions  Vu. 

e)  x  =  some  abcissa  such  that  there  is  an  ordinate  y 
satisfying  the  condition:   "(x,y)  is  on  the  same 
trajectory  vu  branch  as  (x  ,y  )." 

f )  y  is  to  be  found. 
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*V, 


n+1 


% 


<  K,   K  <  1 


(fc.a) 


n 


However,  it  is  clear  that 


^ 


lim 


n+1 


:'(yf) 


ST  ^ 


(4.9) 


where  the  prime  means  "derivative  with  respect  to  y_ . 


Now  (4.8)  and  (4.9)  imply  that  (4.6)  converges  if  and  only  if 


g'(yf) 


f1!^ 


<  1 


(4.10) 


Therefore,  one  of  formulae  (4.4a,  b)  converges  and  the  other  diverges.  There  is 
no  way  to  know  a  priori  which  one  will  converge,  since  we  would  need  the  solu- 
tion y  of  (4.3)  to  answer  this  question.  However,  assuming  we  start  from  a 
good  initial  guess  y  ,  if  we  calculate  yg  and  y^   from  both  formulae  (4.4a,  b), 
the  initial  tendency  should  be  clear. 

Another  way  would  be  to  differentiate  both  sides  of  (4.3)  with  respect 
to  y,  and  compare  the  two  results  for  the  initial  guess  y^  hoping  that  compar- 
ison at  y  would  yield  the  same  qualitative  result.   Call  f(y)  and  g(y), 
respectively,  the  sides  of  (4.3)  with  the  larger  and  smaller  absolute  value  at 
point  y  ,  and  take  that  of  formulae  (4.4a,  b)  which  conforms  to  (4.7). 

If  eventually,  the  selected  formula  diverges,  then  we  should  try  to 
improve  the  initial  guess  y1   and  repeat  the  procedure  outlined  in  the  previous 

paragraph. 

A  method  of  extrapolation  usually  allows  improvement  of  any  trial  yn 
using  the  previous  results  for  y^  and  y^2 :  from  (4.7)  we  write  (4.1l)  below: 
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Ayn  -  f'^gty  J]  -  y 


(4.11) 


So,  to  yQ_2  and  y^  there  correspond  respectively  the  variations 
Ayn_2  and  Ay^;  the  tvo  points  (y^, Ay^)  and  (y^Ay^)  define  a  line 


(y  ,Ay  ):' 
w  n'  ^  n 


■r-  /X-l         \  /  t%„ 


'n-2       y  VVl         n-2^n-2, 


^n  '  'nl  5~  "  V     +    (Vl  "  ^-2  5^  )  C^J 


Let 


*      -   0 


then: 


v     -  yn-2^n-l  -   Vl  ^n-2  „         % 

yn  Ay" T^ (4.13) 

Jn-1         jrn-2 

Try  substituting yn  into  (4.1l).   Stop  when  Ayn  is  small  enough. 

This  method,  even  though  more  involved,  would  speed  up  the  convergence, 

it  is  more  tolerant  with  respect  to  initial  guesses,  and  stabilizes  the  method 

to  the  point  of  usually  producing  a  convergent  sequence  of  numbers  y  -*  v  even 

n   ''s 

in  a  case  for  which,  if  directly  applied,  equation  (4,7)  would  diverge. 

If  the  problem  consists  in  finding  time  intervals  only,  and  we  are 
not  concerned  with  other  characteristics  of  the  trajectory,  then  the  trajectory 
equations  (3-57)  in  the  time  domain  could  be  used  directly.   They  can  be  written 
as 


^n  is  the  ex"trapolated,  or  expected,  value  of  Ay  . 
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Xf  .  A>T°f  ♦  Ap>T°f  +  x*  (U.lUa) 

y     .  A  X  eV°f  +  AftXRe^0f  (U.lto) 

Jf    a  a  p  P 

where  T   =  t  -  tn,  and  x*  is  the  abcissa  of  the  singular  point  corresponding 
to  the  region;  clearly  (xQ,y0)  is  the  initial  point.   So: 

x  n  =  Ae^f  A.^  +  x*  (4.15) 

n+1    a        P 

and,  by  the  use  of  a  method  like  the  one  expressed  by  equation  (4.13), 

xAT  ,  +  (x  T  .  -  x  ,T  ) 

m  "-1       n  n"1 n-1  n  (4.l6) 

n+1  Xn  "  Xn-1 

Now,  T  ,  would  be  used  in  (4.15).   The  numbers  x    and  x  are 
'   n+1  n_i-      n 

respectively  the  results  of  (4.15)  when  fed  with  T^  and  Tn.   Clearly 

AT  =  T  -  T  n  (4.17) 

n    n    n-1 

When  AT   is  small  enough,  the  process  is  stopped  and  T   corresponding 
n  wx 

to  x  could  be  fed  into  (4.l4b)  to  find  the  corresponding  value  of  y. 

There  remains  the  problem  of  how  to  find  a  fairly  good  initial  value 
for  the  iterative  procedure.   One  good  way  is  to  assume  that,  in  a  crude 
approximation,  in  regions  I  and  III,  P  moves  in  a  straight  line  towards  the 
singularity  (virtual  or  real),  and  that  in  region  II  it  moves  parallel  to  the 
asymptote  of  positive  slope. 

Therefore : 

In  regions  I  and  III,  with  x*  standing  for  the  corresponding 
singularity: 
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(x*  -  xf  )y0 

yi  =  "  x*  -  xQ  (^.18) 

In  region  II,  with  \  being  the  positive  natural  frequency: 

yi =  Vxf  -  V  +  y0  (^.19) 

Both  in  (^18)  and  (k.±9),    (x0,yQ)   is  the  initial  position,  and  we 
wish  to  find  the  first  approximation  yx  to  the  ordinate  yf  corresponding  to  the 
abcissa  x  . 

To  find  a  first  approximation  TQ1  to  the  time  interval  for  P  to  go 
from  (xQ,y0)  to  (xf,yf ),  consider  the  first  approximation  to  the  trajectory  as 
being  a  straight  line  from  UQ,yQ)   to  (x^).   We  know  that,  whatever  the 
trajectory  y(x)  may  be, 


T0f  =J  Vtt  **  (^20) 


xo 


Therefore,  if  y(x)  is  a  straight  line  with  slope  y',  going  through 
(x0,yQ)  and  (xf,y1),  then P  crosses  the  line  x  after  a  time  interval: 


T01  =  ?-^  (^.21) 


and  naturally, 


y'  =  -± 2  ,k   22) 

xf  -  x0  ^.^; 


NOTE:      In  a  general  form,    if  a  trajectory  is   a  straight  line  between  points 
fa:<WandMvV>   Tab  =  tb  -\±s 
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T  1'jfoh  (4.23) 

ab       y*  ya 


and 


v.    =  \U±  (4.24) 

J  x_    -  x 

b  a 


Whenever  x  and  x  are  in  different  regions,  we  have  to  proceed  by 
steps.   Assume  a  I  to  III  transition  with  xQ  in  region  I  and  xf  in  region  III; 
y  is,  of  course,  given,  and  y  is  to  be  found.  Then  (with  the  notation  shown 

in  Fig.  22): 

(i)  With  (x  ,y  )  as  initial  point  on  the  equation,  find 

y  .  at  x  =--l  ,    by  iteration. 
Ja}  a  7   _ 

(ii)  Use  equation  (3.6la.i)  to  find  yfe,  at  x^  =  -  -j  . 
(iii)  With  (x^y,  )  as  initial  values  on  the  equations,  find 


V   ,    at  x     =   +  — 
Jc'  c  7 


,    by  iteration. 

1 


(iv)  Use  equation  (3.6lb.i)  to  find  yd  at  xd  =  +  y|  . 
(v)  With  (x,,y,)  as  initial  values,  find  the  point  yf>  at 

x  =  x  ,  by  iteration,  assuming  that  P  crosses  x  =  xf 

while  the  trigger  is  ON. 
(vi)  Calculate  all  time  intervals  by  (3-58).   Then 

i0f    Oa    be    df 
(vii)   Consider  the  trigger  duration  T^,  and  suppose  that 

P  crosses  x  =  x  after  the  trigger  is  turned  OFF, 

either  for  the  first  or  second  time. 

(viii)  With  (x,,y,)  as  initial  values  of  the  coordinates, 

find  (x  ,y  ),  by  iteration,  where,  after  a  time  interval 
e  e 

T  =  T  -  (T   +  T   ),  trigger  turn-off  occurs. 
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(ix)  Find  ye  -  ye  -  yQ. 
(x)  With  (x  ,y  )  as  initial  values,  find  yf,  at  x  =  xf, 
by  iteration;  there  may  be  none,  one,  or  two  values, 
(xi)  The  time  intervals  can  be  calculated  by  (3-58),  and, 
in  this  case,  where  P  crosses  x  =  x  after  trigger 

t 


turn-off, 


T   =  T  +  T 
X0f    W    ef 


Of  course,  this  algorithm,  with  slight  modifications  of  detail,  can 
be  used  to  find  optimum  values  for  T„,  instead  of  having  Ty  as  part  of  the  data. 

This  same  algorithm  can  be  applied  to  a  III  to  I  transition,  after 
the  obvious  interchange  of  reference  to  regions  I  and  III. 

See  illustrative  examples  in  Chapter  6. 

k.k     Graphical  Constructions 
k.k.l     The  ^,X)  Plane  Method 

Equations  (3,15)  and  (3.U8)  define  the  transformation  of  variables 
from  (x,y)  to  ($,x),  and  can  be  written  equivalently  in  a  single  pair  of 
expressions : 


> 


<&  =  .   -1  .   [+v(x-x*)  -  y] 


p   a 


>  (^.25) 


NOTE:   Considering  that  y*  =  0, 


'  Of  course,  if  P  crosses  x  =  x  twice  after  trigger  turn-off,  we  must  know  up 
to  which  passage  of  P  by  x  =  x  we  wish  to  calculate  Tgf . 
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Here 


equation  (3-59)  and  (3.5l)  are  repeated  for  convenience,  with 


(t  -  t  )  replaced  by  T: 


VT 
$  =  $  e  a 


> 


X  =  Xne  p 

0     J 


(4.26) 


Xa       °o   XP    *o 


(4.27) 


We  see  that  it  is  possible  to  define  a  new  pair  of  variables,  say  $  and  %,    as 
follows : 


$  = 


X  = 


0 


"\ 


> 


(4.28) 


where  <t  Q   and  xQ  are  the  initial  values  of  €    and  X,  so  that  $   =  x  =  0  by 


0   *0 


definition. 

Therefore,  (4.27)  becomes 


a 


(4.29) 


This  means  that  in  the  (0,  x)  plane,  the  trajectories  are  straight 
lines,  and  intervals  can  be  marked  as  a  linear  scale,  either  along  a  trajectory 
or  along  a  vertical  or  horizontal  axis,  since  time  is  linear  with  either 
variable. 
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On  the  other  hand,  (4.25)  means  that  any  straight  line  in  the  (x,y) 
plane  is  also  a  straight  line  in  the  ($,x)  plane. 

From  (4.25),  it  is  clear  that  the  transformations  depend  on  the 
region  and  trigger  indices  V\i ,.   This  can  be  indicated  by  attaching  these  indices 
to$,  X,  and  also  to  $  and  x,  i.e.,  (4.25)  becomes: 


*\ 


$ 


1 


w  w[+^'(x"x^"yl 


>» 


K    pV|a    cuvu      ^ 


(4.30) 


Also  (4.27)  becomes: 


<S> 


£w 


VQ 


^avu   °vuO 


XVuO 


(4.31) 


And,  of  course,  (4.28)  becomes 


<& 


Vu 


VuO 


v 


"VU 


XvuO  , 


(4.32) 


And  finally,  from  (4.29) 


JV>u_  _  *JVu_ 


V^   W   W 


(4.33) 


Now 


,  for  each  value  of  Vp.  that  occurs  in  the  problem  (this  must  be 


known ) : 
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a)  Draw,  on  a  linear  graph  paper,  with  $   and  x   axes  marked  on  it,  the 

VfJ.         V|_l 

lines  corresponding  to  the  following  phase  plane  lines: 
(i)  The  coordinate  x  axis  with  a  scale  on  it. 
(ii)  The  direction  of  the  vertical  lines, 

b)  On  a  log-log  graph  paper,  mark  the  coordinated   and  x   axes,  and  the 
direction  of  the  trajectory  lines,  as  well  as  a  time  scale,  and  also  the 
scaled  curve  corresponding  to  the  x  axis, 

c )  Given  a  trigger  amplitude  W,  draw  the  corresponding  lines  x  =  +  -  in  the 

-  7 
linear  graphs  where  u  =  1,  and  from  these  graphs,  draw  the  corresponding 

curves  in  the  log-log  graph,  by  means  of  a  point -by-point  transportation. 

d)  The  log-log  graphs  representing  the  various  (ff    x"v  )  planes  provide  a 
means  for  very  fast  calculation  of  times  over  the  trajectory  corresponding 
to  the  given  trigger  amplitude. 

e)  The  set  of  log-log  papers  plus  the  linear  graphs  allow  a  fast  calculation 
of  trajectory  times  for  any  trigger  amplitude  (within  the  bounds  of  the 
graph  papers,  of  course). 

It  is  clear  that  this  method  is  advantageous  mostly  in  the  case  where 
several  calculations  must  be  performed  for  the  same  system. 


^.2  A  Simple  Method  on  the  (x.y)  Plane 

This  is  a  less  accurate,  but  faster  method,  and  more  versatile  in 
solving  problems  for  several  different  trigger  amplitudes.   It  will  be  called 
"the  phase  plane  method  A," 

It  consists  in  approximating  the  system  phase  plane  portrait  with 
four  straight  line  segments,  under  the  following  assumptions: 
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a)  We  assume  that,  in  regions  I  and  III,  P  moves  in  a  straight  line  towards 
the  corresponding  singularity  x   ,  V  =  I,  III,'  (1  =  0,  1;  and  that  in 

V(J. 

region  II  it  moves  parallel  to  and  in  the  same  direction  of  the  asymptote 
nearest  to  it. 

b)  The  discontinuities  of  y  at  x  =  -  -  and  at  x  =  '+.  -     can  be  calculated  by 
equations  (3.6la.i)  and  (3.6lb.i)  respectively, 

c)  As  a  very  fast  method  at  the  cost  of  obtaining  a  somewhat  poorer  accuracy, 
these  discontinuities  of  y  at  the  boundaries  of  region  II  (i.e.,  at 

x  =  +  — )  can  be  completely  ignored. 

d)  So,  as  illustrated  in  Fig.  23,  we  take  a  linear  graph  paper,  mark  on  it  the 
scaled  x  and  y  axes,  vertical  lines  x  =  +  -,  and  singular  points  XvQ,  and 
also  the  direction  of  the   asymptotes  related  to  the  saddle  point  of 

region  II. 

Then,  given  a  trigger  amplitude  W,  we  mark  the  points  Xvl  and  also 
P  .   Suppose  a  I  to  III  transition;  then  PQ  =  P:(xIQ, yQ). 

(i)  Draw  the  segment  of  the  line  PqX-^  contained  in  region  I. 


a 


(ii)  The  intersection  of  this  line  with  x  =  -  -  is  y 

(iii)  Find  y  by  (3.6la.i). 

(iv)  From  y  draw,  inside  region  II,  a  line  segment  with 

slope  \__T  (the  positive   asymptote  slope);  its 
pli 

intersection  with  x  =  +  —  is  y  . 
(v)  Find  yd  by  (3.6lb.i). 
(vi)   Consider  Pd  =  (+  r,  yd);  draw  the  line  P^j-^. 

Let  \,  be  the  slope  of  this  line,  and  call  it  X 


line 


i 


'  In  the  descriptions  of  these  graphical  methods,  references  to  a  line  shall, 
in  general,  be  made  using  the  symbol  for  its  slope. 
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W"< 


1 


**< 
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(vii)   Consider  P„:(x  ,y  )  be  the  intersection  of  lines 


•  f.v-f*jf 


PdXIIIl  aM  X  =  Xf ' 
(viii)   Calculate  Tjf  -  TQa  +  Tbc  +  Tdf  by  (U.2l),  supposing 

that  P  crosses  x  =  x  while  the  trigger  is  ON. 

X 

(ix)   Consider  the  trigger  duration  T  ,  and  suppose  that 
P  crosses  x  =  x  for  the  second  time'  after  the 
trigger  is  turned  OFF. 
(x)  Find 


W^Oa^1 


y  =  y-,e 

Je    d 


and  mark  P  : (x  ,y  ),  the  point  immediately  before 
e   e  e 

trigger  turn-off  occurs, 
(xi)  Find  yg  =  yg  -  yQ,  and  mark  pe-(xe^e)^  the  Point 
immediately  after  trigger  turn-off. 
(xii)  Draw  the  line  P~X~III(y  and  mark  Point  pf  :(xf>yf) 


of  intersection  of  the  lines  P JLjjjq   and  x  =  xf 
(xiii)  Find  T  _  by  (U.2l),  and  then  in  this  case  where  P 
crosses  x  =  x  after  trigger  turn-off, 


T0f  =  TW  +  Tef 


Of  course,  after  the  obvious  changes,  the  same  algorithm  can  be 
applied  to  a  III  to  I  transition. 


This  method  does  not  apply  with  acceptable  accuracy  if  P  crosses  x  -  xf  for 
the  first  time  after  trigger  turn-off. 
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4„^„3  An  Approximate  Method  on  the  (x,y)  Plane 

This  is  a  slightly  more  sophisticated  method  then  method  A;  it  will 
be  called  "the  phase  plane  method  B." 

We  will  use  the  phase  plane  equation  (3.1l)  which  we  repeat  here  for 
convenience : 


,    d  ,   -  c  x   b 

dy    Vu    Vu     Vu  a  .  , 

ta  =  '   avy  ^  "  if'        %  -   dv  +  ^nW  (3.11) 


To  plot  a  I- III  transition  in  the  phase  plane,  we  do  as  follows 
(i)  Find  y  f    and  x  ,  for  all  Vu  conditions, 
(ii)  Draw  the  lines  y  =  \   (x  -  x*)  for  x^,  *1TL>    xIII0-> 
and  x    ,  and  call  them,  respectively,  the  (3 __, 

pm>  piiio>  and  PIII1  lines° 

(iii)  By  (3.11),  we  find  ^  at  ?Q:{*1Q,yQ)t    and  call  \Q 
this  derivative, 
(iv)   Draw  a  line  with  slope  \     through  P  ,  and  call  this 
the  \     line. 
(v)   Check  if  the  \     line  intersects  the  3   line  inside 
region  I. 

If  this  is  so,  call  the  intersection  P  : (x  ,  y  ) 


and  call  P  : (-  — 
a  v   7 

PTn  and  x  =  -  — . 
II  7 


, y  )  the  intersection  of  the  lines 


Otherwise,  ignore  the  intersection  of  lines  Xn 


and  B__.  and  call  P  :(-  — 
II'  a  v   7 

lines  p  and  x  =  -  — . 


, y  )  the  intersection  of 


(vi)  Find  yb(x  =  -  i 


)  by  (3.6la.i) 


■13^- 


(vii)     Find  y   (x  =  +  — 


)   by  the   iterative  numerical  procedure 


c*      7 

using,  for  example,  one  of  equations  (k.k),    or, 


instead,  assume  the  trajectory  in  region  II  is  a 
straight  line  of  slope  X RTT°   (The  choice  depends 
entirely  on  a  compromise  between  accuracy  and  com- 
putation time.) 


(viii)  Find  yn(x  -  +  — 


)   by  (3.6ld.i). 


(ix)   By  (3.1l),  calculate  the  slope  — -  at  P -,:(  +  ~ 


+>yd} 


(for  region  III),  and  call  3  this  derivative. 

(x)   Draw  a  line  with  slope  X     through  point  P,  and  call 

it  the  A.,  line, 
d 

(xi)  Find  the  intersection  P  : (x  , y  )  of  lines  X     and 


III1' 


The  trajectory  in  region  III  with  a  very 


long  trigger  is  taken  as  segments  P^P^  and  P  X 

of  lines  A,,  and  3-r-r-r-i  •   Call  this  the  line  P,XTTT„. 
d      IIIl  d  III! 

(xii )  Do  as  directed  in  (vi )  to  (xiii)  of  method  A,  but 
modify  instruction  (x )  of  that  method  to: 
(x)  Find  P  :(x  ,y  ),  the  point  immediately  before 
trigger  turn-off  occurs  by: 


yel  =  yae 


WTOa+V] 


< 


ye2  =  yele 


VlIl[V(TOa+Tbc+V 


r — 


ye  =  yel>    if   yel  >   y3 


ye  -  ye2>    if   yel  <  y; 
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Again,  after  the  obvious  modifications,  this  applies  to  a  III  to  I 
transition. 

An  illustrative  example  is  presented  in  Chapter  6. 
Obs.:   Notice  that  this  Method  B  can  be  a  hybrid  numerical  and  graphical 
method.   Various  such  combinations  can  be  made  and  we  feel  that,  some  of  these 
combinations  may  be  good  compromises  between  speed  and  accuracy. 

4.5  Approximate  Analysis  of  Waveforms 

One  of  the  most  important  aspects  of  transition  waveforms  is  the  time 
duration  of  the  various  phases  of  a  transition  as  defined  in  section  4.2.   The 
exact  shapes  of  a  particular  variable  (voltage  or  current)  as  a  function  of 
time  is  less  important  than  its  general  characteristics,  such  as  delay  time, 
rise  time,  average  form  in  each  region,  minimum  and  maximum  values,  etc.   The 
exact  shape  is  important  insofar  as  it  influences  the  calculations  of  these 
characteristics,  especially  the  various  time  intervals  elapsed  between  the 
definite  changes  in  character  of  the  curve,  generally  described  by  changes  in 
the  values  of  the  pair  of  indices  Vu. 

Furthermore,  even  if  we  do  have  the  exact  (analytic)  solution  of 
(2.103),  it  will  not  do  us  much  good. 

We  can  solve  the  problem  for  the  waveforms  of  all  variables  based  on 
the  solution  of  (2.103).   But  then--besides  the  fact  that  (2.103),  and  therefore 
any  solution  based  on  it,  is  already  an  approximation  to  the  real  problem--the 
important  general  characteristics  of  the  waveforms  are  hidden  in  a  fairly 
cumbersome  analytical  formula,  which  would  take  a  considerable  time  of  tedious 
labor  to  plot . 

Our  aims  consist  mostly  in  analyzing  and  evaluating  an  existing  flip- 
flop  or  improving  its  design,  selecting  a  better  trigger  or  loading  circuits, 
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determining  optimum  trigger  duration  and  better  waveforms,  and  better 
understanding  the  operation  of  bistable  circuits. 

For  these  purposes,  an  approximate  plot  of  the  several  variables 
which  could  be  obtained  in  a  reasonably  short  time  would  be  far  more  useful. 

In  this  section  we  will  suggest  some  methods  by  which  such  graphs 
can  be  obtained. 

U.5.1  Collector  and  Base  Currents 

Assume  that  an  approximation  to  y(x)  (phase  plane)  has  been  obtained,, 
consisting  of  four  line  segments,,  one  for  each  value  of  Vu  (three  line  seg- 
ments in  the  case  of  optimum  trigger  duration). 

This  can  be  obtained  either  by  the  second  graphical  method  described 
in  section  k.k.2,    or,  if  time  durations  are  extremely  important,  by  calculating 
time  durations  with  one  of  the  iterative  techniques  described  in  section  k.3, 
and  then  using  (4.21 )  and  (4.22)  to  determine  the  position  of  line  segments 
which  would  result  in  the  same  trajectory  time  intervals. 

a)  With  (4.21),  several  points  can  be  marked  over  this  approximate  trajectory 
constituting  indeed  a  (nonlinear)  time  scale. 

b)  Or  else,  considering  that  x(t)  has  the  form 


x  =  A  +  Be 


(4.3*) 


A,  B  and  \   can  be  found  for  each  value  of  Vu. 

By  a)  or  b)  above,  or  any  equivalent  method,  plot  x(t)  and  y(t). 

(i)  Collector'  Currents 

If  (2.93)  is  assumed,  an  approximate  graph  of  the  collector  current 
variables  w,  k  -  1,  2,  is  immediate  for  they  will  be  constants  outside  region  II 
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(either  0  or  1,  whatever  the  case  may  he)  and  will  be  linear  functions  of  x 
inside  region  II,  so  that,  by  just  assigning  new  scales,  the  two  curves  w  (t) 
can  be  obtained. 

More  accurate  curves  can  be  obtained  by  using  equations  (2,56)  and 
(2.57),  which  will  yield  results  in  closer  approximation  to  the  real  transistor 
currents,  than  the  model  represented  by  equation  (2.IO3).   Use  of  a  graph  of 
the  tanh  x  would  allow  a  completely  graphical  procedure . 

(ii)   Base  Currents 

Here  use  of  equation  (2.9k)   yields  graphs  of  z  ,  k  =  1,  2,  with 
almost  equal  ease.   The  first  term  is  directly  proportional  to  the  corresponding 
collector  current  w  ,  and  the  second  term  is  proportional  to  y(t)  in  every 
region,  since  cp'(x)  is  a  constant  in  every  region. 

In  this  case,  use  of  (2.58)  and  (2.59)  to  improve  accuracy  would 
hardly  be  justified. 

J+.5-2   Collector  Voltages 

By  inspection  of  Fig„  2k   we  get  immediately: 

k  =  1,  2 
dvk 

voi  =  \  +  Rikcik  IF  "  Rik^k  +  Sk  +  V^    y-i>e         fr-36) 

I  t  k 

where  t'  is  the  nonnormalized  time  variable.   Normalizing  as  before  and  setting 

ui  =  \  »»▼  v   J  =  I*  2  (4.37) 


J   2   oj- 


we  get 
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FIGURE  24: 
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PASSIVE     NETWORK 
YIELDS   THE    EQUATION 
FOR  THE   OUTPUT  VOLTAGE. 
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VVTV  2Pk  F^  K  +  \)  -  ^k  f#     (M8) 


ok  ok 


,TE 


with  k  =  1,  2}   I   =  1,  2}   $>  f   k.   As  before,  we  can  interpret  0   in  terms  of 

an  equivalent  current  source  of  strength  s  and  a  parallel  conductance  G  , 

k  k 

according  to  (2.79),  repeated  here  for  convenience .   Assuming  a  rectangular 
trigger,  and  making  use  of  the  index  |i  =  0,  1; 

•*- **-¥*"*    k  =  1'2  (4-39) 

k 
Substituting  (4,39)  into  (4,38),  and  using  y  =  x  ,   we  get: 

u*  =  ^  +  RiAA  +  ^  yk  "  2?k  FT  <sk  +  zk>  "  sPk  f#    ^°) 

ok  ok  E 

with  k  =  1,  2;  i  =  1,  2;  i  /  k. 

This  equation  is  very  general,  and  allows  one  to  find  both  collector 
voltages  of  a  general  Eccles -Jordan  flipflop  (symmetrical  or  nonsymmetrical) 
if  x1   and  xg  are  known  and  also  holds  for  the  asymmetrical  flipflop,  by  dropping 
the  indices  k  and  &. 

For  the  moment  we  shall  focus  our  attention  on  the  asymmetric  flip- 
flop  and  on  the  symmetric  Eccles -Jordan. 

In  the  first  case  (asymmetric  flipflop)  we  get: 

T  R  R  T 

u  =    (1  +  R.G  )x  +  ^  y  -  2p  ^  (s  +  z)  -  2p  -*JL        (kM) 

o  o  E 

And  in  the  latter  case  (symmetric  flipflop)  we  get  by  subtraction, 
and  setting 


u  -  «g  -  ^  (^2) 


u  -  (1  +  R.G^x  +  ^  y  -  2P  J  z  -  2P  J  ,W  +fe^)    (^3) 


with 


M 


0,  if  trigger  is  OFF 

1,  if  trigger  is  ON 


These  equations  immediately  suggest  the  procedure  for  obtaining  the 
collector  voltage  variable  u(t)  from  x(t),  y(t),  z(t)  and  Wj  it  is  clearly  a 
very  easy  graph  to  obtain  from  the  preceding  ones,  since,  in  each  region,  it 
consists  of  a  constant  plus  a  linear  combination  of  the  previous  curves,  with 
only  the  constant  and  possibly  the  coefficient  of  x(t)  having  different  values 
for  the  two  distinct  trigger  states. 

k.6     The  Influence  of  Parameters  on  Transition  Times  -  -Simplified  Equations 

We  would  like  to  have  some  qualitative  notion  about  the  effects  of 
the  various  parameters  upon  the  overall  transition  time.   We  are  also  interested 
in  learning  something  about  the  total  charge  fed  into  and  removed  from 
transistor  bases  and  capacitors,  and  their  relation,  if  any,  with  transition 
times.   Besides  that,  some  characteristics  of  waveforms,  such  as  maximum, 
minimum  and  settled  levels  of  collector  voltages  and  peak  base  currents  also 

interest  us. 

At  this  point  we  must  stress  that  we  are  searching  for  more  qualitative 
criteria,  i.e.,  first  order  approximation  formulae  which  could  help  considerably 


t  Notice  that  x  =  x±  -   xg,  i.e.,  the  order  of  the  indices  is  reversed  in  the 
two  definitions. 
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in  the  evaluation  and  understanding  of  flipflops,  and  not  for  exact  (or  good) 
engineering  design  formulae .   In  this  respect  the  character  of  this  section 
is  entirely  different  from  the  general  character  of  this  dissertation. 

4.6.1  The  Optimum  Flipflop 

Let  us  assume  (since  it  is  possible  in  principle)  that  a  flipflop  has 

been  constructed  such  that,  if  the  trigger  duration  is  optimum.,  i.e.,  if 

Xe  =  XIIIO>  then  ye  =  °>  where  xe  and  ye  are  the  coordinates  of  P  ,  the  position 
of  P  immediately  after  trigger  turn  off. 

For  this  flipflop,  we  can  say  that,  in  a  first  order  approximation, 
the  transition  time  is  given  by: 


yn       xTn  -  xx^         y^         X„,.,  -  x. 


'IIIl 


1 


0 


■II    10 


0 


(4.44) 


■mi   ino 


Let  us  assume  that  xim  -  xIII0  =  x^  -  *       =  Ac.   Then 


TR   y0     (Ac)2 


&M) 


From  Table  II: 


Ac  = 


1  +  R5, 
s  1 


/    1  +  R  G,\ 


(4.46) 


■II    IIIl   Vl  +  R  G, 

s  1. 


)2[(n'¥  +  J)2  -  (Htf] 


(^7) 


This  comes  from  assuming  a  straight  line  approximation  to  the  trajectory, 
neglecting  the  active  time,  and  using  equations  (4.23)  and  (4.24)  in  regions 
I  and  III.   Notice  that  state  symmetry  is  not  required  for  this  approximate 
equation  to  be  valid. 


So, 


-1*4-2- 


where : 


nfW  +     1 


1  +  RsGl 

1  +  RsV 


l(J  +  Hty) 


"|2 


(n'W  +  of   -    (Ht)2 


(k.kQ) 


,)     n"    =  2 


R 

c 


b)     cpk  =    [1  +  2Bk  +  p] 


o 


d)     J  = 


% 


if  it   is   an  asymmetric  flipflop 


Ap  =  cp     -   cp   ,,      if  it   is   a  symmetric  flipflop 


)  H  = 


1,    if  it  is  an  asymmetric  flipflop 
2.3        if  it  is  a  symmetric  flipflop 


Assume:  a  =  axxi   =  a^  and  mI  =  mIII  ~  m° 


m_ 

y    =  W  =  -  W 

J0       a      a 


(^9) 


where  a  and  m  are  given  in  Table  I; 


a)  m  =  2p 


T. 

10 


b) 


a  = 


T.T 
1  O 


And  therefore 
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y     =  2P1  w  =  2pr_ 

^0        T_.    W        R  C. 

Oi  o   1 


(^.50) 


So 


y0    Td  ♦  rsGi)  Lr0    I    - 1  +  RsgJ    2W  . 


(U.51) 


It   is   clear  that 


f  ~  1,        if       (1  -  a)R    «  aR 

s  o 


so  that  Hty  *»  Ho 

And  assumption  of  state   symmetry  eliminates  biasing  from  the  formula: 

J  =  0 


From   (4„48)   and   (4.51)  we  get: 


R  C. 
o   1 


"TR       T(l+R  g   ) 
si 


rR        / 
u  o       \ 


1  +   R  Gn\  _      „. 
s   1  \  J  +  Hlr 


1   +   R  Gn 
s   0. 


2W 


W 


#n/ 


R  /  1   +  R  GJN  _        " 

__s  A        s_l  \G  +  H 

R  V  '   1  +  R  Gn  /    2W 

■  o  \  s  oy 


v4W  +  J)2 


v2 


(^52) 


And,,    if  t  "  1  and  J  =   0,   we  get 
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R  C. 
o  1 


TR 


T(l  +  R  G,  ) 
v     s  1 


^s   HK 

R  +  2W 
•  o 


•  fou 


\      hk 

R  +  2W 
o 


w 


k2W 


(*.53) 


where 


K  =  1 


1  +  R  Gn 
s  1 

1  +  RG, 
s  0 


As  a  final  simplif ication,  if 


Go  =  Gi 


(equivalently,    K  =  0) 


we  get: 


T, 


R  C. 
s    1 


1  T(l  +  R  G.) 
s    1' 


bn> 


1  - 


IR    \2 
o   \ 


2R  ¥, 
•v     s  / 


{h.5k) 


Let  us  suppose  that  we  have  a  fairly  large  trigger,  and  that  p  is 

i 


also  sufficiently  large1  so  that 


«  1 


(h.55) 


Then  taking  the  first  two  terms  of  the  series  expansion  of  the  argu- 
ment of  the  natural  logarithm,  and.  then  taking  the  first  term  of  the  series 
expansion  of  the  logarithm  itself,  we  obtain: 


f  These  conditions  are  not  a  property  of  f lipf lops .   In  fact  HRQ/2R  W  close  to 
1  is  practical,  since  W  «  10"1  is  practical.   However,  we  assume  that  W  has 
been  chosen  large  to  speed  up  the  transition. 
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T~~  -      -*-2i-  (^56) 


s         S  X 

Equation  (4.56)  should  be  an  acceptable  first  order  approximation 
for  the  transition  time  whenever  the  flipflop  and  trigger  satisfy  all  of  the 
assumptions  leading  to  it.   Observe  that  there  is  a  hierarchy  of  equations  with 
more  and  more  restrictive  assumptions;  all  of  them  assume  the  optimum  flipflop 
described  before.   Then: 

Equation  (4.52)  is  very  general;  the  only  assumption  is  that 

the  flipflop  is  either  symmetric  or  asymmetric; 

Equation  (4.53)  assumes,  further,  that  ty  =  1,    and  J  =  0; 

Equation  (4.54)  assumes,  still  further,  that  G  =  G  , 
i.e.,  that  K  =  0; 

Equation  (4.56),  besides  the  above,  assumes  (4.55)  to  be  valid. 

Also  remember  that  all  times  are  normalized  with  respect  to  T,  so  that 
nonnormalized  times  would  not  include  T  in  the  denominator.  For  example,  (4.56) 
would  read : 


H2 


T2 


Mi  +  r  o.  )vr     l  i       oi-' 


Remembering  that  W  =  — — ,  we  get  : 


»in      4R  C.(l  +  R  Gjr 
s  iv     s  1' 


tTr  = ; ; 72  (^-58) 


Another  expression  for  TTR  can  be  obtained  as  follows:   Consider 
that  the  charge  variation  A},  of  C.  between  the  two  stable  states  is  given  by 
(see  O.69)): 

*k  "  ^IIIO  -  V  '  Ci  =      I^         I  ,(1  +  H  o  )   (1  +  W 

2  '         2       s  0 

(^.59) 


Therefore, 


/»AA  2    (1  +  RSGQ)2 


R  C. 
s  1 


(U.60) 


If  G  =  G  =  0  (i.e.,  the  trigger  circuits  are  perfect  current 
sources),  and  if  ty  ^  1, 


R  C. 
s  1 


(4.61) 


This  equation  should  give  us  a  somewhat  crude  but  satisfactory  first 
approximation  to  the  transition  time. 

We  insist  that  use  of  (4.59)  should  always  be  cautious,  since  some  of 
the  assumptions  made  in  its  derivation  are  somewhat  vague,  and  others,  if 
ligitimate,  will  seldom  be  fulfilled.   So,  (4.59)  is  usable  for  estimating 
results,  i.e.,  as  a  kind  of  figure  of  meritj  it  is  definitely  not  a  design 
formula. 


t  This  approximation  is  made  under  our  assumptions?  ?  -  0,.  ♦  =  1,  GQ  ==  G±,    s< 
there  is  a  cancelling  out  in  the  second  term  of  (h.bO). 


At  first  glance  it  seems  strange  that  the  transition  time  does  not 
depend  on  T  in  a  first  approximation.   In  other  words ,  it  seems  strange  that  T 
is  not  a  factor  of  prime  importance  in  the  transition  time.   The  following 
discussion  should  account  for  this  observation „ 

First  of  all,  since  the  beginning,  we  have  completely  ignored  the 
active  region,  and  second,  we  have  assumed  that  y,,  the  ordinate  of  P  when 
entering  region  III,  would  be  such  that  y  =0,   Of  course,  if  T  has  some 
influence  in  the  active  region,  it  will  determine  the  value  of  y  }    in  assuming 
yd  to  have  a  convenient  value,  we  have  ignored  the  effects  of  the  transistors, 
or  in  a  better  way,  we  have  assumed  that  there  is  a  relationship  between  the 
transistors  and  the  passive  network  such  that  the  optimum  flipflop  assumption 

is  verified.   In  this  sense,  T  should  be  related  to  T  .,  and  therefore 

01' 

(especially  if  this  relation  were  found  to  be  linear)  T  .  could  be  replaced  by 
its  expression  in  terms  of  T „   Then  T  would  be  the  prime  factor  in  all  those 
equations,  and  T   would  not  appear  at  all.   We  could  also  have  a  linear  com- 
bination of  both  parameters.   That  we  have  started  using  T  .  was  a  question  of 

01       * 

convenience;  the  assumptions  made  establish  a  certain  relationship  between  T 

oi 

and  T .   We  conclude  that,  after  all,  T  is  a  very  important  factor  in  the 
transition  time. 

Even  more  important  than  the  approximation  of  Tmo  furnished  by  these 

TR 

formulae  in  the  case  of  an  optimum  flipflop,  is  the  following  consideration: 
(i)  Even  if  the  flipflop  is  not  optimum,  the  transition 
time  should  not  be  substantially  different  from  the 
results  obtained  by  the  use  of  these  formulae.   They 
would  be,  at  any  rate,  a  first  order  approximation  to 
the  transition  time. 
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(ii)  They  would  certainly  "be  true  in  a  qualitative  sense, 
i.e.,  as  indications  of  the  relative  effects  of  the 
various  parameters,  as  well  as  the  order  of  magnitudes 
and  directions  of  change. 


k.6.2     The  Total  Charge  Interchanged  Between  the  Transistor  Bases 
These  were  established  in  Chapter  2,  equations  (2.20)  as 


Therefore,  the  total  charge  variation  is 

^B  =  m(Xh 
Some  relations  can  be  established  here,  such  as 


mh  -  ~f  =  —   &  +   KsG0 
01 


But  let  G  =  GQ;  then,  from  (k.6o) 


*U  "  ^  '  *  •   C1  +  R,Gn) 


Also, 


TR 


(^•63) 


1  (1  +  R  Gn)  (h.6k) 


TTR  "  1^2 It/         R^T        (1  +  W) 


i  2I.T  J  (1  +  RsWi 


Since  the  charge  that  enters  one  base  is  equal  to  the  charge  that  leaves  the 
other,  we  can  talk  about  the  charge  transferred  between  the  bases,  though 
this  transference  is  really  only  a  mathematical  cancellation,  not  physical 
transference. 
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If   G1   =   0, 


T™=1 


HA}  _T    A  2 

B  oi N 

21  T 

Li 


R  C. 
s    1 


(^.67) 


Also,    it   is   clear  that,    if  G  =  0, 


1  Ol  ^1 


(4.68) 


which  is  very  illustrative  of  the  type  of  condition  relating  the  transistor 
parameter,  the  passive  network  parameters,  and  the  two  stable  states. 
Notice  that  for  the  symmetric  flipflop, 


qB  =  qBl  "  qB2 


qi  =  qil  -  qi2 


%  =  **B1  "  ^B2 


.''.  Aq.  =  Aq    -  £q 

l    ^ll    ^i2 


^ 


r  (^.69) 


j 


W  -  W1  -  W£  but  the  triggers  W^  and  W2  are  assumed  to  occur  simultaneously 


^.6-3   Collector  Voltages --Maximum,  Minimum  and  Settled  Values 


(i)  Maximum:   v^^  =  V^  +  R^  +  q^) 


(ii)  Minimum:  v  .  .   =  v   +  R  (l.fl  +  (l  -  a)L 

kmin    Ck    ov  tl/,  J& 


(iii)  Settled:  v„,w  -•■   V„n  +  R  aL 


> 


C1I0    CI    o 


(iv) 

(v) 

(vi) 


1III0  =  V01  +  Ro(l  -  aK 


VC2I0  =  VC2  +  V1  -  ^h 


VC2III0  =  VC2  +  RoaIE 


y  (^.70) 


j 
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where:    k  =  1,  2;  I   =  1,  2;  I  t   k 

v     =  collector  voltage  of  transistor  T,  when  x  =  x  n 
CkVO  K  vu 

V   =  collector  supply  voltage  of  transistor  Tk 

In  case  of  the  asymmetric  flipflop,  drop  the  indices  k  and  k,    and 

v     =  V   . 
ClvO    CI 

k.d.k     Peak  Values  of  Base  Current 

Considering  the  optimum  flipflop,  and  the  approximate  model  whose 
equation  is  (2.103),  it  is  clear  that  the  peak  base  current  would  be  given  by: 


(Ml) 


z*          = 

peak 

7yd 

d  ==  y0 

XIII1 
Ax 

z*  f  if  the  flipflop  is  asymmetric 

with  z*  =  i 

*  -   z*,   if  the  flipflop  is  symmetric 


z 


1    2 


and  the  symbol  "*"  means  the  component  of  the  current 
corresponding  to  base  charge  variation 

y 

The  approximate  form  of  ~,   from,  say,  (^.5*0,  and  the  approximate 


form  of  x    ,  yield 


t(1  +  R  G,  )       nW  ,.         . 

s  1       nw  (1+  72) 

yd  :      R_C,   '  '  (1  +  \\) 


s  1 


so  that 


Notice  that  n  =  pn' . 
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T   «*,  .  „   T 


Zpeak  -  2  —   7W     °r     Speak  =  2  ~  7\  ^) 

01  01 

We  stress  that  this  value  of  z     refers  to,  not  only  the  approximate 
model  of  equation  (2. 103),  but  to  this  model  with  all  the  restrictions 
implicitly  imposed  in  the  evaluation  of  y . . 

Therefore,  such  an  expression  is  specially  meant  to  give  us  an 
acceptably  close  idea  of  the  values  of  the  base  current  in  any  given  case,  when 
just  a  fast  estimate  is  required „ 

^-1     The  Problem  of  Circuit  Optimization 

Whenever  one  tries  to  state  a  problem  of  optimization,  besides  a  clear 
statement  of  what  is  to  be  optimized,  two  basic  questions  must  be  answered. 
First:   "Under  what  criterion?" 
Second:  "What  are  the  constraints?" 

The  amount  of  material  written  on  these  optimization  questions  is 
very  large .  We  shall  not  try  to  find  complete  answers  here,  but  rather,  to 
open  the  discussion  by  some  pertinent  approximations. 

The  first  question  is  what  characteristics  we  could  wish  to  optimize: 
trigger  duration  and  amplitude  (if  not  its  waveform!),  circuit  parameters,  or 
the  transistor  characteristics. 

This  first  question  being  decided,  we  could  go  on  to  the  second 
question,  and  try  to  be  specific  about  stating  an  optimization  criterion,  i.e., 
an  interpretation  of  the  word  "improvement!" 

Of  a  host  of  possibilities,  we  can  state  the  following  three  as 
examples : 
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a)  The  time  interval  between  the  moment  when  the  trigger  is  turned  ON  and  the 
moment  the  collector  voltage  is  settled  (under  what  criterion  to  decide 
this?)  is  to  be  minimized.   Call  this  time  the  "collector  voltage  switching 
time.,  "  T   . 

b)  The  time  interval  between  the  moment  when  the  trigger  is  turned  ON  and  the 
moment  the  base  (or  base-to-base)  voltage  is  settled,  that  is,  the 

"base  voltage  switching  time/'  T  ,  is  to  be  minimized. 

c)  Instead  of  minimizing  switching  times,  one  might  wish  to  have  a  given  delay 
time  and  a  minimum  active  time,  or  a  minimum  switching  time  with  a  given 
delay. 

And  so  on  I   The  above  illustrates  the  point. 

We  have  already  attempted  to  approach  question  number  one,  in  a  very 
tentative  way,  with  respect  to  the  variable  x  (see  3°5d  and  e)  in  defining,  for 
a  special  purpose,  a  concept  of  "optimum  trigger  duration, "  which  was  related 
to  the  minimization  of  a  defined  "transition  time"  TTR,  for  the  variable  x. 
The  difficulties  were  apparent  and  that  discussion  stands  as  a  good  example  of 
the  issues  involved. 

The  second  question  is  usually  easier  to  settle,  since  constraints 
are  naturally  stated  either  as  inequalities  or  as  relations  between  the 
variables,  or  some  other  mathematical  statement.   To  incorporate  constraints  in 
an  optimization  algorithm  is  still  another  thing;  but  it  has  been  done  success- 
fully for  several  problems,  and,  once  stated,  there  is  no  a  priori  reason  to 
expect  the  problem  to  be  intractable.   The  theory  presented  so  far  suggests  a 
number  of  techniques  to  approach  optimization  problems,  once  they  are  stated  in 
a  mathematical  form. 

As  a  last  observation,  it  is  worth  reminding  ourselves  that  problems 
of  optimization  tend  to  raise  questions  of  existence  of  solutions  (realizability) 
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and  that  nothing  has  been  said,  for  example,  about  the  realizability  of  our 
hypothetical  "optimum  flipflop"  so  liberally  used  (as  an  approximation  device) 
throughout  this  Chapter  4. 

k.Q     Summary 

After  defining  a  nomenclature  for  time  intervals  over  a  phase  plane 
trajectory,  we  have  presented  some  methods  for  the  calculation  of  points  and 
time  intervals  for  a  given  trajectory., 

An  iterative  numeric  procedure  allows  the  exact  calculation  of 
y  (x  )  and  t,  (x  )  if  y  (x  )  and  t,  (x,  )  are  known,  for  any  given  pair  of  abcissae 

x  and  x.,  . 

a      b 

Similarly,  a  fairly  sophisticated  graphical  construction  using  two 
variable  transformations,  (x,y)  -*  ($,x)  "*  ($jX)>  was  also  presented,  and  shown 
to  yield  accurate  results,  given  the  limitations  of  a  graphical  construction „ 

A  more  naive  construction  on  the  phase  plane  was  described,  which 
yields  somewhat  less  accurate  results,,  but  is  extremely  simple  to  apply. 

It  was  also  suggested  that  some  hybrid  constructions  graphical  and 
numeric,,  might  be  ideal  for  accuracy  and  practicability  of  use„ 

A  graphical  procedure  to  obtain  fairly  good  plots  of  collector  and 
base  voltage  and  current  waveforms  was  described „ 

Engineering  interest  in  simple-minded  formulae  which  can  work  as  rules 
of  thumb  for  the  rapid  evaluation  of  circuit  characteristics  has  led  us  to 
discuss,,  by  means  of  an  ultra- simplified  model,  a  set  of  such  relationships  „ 

Finally,,  the  optimization  problem  was  proposed  in  a  first  approach 
discussion,, 


5.   EXTENSION  OF  THE  THEORY 

5.1   Introduction 

We  have,  so  far,  confined  ourselves  to  the  asymmetric  and  the 
symmetric  flipflops  subjected  to  a  rectangular  trigger,  and  also  we  have 
implicitly  assumed  that  neither  C±,    Cq   or  T  is  zero. 

In  this  chapter  we  shall  discuss  the  problems  involved  in  applying 
this  theory  to  other  situations,  and  indicate  the  methods  and  modifications 
involved. 

5.2   Case  When  T  Is  Negligible 

This  is  a  very  unlikely  possibility,  but  it  may  happen.   In  case  it 
does,  we  can  take  T  =  0  as  a  good  approximation.   Then,  the  coefficients  of  the 
equilibrium  differential  equations  apparently  are  meaningless! 

However,  looking  back  to  how  these  equations  were  established,  we 
will  see  that  T  was  used  only  as  a  convenient  time  normalization  constant.   Of 
course,  if  it  is  too  small  (or  too  large,  as  we  shall  seel)  it  ceases  to  be 
convenient,  and  some  other  time  interval  T  (such  as  TQ±,    for  example)  could  be 
used  as  a  time  normalization  constant. 

In  performing  this  renormalization  of  time,  we  replace  T  with  T  in 
equation  (2.3*0  and  on  all  related  equations  from  then  on.   By  letting  T  =  0  in 
equations  (2.29)  and  (2.30),  lQ±   and  l^   disappear  from  the  expressions  for  i^ 

and  iB2. 

The  result  of  this  is  that  the  charge  storage  in  the  base  along  with 
its  related  current  will  be  negligible,  and  only  the  recombination  component 
of  the  base  current  needs  to  be  considered.   Then,  the  equilibrium  equations 
will  not  contain  terms  like  zQ  and  §Q.   Except  for  this,  the  theory  is  exactly 
the  same,  and  applies  exactly  in  the  same  way. 
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5.3  Case  of  Negligible  External  Capacitances 

Again  we  have  a  possible,  although  unlikely  situation,  which  becomes 
important  especially  because  it  can  be  solved  in  a  special  way,  i.e.,  not  just 
an  extension  of  the  general  theory. 

By  making  t±t   1 ±q}   t q±   and  Tq  all  zero  in  equations  (2.87),  (2.88) 
and  (2.90),  we  would  get,  respectively: 

For  the  asymmetric  flipflop,  from  (2.87) 

y  =  r2  j  (!  -  PHanhx  -  -  (1  +  RsG(i)x  +  (l  +  2B  +  p)  +  2  -2.  s  I   cosh2x 
s  o  -^ 

(5.1) 

For  the  symmetric  Eccles -Jordan  flipflop,  we  get  from  (2.88) 


.2 

cosh  x 


of  1  R   1 

y  =  r"  {(!  "  PHanhx  -  —  (1  +  R  G  )x  +  (^  -  BQ)   +  ^  sf 
s  0  J 

(5.2) 

For  the  nonsymmetric  Eccles-Jordan  flipflop,  i.e.,  the  most  general 
case,  directly  from  (2.90): 

(1  +  RskGku)xk  =  Pk  {C1  +  2Bk  +  Pk)  -  C"1^^  "  ^k)tanh  X 

(5-3) 

(    i  >k  Rsk      ,2    0  Rsk 
+  (-1;  - —  y  sech  x  +  2  - —  3. 

ok  ok 

so  that,  since  x  =  x  -  x  ,  we  get: 
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P  f  TOT 

=  (i  +  r1  G  )  {(1  +  2Bi +  pi}  +  (1 '  pi^tanhx  •'  -r.  ysech  x  + :!  f: 


bl 


bl 


-Trrra{(1  +  aBa  +  p2)-(1- 


s2  2u 


Rs2      2      Rs2 
p  )tanhx  +  ^ —  ysech  x  +  2  - —  s^ 

o2  o2 


(5A) 


w 


ith  the  result  that: 


.2 

cosh  x 


y 


■*!    Si 


P2Rs2 


"t1  +  RslGm^oT  +    (1  +  Rs2G2^ 


p1(l   -   Px)       p2(l   -   P2)" 

1  +  R  _Gn      +  1  +  R  0G0 

si  lu  s2  2|-iJ 


tanhx  -  x 


Pl(l  +  2BX  +  Px)        p2(l  +  2B2   +  P2) 


1  +  R  ,Gn 
si  lu 


1  +  R  0G0 
s2  2u 


+  2 


PlRslSl 


P2Rs2S2 


L(l+RslG^)R01 


^1+Rs2S^R02- 
(5»5) 


So,    in  every  case  we  have  y  =  f    (x),°    of  course  we  are  assuming  that 
s  (t)   is   a  rectangular  function.      Therefore,    we   can  find  x(t),    or  better  t(x), 
by  the  formula: 


t   -  t, 


Ve5 


St,       where       v  (|)  =  f  (!) 


(5.6) 


and  we  mean  that,  if  u  changes,  at  a  certain  point,  we  must  find  its  abcissa 

x  ,  and  continue  the  integration  after  x  with  the  new  function  of  x. 
a  a 

It  is  easy  to  see  from  equations  (2 087)  and  (2.88)  that  these  cases 

are  still  exactly  solvable  even  if  only  CQk  =  0  and  Cik  /  0,  in  the  same  way  as 

when  C   =  C   =  0.   The  only  difference  is  that,  in  (5.1)  and  (5-2),  instead 
ik    ok 

R     P 
of  _2  cosh  x  as  a  factor  on  the  right-hand  side,  we  shall  have  the  following 
R 
s 

modifications  s 
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For  the  asymmetric  case,  from  (2.87),  replace, 

R 
in  (5.1),  ~  cosh  x   by   L _  (5o?) 

s  sech  x    i 

R   — « +  — 

s    R      T 
o 

For  the  symmetric  case,  from  (2.88),  replace, 

Ro     2              1 
in  (5.2),  —  cosh  x   by   ± —  (5.8) 

s  sech  x    i 

Rs  —R  +  2T 

o 

Note  in  (2.90)  that,  if  CQk  =  0,  even  the  nonsymmetric  case  is  con- 
siderably simplified,  since  it  will  be  reduced  to  a  second  order  case,  i.e., 
two  first  order  equations.   Then,  if  R   =  R    the  system  can  be  exactly  solved, 

S-L      Sd  * 

just  like  the  symmetric  case.   Otherwise  it  would  be  approximately  solvable,  like 
the  nondegenerate  symmetric  system. 

5 ok  Nonsymmetric  Eccles-Jordan  Flipflops 

The  difficulty  in  the  case  of  the  nonsymmetric  Eccles-Jordan  flipflop 
is  that  there  is  no  way  (except  for  some  extremely  fortunate  coincidence)  to 
reduce  the  two  equations  (2.7*0  in  x±   and  x£  into  a  single  equation.   The  fact 
is  that  this  circuit  has  one  more  degree  of  freedom  and  there  is  no  possible 
reduction  to  the  previous  cases.   Nevertheless,  we  can  do  something  about  solving 
the  system.   Suppose  that  we  carry  out  an  approximation  of  equations  (2.7*4-), 
taking  cp(x)  instead  of  tanh  x  just  as  we  have  done  to  obtain  equation  (2.IO3). 
The  result  will  be  the  pair  of  equations  expressed  by  (2.104),  whose  coefficients 
are  shown  in  Table  1.2,  and  which  is  repeated  below  for  convenience: 
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DO  O  i    OO         .,    ,     O  , 

akvxk  +  bkx  +  ckx  =  akvx  +  \vx  +  ckvx  +  dkv 

+  fk(x)[6(x  +i)   -   6(x   -i)]    +  Vk  +  Vk 


(2.104) 


with 


k  =  1,    2 

These  two  equations  are  coupled  only  in  the  active  region  II  (here  the 
three  regions  are  still  defined  in  terms  of  the  base-to-base  voltage  variable 
x  =  x   -  x  ).   Except  for  region  II,  each  equation  is  of  the  same  form  as 

(2.103)1 

Thus,  we  can  define  another  plane  where  x  and  x^   are  represented 

independently  but  on  the  same  horizontal  axis.   Call  it  the  x^.  axis. 

In  this  plane,  y  and  y  would  also  be  represented  independently  but 
on  the  same  vertical  axis.   Call  it  the  yfc  axis. 

We  will  still  divide  this  plane  into  three  regions,  but  the  region 
boundaries  will  be  determined  on  the  (x,y)  plane,  rather  than  on  the  (\?Y^) 

plane . 

That  is  to  say:   If  x  is  in  region  I  or  HI  of  the  (x,y)  plane  then 
x  is  in  its  region  I  or  113^  and  x^   is  in  its  region  Ig  or  Illg  of  the 

(*k>yk)  Plane. 

If  x  is  in  region  II  of  the  (x,y)  plane,  then  both  x±   and  x^   will  be 

in  their  respective  regions  II  and  II2  of  the  (xk^yk)  plane° 

Therefore,  region  II,  which  is  nothing  but  the  representation  of  the 
active  region  of  both  planes,  in  the  case  of  the  (xfe,yk)  plane,  will  correspond 
to  two  regions,  one  for  x±   and  another  for  x^      These  regions  are  determined  by 
the  values  of  x±   and  x^   when  \x±   -  xg |  =  -  (see  Fig.  25 ). 
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In  turn,  these  values  depend  only  on  how  they  start,  i.e.,  the 
relative  values  of  their  respective  initial  ordinates  y   and  y  -,  and  which 
one  starts  first  (receives  a  trigger  first).   So,  the  regions  for  x1   may  not 
coincide  at  all  with  the  regions  for  x  ,  and  besides  they  have  a  certain  con- 
figuration only  for  a  given  transition:   i.e.,  in  the  (xk>yk)  plane  the  region 
configuration  is  a  function  of  the  system  and  of  the  triggers. 

Another  plane  is  very  helpful,  and  can  be  used.   It  is  the  (x  ,xg) 
plane,  in  which  the  active  region  is  a  strip  of  parallel  lines  going  through 
the  origin,  intersecting  the  coordinate  axes  at  points  (+  —,  +  -)  thus  bisecting 
the  first  and  third  quadrants.   The  representative  point  Q  of  the  system  is 
the  point  of  coordinates  (x  ,x  ),  and  it  is  a  simple  matter  to  go  from  the  time 
scaled  trajectories  of  the  two  points  P  and  Pg  in  the  (*k>yk)  plane  to  the 
trajectory  of  Q  in  the  (x^x^)   plane. 

Use  of  the  (x  ,x  )  plane  makes  it  easier  for  us  to  find  the  points 
(xla,yla)  and  (x2-,  yg-)  where  (x1&  -  x^-)   =  +  |  (the  sign  +  according  to  the 
direction  of  the  transition),  i.e.,  the  points  where  x  enters  or  leaves  the 
active  region. 

Now,  inside  the  active  region,  equations  (2,10*0  form  a  system  of 
two  linear  second  order  differential  equations  in  x  (t),  k  =  1,  2 .   We  can 
easily  solve  this  system  of  equations  for  x  (t )  and  xg(t),  y^t)  and  yg(t),  and 
so,  x(t)  =  x  -  x  and  y(t)  =  y^  -   y2  can  be  found,  and  from  these,  the  points 
(xlc,ylc),  (x2^,y2^)  where  x  comes  out  of  the  active  region. 

From  then  on,  the  equations  (2.10*0  are  again  independent,  and  the 
remaining  trajectories  y^x  )  and  y2(xg)  can  be  found.   Figure  25  illustrates 
this  discussion. 

The  case  where  both  Cik  and  CQk  or  just  CQk  are  negligible  has  already 
deserved  special  mention  in  the  previous  section,  for  it  is  exactly  solved  by 
equations  (5»l)  to  (5.8). 
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5.5  Other  Types  of  Trigger 
5.5«1  Introduction 

We  have  concentrated  our  efforts  on  a  theory  using  a  rectangular  trigger 
for  two  main  reasons:  the  wave  form  can  often  be  approximated  t>y  a  rectangular 
form,  and  a  rectangular  trigger  lends  itself  easily  to  a  phase  plane  treatment. 

We  feel,  however,  that  some  comments  are  necessary  on  the  most 
common  nonrectangular  trigger  waveforms,  such  as  those  mentioned  in  2.k. 

5°5°2   Impulse  Trigger 

A  trigger  can  be  considered  as  an  impulse  if  the  two  approximate 
conditions  hold: 


(i)  W   »  W  . 

av     mm 


(ii)   q^  «  Aqi  +  £q 


f  (5.9) 


where 


(i)  W    =  average  trigger  current  variable 

cLV 

Wmin  =  minimum  rectangular  trigger  amplitude  necessary 
for  a  transition 
(iii)   q.  =  W    '  T  ,  is  the  charge  transported  by  the 
trigger 
(iv)  T,  =  trigger  duration,  assumed  here  to  be  well  defined 
(v)  &l±,    &lB,    as  defined  in  (4.6,2),  are  the  total 
variations  of  charge  between  the  two  states,  of, 
respectively,  the  input  capacitors  C,  and  the 

IK 

base  storages  <> 
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(vi)   To  simplify  matters,  we  shall  reason  only  with  the 
asymmetric  or  symmetric  flipflops  in  this  section. 
If  conditions  (5.9)  are  met,  there  are  two  ways  to  compute  transition 
times  (in  the  case  of  an  impulse  trigger  the  transition  waveforms  are  meaning- 
less); we  will  assume  that  T  =  TTR. 

As  the  crudest  possible  transition  time  evaluation,  assume  a 
rectangular  trigger  of  amplitude  Wav  and  duration  Tt  =  TTR.   Then 

2  w   =§PIW  (5.10) 

J0       a  av   T  .   av 
oi 


rXIH0  ^  _  a(xIIIQ  -  xIQ) 

TR  "./       V~        mW 

±n  d  J0  av 

XI0 


(5.11) 


We  find  T'   from  its  normal  form  TTR: 


HT  . 
T.   =  2i_ if    T   =  Tmp  (5-12) 

TR    (1  +  RGJ   '         t    TR 


s  0  av 


As  a  less  crude  method,  assume  the  transition  is  complete  when  the 
charge  fed  by  the  trigger  into  the  input  capacitances  and  the  bases  is  equal 
to  the  total  charge  variation  between  the  stable  states: 

^  =  AqB  +  Aq1  (5.13) 

We  are  implicitly  assuming  that  the  charge  lost  through  both  recom- 
bination inside  the  bases  and  the  input  resistances  during  the  transition  is 
small  compared  to  the  variation  of  stored  charge.   The  trigger  duration  is  again 
assumed  to  be  optimum. 
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From   (k.6k), 

AqB  =   HTaIE  (5.14) 

^i  =  ^V1  +  RsG0}  (5.15) 

qt  =aIEWavT;   =aIEWavTiR  <5.1*) 

From   (5.^3),    after  denormalizing  T        into  T'    , 

TR  TR 

Tm  "  vT"  [T  +  V1  +  R.°o>]'      lf     Tt  =  ttr  (5-"' 

av 


where  H  =  1,  2,    is  the  symmetry  factor. 

Further  simplification  in  (5-12)  and  (5-17)  is  possible  if  G  =  0. 


We  get  from  (5.12) 


HT  . 

*JR  -  -^  (5-18) 

av 


and  from  (5.17) 


TiR  "  v5"  [T  +  Tol]  (5.19) 

av 


And  we  see  that  assumption  of  a  constant  value  y  of  y(x)  is  equiva- 
lent to  neglecting  the  transistor's  collector  time  constant  T  with  respect  to 
T   ,  which  may  be  warranted  or  not.   From  (5.19),  we  conclude  that 


TV  =  f-  aL  (5.20) 

mm    IE  \s         * 


is  the  minimum  transition  time  that  can  be  obtained  from  the  given  transistors 
and  trigger  (by  making  C.  =  0).   As  a  result,  we  can  use 
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qB  -  TaL  (5.2D 


as  a  figure  of  merit  of  a  transistor  for  use  in  switching  circuits . 
And  for  the  flipflop  and  trigger  we  have: 

IT'.   =  Hq^  (5-22) 

t  mm     B 


5.5.3  Exponential  or  Sinusoidal  Triggers 

Here  the  trigger  waveforms  are  continuously  changing  functions  of 
time,  and  therefore  phase  plane  treatment  is  not  indicated,  since  the  time 
variable  appears  explicitly  in  the  differential  equation(s)  and  cannot  he 

eliminated. 

We  have  to  work  either  in  the  time  domain  or  in  the  frequency  domain 
by  means  of  integral  transforms.   It  is,  in  general,  easy  to  solve,  directly  or, 
for  instance,  by  Laplace  transforms,  the  three  region  second  order  linear 
differential  equation  under  an  exponential  or  sinusoidal  forcing  function,  so 
that,  in  any  given  problem,  a  numerical  solution  can  always  be  found  for  wave- 
forms, transition  times,  etc. 

A  theory  covering  these  and  other  time- varying  trigger  waveforms, 
i.e.,  finding  analytical  expressions,  relationships,  approximate  formulae  and 
methods  for  the  fast  calculation  of  transition  times,  waveforms,  etc.,  would  be 
an  entirely  new  proposition  altogether,  and  clearly  outside  the  scope  of  a 
phase  plane  theory  of  flipflops  such  as  the  present  work  proposes  to  be. 

5 .6  Use  of  Integral  Transformations 

In  any  of  the  three  regions,  (2.103)  is  a  linear  second  order  dif- 
ferential equation,  and  (2.10U)  is  a  system  of  two  linear  second  order 
differental  equations. 
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Therefore  integral  transform- -or  operational  methods --except  for 
other  reasons --can  be  used,  with  whatever  advantage  one  might  have  from  them. 
In  particular,  Laplace  transform  methods  could  be  used.  The  main  advantage  of 
these  transform  methods  is  that  they  simplify  the  solution  of  the  differential 
equation  under  an  arbitrary  transformable  forcing  function,  in  our  case,  the 
arbitrary  trigger „   One  added  advantage  of  these  methods  is  that  they  make  it 
easy  to  solve  (2.104)  for  x,  which  is  the  system's  state  variable. 

The  results  are  presented  below,  for  completeness : 

From  (2.103)  one  gets: 


Wt(a)  -*[s(t)] 


with  J     x(o)  «*[x(t)] 


a  =  a  +  j  a,  is  a  complex  variable 


X(a)  - 


mya  +  na  +  dv 
*(ava2  +  \a  +  c  ) 


•  Wt(a) 


X0a  +  y0 


+  — 


(5=23) 


where  P0:(xQ,y0)  is  the  initial  point  under  each  Vu-condition, 
From  (2.104)  one  gets: 


with 


K* 


(a)  =*[s.(t)] 


Xk(a)  =^[xk(t)] 


And  also: 


X  =  Xl  "  X2 


X0  '~~~   X10  "  X20;  y0  =:  yl0  '  y20 


x  =  x±  -  x2 


Parameters  are  as  given  in  Table  1.2. 
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(mlva2+n1a+dlv)(a2g2+b2ua+c^)Wtl(  a)    -    (^a^n^^K^+^a+c^)^ 
a{(a.a2+b.    a+c      )(a   CT2+b      a+c2    )   -    [  (*^-*^)o2  +  (^v"^V)a  +  (clv-C2V)]} 


llu  ™luu^lu'^2u  '~2uu'~2u 


2 
a 


where  P  : (x  ,yn)  is  the  initial  point  under  each  Vu  condition. 

This  makes  it  obvious  why  the  approximate  and  graphical  methods  are 

important! 

We  can  also  find  X-^a)  and  X2(a),  hy: 

K°2  +  Kva  +  ckv}   '   (a2x(a)  '  xoa  -  yo}  +  a  '    Kva2  +  \a+^vh\ia) 

Xv(o)  = 2~ 2 


k' 


a  (mkva  +  V  +   d^) 


(5-25) 


XkOJ  +  yk0 

2 
a 


where  P   • (x   ,y  n)    is  the  initial  point  under  each  Vu  condition. 
™     kO  v  kO'^kCr 

Since  X(a)  would  have  to  be  found  first,  this  makes  it  doubly  obvious 
why  the  approximate  and  graphical  methods  are  important. 

5.7  Summary 

In  this  chapter  we  have  shown  how  two  degenerate  cases  (t  =  0  and 

C   =  C   =0)  relate  to  the  theory  presented  so  far.   The  case  of  1   =   0  was 
ik    ok 

shown  to  be  essentially  included  in  the  theory,  since  T  has  been  used  as  a 
normalization  constant  for  no  other  reason  than  that  of  convenience.   The  other 
(,flqp  0f  r       =   C       =   0.  or  C   =0,  have  been  shown  to  be  exactly  solvable,  the 
first  even  for  the  nonsymmetric  flipflop,  and  the  second  for  at  least  the 
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asymmetric  and  symmetric  flipflops,  and  possibly  (if  R   =  R   )  for  the 

si    s2 

nonsymmetric  case  which,  at  any  rate,  is  reduced  to  a  system  of  the  second 
order. 

A  phase  plane  method  suitable  for  the  general  nonsymmetric  Eccles- 
Jordan  flipflop  was  given,  and  a  discussion  showed  that  there  are  areas  where 
the  nonsymmetric  flipflop  is  equiA^alent  to  two  independent  asymmetric  flipflops 
(regions  I  and  III);  the  trajectories  in  the  active  region  (region  III)  must  be 
found  by  solution  of  the  system  in  the  time  domain. 

Finally,  we  have  discussed  other  types  of  trigger  waveform.   Besides 
the  almost  trivial  impulse  trigger,  for  which  some  relationships  have  been 
established,  the  other  cases,  such  as  exponential  or  sinusoid,  cannot  be  treated 
by  a  phase  plane  theory.   They  can  be  treated  analytically  or  numerically; 
however  no  general  results  are  available.   The  equations  must  be  solved  in  each 
specific  case. 

Approximations  (2.103)  and  (2.104)  are  also  important  in  allowing  a 
phase  plane  treatment  of  the  most  important  cases  (second  order),  besides 
allowing  treatment  in  the  time  domain  (directly)  or  in  the  frequency  domain 
(integral  transforms)  for  any  case. 

Finally,  we  have  briefly  discussed  the  application  of  Laplace  trans- 
form methods  to  (2.103 )  and  (2.104),  and  presented  special  formula  (5.23)  and 
entirely  general  formula  (3.2k),    thus  covering  all  possibilities. 


6.   EXPERIMENTAL  EXAMPLES 

6.1  Introduction 

The  present  chapter  has  a  double  purpose.   We  wish  to  illustrate  the 
application  of  some  of  the  described  procedures,  and  also  to  test  the  accuracy 
of  the  theoretical  results  as  compared  to  experimental  fact.   No  extensive 
program  of  experimentation  is  intended;  only  a  few  examples  were  treated  which 
should  suffice  to  provide  seme  feeling  for  the  quality  of  the  theory. 

The  experiments  we  have  carried  out  consist  in  triggering  a  flipflop 
with  a  rectangular  current  trigger  from  what  was  practically  a  current  source, 
i.e.,  the  collector  of  a  transistor.   The  trigger  had  a  reasonably  good  waveform 
but  we  did  not  attempt  to  obtain  an  exceptionally  good  rectangular  shape. 

As  for  the  flipflops  themselves,  we  took  two  classes:   one  was  a 
slowed-down  flipflop  where  relatively  large  capacitors  were  paralleled  with 
the  T  base-to-ground  and  T  collector-to-ground  terminals;  the  other  had  just 
parasitic  capacitances,  which  were  carefully  measured.   Only  the  asymmetric 
structure  was  used.   In  each  case  transition  times  and  waveforms  were  measured 
and  recorded  for  different  values  of  trigger  amplitude  and  various  values  of  R± 

and  R  . 
o 

Corresponding  calculated  values  were  found  and  comparisons  between 
theoretical  and  experimental  values  are  presented  in  the  tables.   The  transistors 
used  were  the  same  for  all  flipflops,  2N1309's. 

6.2  Measurement  of  T,  C.  and  C 

The  collector  time  constant  T  determines  the  influence  of  the  base 
current  terms  upon  the  solution  of  the  flipflop  equation. 

Whenever  T  is  negligible  compared  to  the  other  time  constants  of  the 
system,  it  becomes  irrelevant  and  the  base  current  terms  may  be  ignored,  as 


■168- 


-169- 

in  Chapter  5.   But  if  T  is  comparable  to  the  other  system  time  constants,  it 
becomes  critical  and  must  be  carefully  measured. 
The  system  used  here  was  as  follows: 

a)  The  transistor  pair  whose  collector  time  constants  (assumed  equal)  are  to  be 
measured  were  assembled  into  a  switching  amplifier,  with  no  collector  lead, 
and  a  current  of  1  ma  fed  into  the  parallel  emitters . 

b)  A  (periodically  repeated)  step  voltage  with  amplitude  just  enough  to  switch 
the  current  from  one  transistor  to  another  was  applied  to  the  base  of  T  , 
and  the  base  current  of  Tg  was  recorded  and  integrated  with  respect  to 
time  (see  Fig.  25). 

c)  Since  there  is  only  a  negligible  voltage  variation  at  the  base  of  T 
(grounded  lead)  the  parasitic  capacitances  have  only  a  negligible  effect 

on  the  measurement.   Recombination  current  can  also  be  neglected  in  compar- 
ison to  the  storage  current. 

From  Fig.  25  we  obtain  by  integration 

T  =  15.1  nsec  (6.1) 

The  parsitic  capacitances  have  to  be  measured  in  situ.    This  can  be 
done  by  measuring  the  time  constants  of  voltage  curves  under  applied  step 
currents.   So,  Figs.  27  and  28  yield  C.  and  C  in  all  regions .'  C.    is  found  to 
vary  slightly  from  one  region  to  another  (Figs.  6.5a  and  6.5b)  but  C  remains 
essentially  the  same  in  all  regions.   In  calculating  C  it  is  necessary  to  sub- 
tract the  injected  current  time  constant  (t   =  15  nsec)  from  the  total  collector 
voltage  time  constant  (t   =  6l  nsec)  in  order  to  obtain  the  true  collector 

circuit  time  constant  (t   =  R  C  =  k6   nsec). 

000  ' 


C^  is  the  base-to-ground  capacitance  of  T^;  C  is  the  collector-to-ground 
capacitance  of  T   (see  Fig.  3  and  also  Fig.  5  for  comparison). 
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FIGURE  28:   COLLECTOR  VOLTAGE  RISE  UNDER  INJECTED  CURRENT 
Curve :  v  vs .  t 

Scales :  Vert . :  90  mv/div 

Horiz.:  20  nsec/div 

Total  time  constant:  ^^  =   6l  nsec 

Time  constant  of  the  injected  current:1  Tj  =  15  nsec 

Time  constant  of  the  collector  circuit: 

r  q  =  t  =46  nsec 

o  o    o 

Since  RQ  =  1  KB,  Cq   =  h6   pf 


.1 


Not  shovn  here. 
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This  is  not  done  in  the  measurement  of  C,  since  in  normal  operation 

the  trigger  circuit  does  contribute  to  C,  whereas  it  does  not  contribute  to 

C  . 
o 

We  get  the  values  presented  in  (6.4). 

By  way  of  approximation  we  have  used  as  a  value  of  C.  in  region  II 
the   average  of  its  values  in  regions  I  and  III,  although  in  reality  it  is  a 
continuously  changing  value. 


6.3  Equation  Parameters 

We  have  considered  two  possibilities: 

1.  A  flipflop  loaded  with  relatively  large  capacitors 

2.  A  flipflop  with  only  parsitic  capacitances. 

For  case  1  we  used 


and  had 


C.  =0.02  uf   and   C  =0.01  uf  (6.2) 


R.  =  2  kft     and     R  =  1  kft  (6.3) 

1  o  \   -»/ 


In  case  2  there  are  only  parasitic  capacitances 


c±-< 


56.7  in  region  I 

6l.7  in  region  II       and      C  =  46.0  pf       (6.4) 

66.7  in  region  III 


The  two  resistors  were  chosen  to  be 


R.  =  2  kfi   and    R  =  1  kfi  (6.5) 

1  o 


Besides  this  we  had 


Case  1:   Ig  =  3  ma;  1^   =  1   ma,    p  =  10 

(6.6) 
Case  2:   I  =  5  ma;  L,  =  0.9  ma,   p  -  10 

with  the  values  of  p  calculated  on  the  assumption  that  —  =  h0   v 

In  both  cases  the  value  of  E  was  adjusted  to  make  the  system  state 

c 

symmetric,  i.e.,  to  have  the  stable  values  of  the  base  voltage  of  T  symmetric 
with  respect  to  ground. 

Table  XI  presents  these  parameter  values  in  a  convenient  form. 

The  equation  coefficients  are  obtained  from  Tables  I  and  II,  and 
presented  in  Table  XII. 

Notice  that  all  three  cases  are  normalized  with  respect  to  T. 

6.k     An  Illustrative  Example 

In  order  to  illustrate  some  of  the  techniques  described  in  the  previous 
chapters,  we  shall  consider  case  2,  under  a  rectangular  current  trigger  of 
amplitude 

W  =  0.667,    i.e.,    I  =  0.6  ma 

and  we  will  calculate  the  delay,  active,  balance  and  transition  times  by 
three  different  methods,  and  compare  the  theoretical  results  with  the  experi- 
mental ones . 

a)  Graphical  Method  A:   see  Figs.  32  and  33 

b)  Approximate  Method  B:   see  Fig.  3^ 

c)  Iterative  Numeric  Method. 
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TABLE  XI.   PARAMETERS  FOR  THE  TWO  EXPERIMENTAL  FLIPFLOPS 


Parameter 

Case  1 

Case  2 

C. 

1 

0.02  uf 

61.7  +  5 

--1- 

0 

L+i. 

pf 

co 

0.01  (if 

46.0     pf 

R. 

1 

2  kfi 

2         kfi 

Ro 

1  kft 

1         kft 

h 

3  ma 

5        ma 

h 

1  ma 

0.9      ma 

Ec 

adjusted  for  state  symmetry 

T. 

1 

40  (isec 

123.4  +  10 

--1- 

0 

u+lJ 

nsec 

T. 

io 

20  usee 

61.7  +  5 

0 
-+1- 

nsec 

T  . 
01 

20  usee 

92.0       nsec 

T 
O 

10  usee 

46.0       nsec 

Tn 

T 

T 

T 

15.1  nsec 

7 

0.7 

0.7 

P 

10 

10 

1  -  a 

0,007 

t  («) 

0.979 

.0-979 

cp  (») 

0 

■ 
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TABLE  XII.   PARAMETERS  AND  CONSTANTS  INVOLVED  IN  THE  EQUATIONS 
REPRESENTING  THE  TWO  EXPERIMENTAL  FLIPFLOPS 


Coefficient 


a^, 


vn 


'vn 


a. 


f(S) 


\o 


*, 


III 


A 


*EI 


\ 


°al 


^1 


\ 


all 


"bit 


A 


"alii 


Bill 


Case  1 


6  4 

1.755    •    10     +  10 


o- 

1 
•0- 


4.63  •  103  +  86 


rO- 
1 
-0 


1 
-6 
iJ 


[1 

L+10-J 


26.5   '  103 


60 


■9.275   ■  10 


3 


15.1 


10 


3 


w 


5.26  •  10"3 


10 


-1 

0 
L+1J 

60  •  w 


-10  ■  w 


-3 


-2.40    •    10 

-0.236  •  10"3 


-3.611  •  10 

+0.913    '    10" 


-3 


-2. 40 


10 


-3 


-0.236   •  10"3 


Case  2 


r-22.9' 

29.0 

L26.9- 


1-16.6 

38.3 

Ll8.o-J 


-10 
0 

+10J 


75.2 
81.8 

488. 4 J 


6 


60 


-28.65   ■  y 


3.28  ■  W 


O.988 


10 


1 

0 
L+iJ 

60    •    W 


-10    •    W 


-0.658 
-O.O67 

-1.46 
+0 .  iks. 

-0.608 
-0.0610 
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a)  Upper  curve:   v.  vs.  t 

Lower  curve :  v  vs .  t 
o 


Scales 


Vert. :      0.25   v/div 
Horiz.:      10  p.sec/div 


b )     v .    vs .    t 


v .    vs .   t 

i 


Long  trigger  Optimum  trigger 


v     vs .    t 
o 


v     vs .    t 
o 


Long  trigger    Optimum  trigger 

Scales:  (  Xert' :  W^/ 

[  Horiz. :   0.5  msec/div 

c)  v.  vs.  v.  (approximate) 

J  Upper  curve:   long  trigger 
[Lower  curve  :  optimum  trigger 

Scales:  j^ert.:  12-5  (v/msec)/div 

'  Horiz.:   0.25  v/div 


FIGURE  29:   TRANSITION  CURVES  FOR  CASE  1 


Comments : 

Part  a)  illustrates  the  relationship  between  the  transition  curves  for  T 
base  voltage  v.  and  T  collector  voltage  v  .   Delay,  active  and 
balance  times  are  apparent.   The  time  duration  of  the  transmission 
phases  can  be  measured  from  it. 

Part  b)  illustrates  the  effects  of  trigger  duration  upon  the  waveforms  of 

both  v.  and  v 

1      o. 

Part  c)  is  an  approximate  phase  plane  portrait  of  the  transition,  both  for  the 
long  trigger  (upper  curve,  showing  the  return  of  P  from  X   to  X  , 
after  trigger  turn-off)  and  optimum  trigger.   Both  curves  are  slightly 
tilted  to  the  right  due  to  imperfect  differentiation  of  v.. 
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FIGURE  30:   TRANSITION  CURVE  FOR  CASE  2, 


WITH  W  ==  0.^5;    v.    vs.    t 


Scales 


Vert. :      0.22?  v/div 
Horiz. :      20  nsec/div 
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FIGURE  31:   TRANSITION  CURVE  FOR  CASE  2, 
WITH  W  =  0.667;  v.  vs.  t 


Scales 


Vert.  : 
Horiz. 


0.225  v/div 
20  nsec/div 
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6.4.1  Graphical  Method  A 

We  calculate  ^  x^,  Ax^,  Ax^,  x^,  x^  and  yQ  by  the 
formulae  on  Table  II: 


xio  =  -vf~-v  =  -10 


Ax_  =  nW  =  40.0 


XI1  =  XI0  +  ^1  =  3°"° 


XIII0  =  +P  f^  *  +P  =  +10 


Ax    =  nW  =  40.0 


XIII1  =  XIII0  +  ^III  =  50-° 


mi 
y0=-W=2.19 


We  also  need  I   and  ^STI;  we  find  from  the  note  on  equations  (3.60). 


m 

J  =  -S-   =  O.988 
da 


and,  obviously, 


-bn  +  ^11  -  Sicn     ni)c 

Vl  =  2^ =  °'lk2 


Then  (see  Fig.  34) 
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(i)  We  draw  the  line  P  Xx, 
v  '  o  II 

(ii)  Intersection  of  P  XT,  with  x  =  -  -  is  P 
v  '  o  II  y  a 


1 

XH  +  7   ..   ,  no 


(iii)  Using  equation  (3.6la.i)  we  find  y^  from  yg 
We  have: 


\ -  vt* 


and: 


I   =  0.988 


so: 


yb  =  0.92 


(iv)  Draw,  through  P  ,  the  line  of  slope  kmjj*   whose 

intersection  with  x  =  +  —  is  P  .   Since 

/     c 


xall  "  °-lte 


through 


Ax  =  2.86 


causes 


Ay  =  0.1+06 


-185 


so  that 


yc  =  i-33 


(v)  P  is  found  from  P  by  (3„6lb.i) 


c 


yd  -  (i  +  %c)yc 


We  have 


yc  =  1.33 


and: 


I  =  O.988 


so: 


ya  -  3.oo 

(vi)   (Simplified  in  our  case)  we  draw  the  line  PnX^^Tn 

d  III1 

and  its  intersection  with  x  =  XTTT„  is  P  . 

IIIO     e 

ye  =  yd  • —  =  2o52 

XIII1  "  7 

The  time  intervals  over  the  trajectorycan  be  calculated  by  using 
equation  (4„2l),  which  we  repeat  here,  for  convenience,  in  a  slightly  different, 
but  equivalent  form; 


b    a   Ay     y  v  •  ■  / 

cl 
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Applied  to  region  I  we  get: 


AXT      ya 

T  =  — *  fa    -2   =  1+.31 

d   y0    y0 


1       yc 

A    Vl      yd 


Aacj.^. 


Therefore^ 


-iS  ^  -S  =  3.05 
ye      ^d 


TT  =  TD  +  TA  +  TB  =  9.91 


rp>  =  t   •  T  =  65.1  nsec 
D    D 


mt  _  m   •  t  =  38. 4  nsec 
A    A 


Ti  _  T   .  1   -     46,0  nsec 
B    B 


T»  =  t   0  t  =  1^9.5  nsec 


Similarly  we  can  find  the  trajectory  for  the  other  cases  and  fill 

out  Table  XIII. 1. 

Experimental  values  for  the  time  durations  over  the  trajectories  are 
taken  from  Figs.  29  and  30,  which  are  reproductions  of  pictures  of  oscilloscope 
images  of  the  actual  transitions. 

Observe  on  Table  XIV. 1  the  excellent  agreement  between  experimental 
and  calculated  values  of  total  transition  times  (TT).   Notice  also  that  in 
case  1,  the  agreement  for  the  partial  time  durations  (TD,  Ta,  Tg)  is  also 
excellent.   But  in  case  2  the  agreement  for  the  partial  time  durations,  although 
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fair,  is  less  satisfactory .   There  is,  however,  a  compensation  of  errors, 
especially  between  values  of  TA  and  T  .   This  seems  to  be  a  property  of  this 
approximate  way  of  taking  the  base  current  into  account .   When  the  base  current 
has  negligible  effect,  as  in  case  1,  this  effect  does  not  occur. 

An  apparently  paradoxical  fact  is  that  we  get  a  better  agreement  on 
partial  times,  although  worse  for  the  total  transition  time,  if  we  take  i  =   0, 
i.e.,  if  we  completely  neglect  the  effect  of  the  second  derivative  of  the  base 
current.   The  graph  is  shown  in  Fig.  33,  and  Tables  XIII. 2  and  XIV. 2  present 
the  ordinates  and  times  for  this  very  simple  procedure. 

This  apparent  paradox  can  be  explained  if  we  consider  our  type  of 
approximation  at  the  model  level.   The  piecewise  linear  model  has  a  behavior 
that  differs  in  detail  from  the  nonlinear  real  transistor  pair;  but  the  dif- 
ference is  such  that  it  tends  to  cancel  out  over  the  different  phases  of  a 
trajectory,  yielding  good  overall  results. 

If  an  approximation  at  the  equation  level  had  been  used,  no  impulses 
would  appear.    Clearly,  ignoring  the  impulses  brings  us  "closer"  to  such  a 
type  of  approximation  in  the  sense  of  yielding  a  solution  differing  less  from 
the  exact  one.   The  differences,  however,  although  smaller,  do  not  average  out, 
and  the  overall  result  is,  as  expected,  not  so  good. 

Considering  the  crudeness  of  this  graphical  method  B,  it  is  an 
additional  advantage  that  the  results  tend  to  be  conservative,  yet  fairly  close 
to  the  experimental  values . 

A  last  comment  is  necessary  here.   We  pointed  out  the  fact  that  in 
case  1,  where  we  had  the  (in  practice)  exact  values  of  all  the  circuit  components^ 


This  would  also  be  true  if  the  approximation  at  the  model  level  had  a  smooth 
characteristic  (i.e.,  continuous  with  different iable  first  derivative),  which 
is  not  the  case  of  the  piecewise  linear  function  <f(x). 
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including  known  constant  capacitors  the  agreement  is  as  good  as  one  can 
possibly  have.   So,  ve  infer  that  part  of  the  error  obtained  in  case  2  should 
be  blamed  on  nonconstant  capacitances  and  the  necessarily  incorrect  values 
taken  for  them. 

6.4.2  Approximate  Graphical  Method  B 

Figure  3^  and  the  description  given  in  subsection  k.k-,3   entirely 
complement  each  other. 


w 


x    =  &  (x  )  =  J2W  _  *=  £  _  ^  =  0.07^5 

0   dx  v  I0;   ayQ   a   a   a 


ith  a,  b,  m,  n  calculated  for  V$A  »  II.   And 

x    .  cjZ  /+  1)  „  ^-z^il  .  3!  „  0.0757 
d   dx  v   7      ayd     a 


where  v  =  2,52  is  obtained  from  the  graphical  construction  itself. 
J  & 

\    and  KaTTT1    are  given  in  Table  XfL   Observe  that,  in  order  to 
(311      pIIH 

obtain  a  closer  detailed  approximation,  we  have  taken  I  *  0. 

The  numerical  values  of  the  ordinates  of  the  break  points  are  pre- 
sented together  with  the  graph  in  Fig.  3^  itself. 

Successive  application  of  equation  (4.21 )  to  the  five  branches  a,  b, 
C,  d,  e,  of  the  piecewise  linear  approximate  phase  plane  path  yields  the  time 
intervals  over  these  five  phases  of  the  transition. 

This  is  certainly  not  true,  and  is  actually  a  somewhat  crude  approxi- 
mation, consistent  with  good  results  for  two  reasons:   region  II  is  relatively 
narrow,  so  this  assumption  affects  little  the  value  of  TA,  and  the  error 
propagates  only  to  region  III,  yielding  poor  values  of  Tg,  but  still  good  values 
of  Tm. 
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The  following  variation  of  method  B  is  less  subject  to  this 
restriction: 

a)  Instead  of  assuming  the  path  to  be  a  line  of  slope  LTT1  inside  region  II. 

plxl 

we  draw  the  lines  of  slope  A,  x.rn  and  A.OTT.n  through  point  XTTn  ,  and  we  draw, 

alii      pill      ^  ^      III'  J 

through  P  ,  an  approximate  hyperbolic  arc  consistent  with  the  two  asymptotes 

This  arc  shall  be  taken  as  an  approximation  to  the  trajectory  inside  region 

II.  and  its  intersection  with  the  line  x  =  +   —  determines  the  ordinate  y, 

7  d 

of  P,  (we  take  ,0  =  0). 
d  v 

b)  We  proceed  in  region  III  as  for  the  regular  method  Bo 

The  results  are  presented  in  Table  XVI,,  along  with  corresponding 
experimental  values .   Notice  that  this  method  allows  us  to  distinguish  more 
distinct  phases  of  the  path,  i.e.,  one  phase  for  each  value  of  ~   (over  this 
approximate  broken  line  path). 

We  can  see  that  the  agreement  is  exceedingly  good  for  region  I  and  II, 
even  for  the  detailed  shape  of  the  curve,  but  it  is  not  so  good  in  region  III. 

The  difficulty  is  mostly  due  to  the  calculation  of  the  branch  in  the 
active  region.   Any  small  error  in  the  calculation  of  point  y,  is  magnified 
throughout  region  III,  so,  even  with  a  numerical  method  we  should  expect  larger 
errors  in  region  III.   Even  the  variations  of  temperature  would  cause  errors 
by  changing  the  transistor-pair  characteristic  through  ^  affecting  especially 
the  boundaries  of  region  II. 

So  far  we  have  assumed  that  the  trajectory  is  a  straight  line  of 
slope  ^RTT-|  throughout  region  II. 

c)  Guided  by  the  straight  lines  in  region  I  and  III,  and  by  the  points  P  and 

3, 

P.,  we  draw  the  parabolic  arcs  in  region  I  and  III,  and  join  them  smoothly 
by  a  suitable  arc  inside  region  II,  so  that  these  three  arcs  form  a  curve 
which  shall  be  taken  to  be  an  approximation  to  the  real  trajectory. 
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TABLE  XV o   COMPARISON  OF  TIME  INTERVALS  OBTAINED  BY  METHOD  B 
WITH  EXPERIMENTAL  RESULTS 


Time  Interval  (in  nsec) 

Over  Branch 

Theoretical 
(by  Method  B) 

Experimental 

23°5 

22 

a 

32.6 

30 

b 

56.1 

52 

a  and  b,  i.e., 
region  I 

19.0 

22 

c,  i.e., 
region  II 

18.1 

26 

d 

30.2 

38 

e 

1+8.3 

6k 

d  and  e,  i.e.,, 
region  III 

124.3 

138 

Total  Transition 
Time 
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TABLE  XVI o   COMPARISON  OF  TIME  INTERVALS  OBTAINED  BY  A  VARIANT 
OF  METHOD  B  WITH  EXPERIMENTAL  RESULTS 


Time  Interval  (in  nsec) 

Over  Branch 

Theoretical 
(by  a  Variant 
of  Method  b) 

Experimental 

20.1 

22 

a 

35«2 

30 

b 

55*3 

52 

a  and  b,  i,e.; 
region  I 

12.0  +  8o9 
=  20„9 

22 

c  and  d,  i.e., 
region  II 

23.3 

26 

e 

30.6 

38 

f 

53.9 

6k 

e  and  f,    i„e„^ 
region  III 

130.1 

138 

Total  Transition 
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d)  Now  we  draw  a  convenient  piecewise  linear  approximation  to  this  curve  and, 

by  repeated  applications  of  equation  (^.21 )  to  its  linear  branches,  we  can 

\ 

calculate  the  corresponding  time  intervals .' 

Note,  in  Table  XVI,  that  the  tendency  persists  for  longer  theoretical 
time  intervals  in  region  I,  and  shorter  in  region  III,  but  results  are  somewhat 
better  than  those  for  plain  method  B.   Note  also  that,  since  this  method  makes 
use  of  a  degree  of  arbitrariness,  results  are  bound  to  vary  a  little  with  the 
operator's  judgment. 

6.U.3   Iterative  Numerical  Method 

The  iterative  numerical  method  offers  no  special  difficulty.   It  is 
clear  that  equation  (k.kb)   converges  in  every  region.1'   The  values  of  M^,  Q^, 
M  ,  Q  as  well  as  the  computed  ordinates  yQ,  y  ,  yd,  and  yg  are  shown  in 
Table  XVII.   We  have  again  taken  I   =  0. 

A  direct  application  of  equation  (3-58)  for  each  region  yields  the 
respective  time  durations,  which  are  presented  in  Table  XVIII. 

Notice  that  the  results  are  somewhat  better  than  for  the  other  methods, 
but  they  are  not  perfect.   It  is  even  worse  for  region  II.   However,  this 
numerical  method  is  the  exact  solution  of  the  differential  equation,  s o  we 
conclude  that  an  error  is  due : 


We  might  be  tempted  to  use  equation  (3-58)  but  this  could  be  disastrous. 
This  equation  is  valid  only  for  points  exactly  on  a  true  trajectory,  and  is 
very  critical.   If  the  points  are  only  approximately  on  a  trajectory,  the 
results  obtained  by  use  of  (3-58)  will  probably  be  meaningless. 

A. 

' '  A  rule  of  thumb  is  to  work  always  with  the  large  magnitude  exponent  —  . 
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TABLE  XVII.   PARAMETERS  AND  TRAJECTORY  KEY  ORDINATES 
FOR  THE  ITERATIVE  NUMERICAL  METHOD 


Region  I 

Region  II 

Region  III 

M 
a 

20,68 

-11,82 

24,32 

% 

26,32  -  yQ  =  24,13 

-7»65  -  yfl  =  -9.716 

29.^5  -  yd  =  28,24 

% 

2,105 

1.15 

2,44 

% 

2,68  -  yQ  =  0,49 

0.7^  -  y   -  -1,322 

cl 

2.955  -  yd  =  1.7^1 

70  •   2-19 

=>  y  =2,066 

El 

=>  yd  =  1.214 

=>  ye  =  2,294 

TABLE  XVIII,   COMPARISON  OF  TIME  INTERVALS  OBTAINED  BY  THE  ITERATIVE 
NUMERICAL  METHOD  WITH  EXPERIMENTAL  RESULTS 


Time  Interval  (in  nsec) 

Region 

Theoretical 

Experimental 

58,4 

52 

I 

31  =  3 

22 

II 

61,6 

64 

III 

151  =  3 

138 

Total 
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a)  to  the  piecewise  linear  approximation  use  for  the  most  exact  nonlinear 

differential  equation . 
t>)  to  assumptions  of  lumped  parameters,  constant  in  each  region. 

c)  to  the  measurement  of  circuit  parameters  and  transistor  constants . 

d)  possible  departure  of  the  transistor  characteristics  from  the  ideal  one 
we  have  assumed. 

e)  the  imperfection  of  the  rectangular  trigger  used  in  the  experiments. 

Under  the  above  considerations  the  error  of  about  ten  per  cent  in 
total  transition  times,  along  with  the  good  agreement  in  waveform  (in  the  case 
of  method  B,  for  example)  is  a  satisfactory  result. 


7°   CONCLUDING  REMARKS 

7.1  Summary 

The  purpose  of  this  investigation  was  to  describe  in  detail  the 
operation  of  flipflops  from  a  mathematical  point  of  view,  and  to  devise,  based 
on  this  mathematical  description,  practical  methods  of  analysis,  design  and 
optimization  of  both  flipflop  and  triggering  circuits. 

The  mathematical  description  has  been  accomplished  with  the  establish- 
ment of  equations  (2 087),  (2.88)  and  (2 .90)  in  Chapter  2,  and  with  those 
qualitative  aspects  of  their  piecewise  linear  approximations --equations  (2.99), 
(2.100),  (2 „ 102) --which  clearly  apply  to  the  original  system- 
Methods  of  analysis  and  design  were  devised  by  means  of  a  detailed 
study  on  the  phase  plane  of  the  piecewise  linear  equations,  taken  as  approxima- 
tions to  the  original  nonlinear  equations.   The  singularities  of  the  system, 
the   conditions  for  their  existence  and  the  dependence  of  their  nature  upon 
the  system  parameters,  have  been  thoroughly  described.   The  phase  plane  portrait 
of  the  system  was  described  with  some  emphasis  on  separatrices,  trajectories, 
and  the  influence  of  the  singularity  corresponding  to  a  given  region,  whether 
this  singularity  exists  in  its  proper  region  or  has  a  virtual  image  in  another 
region. 

Based  on  this  study  some  engineering  methods  of  analysis  and  design 
have  been  described  in  Chapter  4,  and  some  simplified  formulae  for  the  rapid 
estimate  of  flipflop  behavior  have  been  presented  in  Chapter  5. 

The  experimental  example  presented  in  Chapter  6  illustrates  the  use 
of  some  of  these  methods,  and  also,  by  comparing  theoretical  with  practical 
results,  some  feeling  is  obtained  for  the  adequacy  of  the  various  methods  and 
for  the  type  of  approximations  (piecewise  linear  at  the  model  level)  used  in  the 
theory. 
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7 ,2   Conclusions 

It  is  apparent  that  we  have  obtained  a  useful  and,,  for  most  practical 

purposes ,  adequate  theory . 

We  feel,  however,  that  there  are  some  questions  to  which  we  do  not 
have  even  unsatisfactory  answers.   A  first  question  is:   why  is  it  that  the 
test  of  all  methods  when  applied  to  the  active  region  yields  a  path  which 
obviously  differs  considerably  from  the  true  path?  Even  the  crude  device  of 
assuming  the  path,  in  region  II,  to  be  a  constant  equal  to  y&  would  produce  a 
result  closer  to  the  true  one  in  that  region! 

Another  question  is :  why  is  it  that  results  are  worse  if  the  impulses 
(second  derivatives  of  cp(x))  are  considered  than  the  results  we  get  when  they 
are  ignored? 

We  feel  that  the  answer  lies  in  a  more  detailed  study  of  the  relation- 
ship between  a  nonlinear  differential  equation  (especially  of  the  second  order! ) 
and  another  equation  which  formally  is  a  piecewise  linear  approximation  to  the 
nonlinear  one„   Specifically,  what  are  the  effects  of 

a)  the  break  points  (error  in  derivatives!) 

b)  the  error  itself 

c )  the  constancy  of  coefficients 

on  the  solution  of  the  approximate  equation  with  respect  to  the  original  one? 
The  present  investigation  gives  the  impression  that  this  type  of  approximation 
should  be  studied  in  detail  and  formalized . 

7„3  Further  Investigations 

There  are  three  directions  for  further  investigation: 
a)  The  study  of  approximate  solutions  to  nonlinear  differential  equations  by 
use  of  solvable  formally  approximate  equations  to  the  original  one,  such  as 
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piecewise  linear  equations  or  other  standard  types  of  equations  with  known 
solutions  o 

b)  The  polishing  of  the  present  theory  by  considering  other  types  of  approxi- 
mation, such  as,,  for  example,  approximation  at  the  equation  level. 

c)  Application  of  the  ideas  we  have  described  to  more  complex  situations,  for 
example, 

(i)   considering  the  nonlinearity  of  parasitic  capacitances, 
(ii)  taking  the  collector-base  junction  capacitances  into 
account, 
(iii)  considering  inductances  in  the  passive  circuit, 
(iv)   considering  the  distributed  nature  of  some  of  the 
parasitic  capacitances. 
Advances  in  one,  some,  or  all  of  these  directions  would  certainly  improve  the 
present-day  techniques  of  switching  circuit  design  for  digital  computers. 
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